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ENERGY  BASIS  OF  DISASTERS  AND  THE 
CYCLES  OF  ORDER  AND  DISORDER 

By 

John  Franklin  Alexander,  Jr. 

August  1978 

Chairman:  Howard  T.  Odum 

Major  Department:   Environmental  Engineering  Science 

A  quantitative  theory  of  cycles  order  and  disorder  was  applied  to 
the  earth  and  evaluated  to  form  an  energy  basis  for  the  global  cycles, 
surges,  and  disasters.  Energy  circuit  language  was  used  to  diagram  the 
world  system  and  show  a  common  pattern  in  the  order-disorder  processes. 

Storms,  floods,  forest  fires,  volcanic  eruptions,  earthquakes,  ur- 
ban fires,  and  wars  were  modeled  as  the  catastrophic  release  of  energy 
previously  converged  and  stored.  Released  energy  disordered  and  re- 
cycled material  available  to  stimulate  a  new  cycle  of  growth.  Cascading 
of  catastrophic  processes  of  disasters  was  modeled  with  a  world  web. 
The  feedback  in  the  global  energy  web  was  provided  by  the  control  action 
of  disaster  pulses.  The  global  model  was  presented  in  both  diagrammatic 
and  differential  equation  form  with  the  energy  flows  and  storages  eval- 
uated. Order-disorder  models  of  the  atmospheric,  oceanic,  biological, 
geological,  and  urban  systems  of  earth  were  connected  to  form  an  energy 
convergence  network. 

The  global  energy  model  was  used  to  calculate  energy  quality 
factors  (ratio  of  energy  of  one  type  generating  energy  of  another  type) 
for  the  earth's  major  energy  transformations.  The  ratio  of  solar  energy 
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required  to  produce  a  catastrophic  event  to  the  energy  released  in  the 
event  (Calories  per  Calorie)  was  found  to  be  1  x  107:1  for  volcanoes, 
5  x  105:1  for  earthquakes,  1  x  105:1  for  urban  fires,  floods,  and  wars, 
2  x  lO^rl  for  forest  fires,  and  2  x  103:1  for  storms.  The  energy  qual- 
ity factors  of  stored  energy  were  used  to  calculate  embodied  solar  energy. 

The  hypothesis  of  selection  for  maximum  power  was  used  to  explain 
the  prevalence  of  systems  that  recycle  structure  by  the  catastrophic 
release  of  energy  pulses  generated  in  disasters.  The  solar  energy  em- 
bodied in  the  disordered  city  structure  was  found  to  be  1.4,  1.3,  and 
1.7  times  the  solar  energy  embodied  in  storm,  flood,  and  seismic  pulses 
of  the  global  web. 

The  dynamic  properties  of  several  configurations  of  models  of  the 
cycles  order-disorder  were  analyzed  by  solving  the  model  in  differential 
equation  form  and  through  computer  simulation.  The  most  suitable  model 
has  a  production  function  supplying  energy  to  a  consumer  that  alternately 
shifts  from  linear  flow  to  a  surge  of  exponential  growth  causing  an 
oscillation. 

Suitability  of  order-disorder  models  were  further  tested  by  appli- 
cation to  three  case  studies:  (a)  effect  of  gamma  radiation  stress  on 
a  tropical  rain  forest  in  Puerto  Rico;  (b)  the  disordering  of  Johnstown, 
Pennsylvania,  by  the  great  flood  of  1889;  and  (c)  the  destruction  of  El 
Prcgreso,  Guatemala,  by  the  earthquake  of  1976.  The  dynamic  effects  of 
the  flood  on  Johnstown  and  the  earthquake  on  El  Progreso  were  modeled 
with  two  order-disorder  models  in  cascade  and  simulated.  The  resulting 
graphs  of  the  destruction  of  the  towns  studied  compared  favorably  with 
observed  results. 
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The  energy  flows  that  produced  order  in  the  radiation,  flood,  and 
earthquake  case  were  13,000,  970,  and  7,700  times  the  energy  flows 
required  to  generate  disorder  when  measured  in  heat  equivalents.  When 
expressed  in  embodied  solar  equivalents,  the  ratio  of  ordering  energy 
to  disordering  energy  was  found  to  be  4:1  for  radiation  disordering 
and  from  25:1  to  35:1  for  pulse  disordering. 

The  theory  provided  suggestions  for  land-use  policy.  Energy  ratios 
that  provide  a  quantitative  basis  for  disaster  planning  can  be  developed 
for  a  local  environment  of  pulsing  energy. 

Possibilities  were  considered  that  cycles  of  order  and  disorder  of 
the  earth  are  synchronized  by  cycles  of  sunspots.  Energy  quality  and 
pulse  amplifier  ratios  of  solar  flares  may  be  high  enough  to  control 
many  global  cycles. 
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INTRODUCTION 

A  major  problem  in  environmental  science  is  the  nature  of  cycles 
of  order,  disorder,  and  disasters.  How  are  these  related  to  the  solar 
energies  that  drive  most  processes  on  earth?  How  do  the  patterns  of 
the  biosphere  limit  the  roles  of  humanity  and  provide  planning  prin- 
ciples? In  this  dissertation,  energy  analysis,  energy  systems  con- 
cepts, and  simulation  methods  are  used  to  test  theories  relating 
order  and  disorder.  Included  are  evaluations  of  earthquakes,  floods, 
ionizing  radiation,  solar  flares,  and  the  world  energy  web.  Specific 
case  studies  include  the  gamma  irradiation  of  a  tropical  rain  forest, 
the  Johnstown  flood  of  1889,  and  the  Guatemala  earthquake  of  1976. 

The  ability  to  converge  energy  in  space  and  time  by  cycling 
matter  may  be  a  fundamental  principle  of  ail  systems  of  man  and  nature. 
As  systems  build  structure  by  concentration  of  energy  in  the  form  of 
ordered  matter,  potential  energy  gradients  are  created.  The  storages 
of  energy  that  are  formed  become  sources  of  even  steeper  potential 
gradients.  Extensive  competition  exists  for  the  use  of  concentrated 
forms  of  energy.  The  principle  of  selection  for  maximum  power  sug- 
gests that  the  system  that  can  process  storages  most  effectively  will 
out  compete  the  other  systems.  When  energy  is  processed  at  yery   high 
rates, pulses  are  generated  that  are  transmitted  against  the  concentra- 
tion gradient  performing  a  control  on  the  surrounding  areas  and  on  the 
energy  convergence  web.  These  surges  are  often  regarded  as  disasters. 
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Theory  of  Order,  Disorder,  and  Material  Cycling 

Figure  la  is  the  basic  configuration  of  the  order-disorder  model 
diagramed  in  energy  circuit  language  (H.  T.  Odum,  1971)  that  was  used 
throughout  this  research  (see  Methods  section).  The  system  converges 
energy  to  power  the  production  of  structure  by  transforming  disordered 
matter  into  ordered  state  (lower  entropy);  however,  the  entropy  of  the 
total  process  increases  as  much  as  the  energy  required  to  power  the  pro- 
duction process  is  degraded  as  waste  heat.  The  structure  built  up  in 
the  production  process  continuously  undergoes  depreciation  in  such  a 
manner  that  the  stored  energy  is  dissipated  as  waste  heat  and  the 
material  recycles  to  form  disordered  matter. 

The  energy  circuit  model  in  Figure  lb  is  similar,  except  an  external 
energy  source  drives  the  consumption-recycle  process.  In  this  model  the 
rate  of  material  recycle  is  proportional  to  the  magnitude  of  the  stress 
energy  source.  In  both  of  the  models  in  Figure  1,  the  matter  is  cycled 
as  energy  flow  through  the  system.  In  one  the  energy  driving  recycle  is 
only  from  its  own  storages. 

A  basic  question,  analyzed  in  this  dissertation  and  illustrated  by 
the  order-disorder  model  in  Figure  1,  is  "what  relationship,  if  any,  exists 
between  the  energy  required  to  order  a  system  versus  that  required  to  dis- 
order it."  In  order  to  address  this  question,  the  model  was  evaluated 
using  data  from  several  examples  and  simulated.  Examples  included  irradia- 
tion of  a  tropical  rain  forest  at  El  Verde,  Puerto  Rico, and  the  catastrophic 
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Figure  1.  Energy  circuit  models  of  order  and  disorder. 

(a)  Energy  circuit  of  basic  process  of  order 
and  disorder  believed  to  be  characteristic 
of  all  systems  of  man  and  nature;  recycle 
uses  stored  energy. 

(b)  Basic  order-disorder  model  with  disordering 
process  driven  by  external  energy  source. 
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destructions  of  Johnstown,  Pennsylvania,  and  El  Progreso,  Guatemala. 
Catastrophic  Consumption  and  Pulsing 

As  systems  build  up  substantial  storages  of  energy  in  the  struc- 
ture form  of  ordered  matter,  massive  potential  gradients  are  formed.  As 
the  gradient  becomes  larger,  the  production  of  order  by  the  system  often 
exceeds  the  loss  due  to  depreciation  through  linear  processes  such  as 
diffusion.  This,  therefore,  increases  the  possibility  and  opportunity 
for  a  consumer  system  to  form  that  will  further  concentrate  the  stored 
energy  and  release  a  pulse  of  energy  in  a  control  action.  To  explain 
this  concept,  the  order-disorder  models  in  Figure  2,  similar  to  the  pre- 
viously introduced  model,  were  developed.  An  autocatalytic  or  self- 
generating  consumer  system  was  added  in  the  disordering  energy  flow 
path  to  further  concentrate  stored  energy.  In  the  system  diagramed  in 
Figure  2a,  a  trigger  pulse,  or  seed,  is  required  to  start  the  consumer 
system;  however,  once  triggered  the  system  will  grow  at  an  exponential 
rate  until  the  energy  storage  in  the  form  of  order  in  the  producer 
system  is  consumed.  This  growth  of  the  consumer  system  produces  an 
energy  control  pulse  that  is  exponential  at  first  and  fades  as  the 
storages  are  exhausted. 

The  model  in  Figure  2b  is  similar  except  the  system  is  capable  of 
self- triggering;  that  is,  when  the  energy  storage  in  the  producer  sys- 
tem reaches  threshold,  the  consumer  system  self-activates.  Consumer 
systems  r.hat   can  process  large  quantities  of  energy  at  high  rates  pro- 
duce catastrophic  pulses  of  energy  when  activated.  The  order-disorder 
models  in  Figure  2  are  suggested  as  a  dynamic  system  that  can  simulate 
the  concentration,  storage,  and  catastrophic  release  process  found  in 
many  systems  of  man  and  nature,  of  which  natural  disasters  are  a  special 


Figure  2.  Energy  circuit  models  of  order-disorder  cycle  with 
catastrophic  energy  pulsing  in  the  disordering 

loop. 

(a)  Externally  triggered  disaster. 

(b)  Self- triggered  disaster. 
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case.  Data  from  the  catastrophic  collapse  of  a  dam  resulting  from  the 
exponential  growth  of  turbulence  were  used  to  evaluate  the  model.  Sim- 
ulations were  performed  and  the  results  compared  to  field  observation 
for  the  Guatemala  earthquake  and  the  Johnstown  flood. 
Secondary  Disorder  and  Cascading 

Two  order-disorder  models  with  catastrophic  consumer  systems  were 
cascaded  to  extend  the  theory  to  explain  and  simulate  the  secondary 
disruptions  created  by  a  major  disaster  pulse  such  as  an  earthquake. 
Once  the  surge  of  consumption  starts  it  may  generate  waves  of  disor- 
dering energy  spreading  into  surrounding  energy  convergence  systems. 
Figure  3  is  an  energy  circuit  diagram  of  this  concept.  In  the  model, 
the  energy  storage  increases  until  the  energy  flow  in  the  disordered 
loop  is  great  enough  to  generate  a  dynamic  structure,  which  feeds 
back  some  of  the  energy  to  pump  more.  This  results  in  exponential 
growth  of  the  storage  of  kinetic  energy  until  the  potential  energy  is 
effectively  exhausted.  When  sufficient  energy  is  coupled  to  the 
structure  of  the  impacted  system,  the  potential  energy  stored  in  the 
impacted  system  is  released  causing  a  secondary  disaster.  Some  notable 
examples  of  secondary  events  triggered  by  earthquakes  are  given  by 
Ayre  (1975):  (1)  failure  of  man-made  structures  such  as  buildings  and 
bridges;  (2)  avalanches,  landslides,  rocks! ides,  and  other  natural  dis- 
ruptions; (3)  fire;  (4)  floods  and  flash  floods  from  tidal  waves  and 
ruptured  dams;  and  (5)  volcanic  eruption.  Note  that  each  of  the  above 
events  involves  the  cascaded  release  of  potential  energy  stored  in  the 
structure  of  the  secondary  source.  The  cascading  of  disasters  is 
believed  to  be  of  considerable  importance  in  understanding  the  dynamic 
behavior  of  systems  under  the  high  stress  of  natural  disaster  pulses. 


Figure  3.  Energy  circuit  model  of  concept  of  cascading  of 
energy  pulses  which  may  be  regarded  as  disasters, 
to  smaller  elements  of  the  systems. 
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The  multiple  hazard  phenomena  must  be  included  in  a  disaster  theory  to 
make  it  plausible  for  application  to  real  world  cases. 

This  theory  was  tested  by  evaluation  of  the  model  with  data  from 
an  earthquake  and  the  secondary  destruction  of  several  towns.  Simula- 
tions were  performed  and  compared  to  observed  destruction. 

Global  Model 


The  global  model  in  Figure  4  connects  energy  concentration,  stor- 
age, and  consumer  subsystems  into  a  web.  This  energy  web  provides  a 
framework  that  helps  to  explain  the  nature  and  function  of  disasters 
in  a  global  context.  The  model  also  outlines  a  conceptual  network  for 
evaluating  the  flows  and  storages  of  energy  that  power  various  global 
processes.  The  web  is  composed  of  self-regulating  units  that  success- 
fully converge  energy  of  the  sun.  Each  module  produces  order  and 
releases  it  with  a  pulsed  consumption  that  forms  a  temporary  high 
energy  system.  To  survive  in  competition  this  pulse  of  energy  must 
couple  with  or  feed  back  to  other  systems.  Pulsing  the  energy  that  is 
fed  back  in  a  control  action  produces  an  amplifier  effect  in  disor- 
dering and  recycling  other  systems.  Systems  that  survive  are  those 
that  become  organized  and  adapted  to  operate  in  a  regime  of  pulses 
and  use  these  pulses  to  power  recycling  processes  and  other  means  to 
achieve  maximum  power. 

The  Global  Energy  Web  and  Natural  Disasters 

Natural  disasters  have  been  a  subject  of  great  interest  to  human- 
ity throughout  recorded  history.  Earthquakes,  volcanoes,  floods,  storms, 
fires,  and  epidemics  take  a  great  toll  in  life  and  property.  Hazards  in 
the  United  States  currently  generate  an  economic  less  in  excess  of  one 
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billion  dollars  per  year,  with  all  data  indicating  that  this  will  become 
an  ever-growing  problem  (White  and  Haas,  1975).  Pulse  disasters  are  a 
major  part  of  all  systems  of  nature  with  catastrophic  energy  releases, 
often  with  cascading  effects,  causing  many  smaller  systems  to  be  disor- 
dered. Energy  concepts  of  order  and  disorder  provide  a  systematic  means 
of  disaster  analysis  and  help  explain  and  generalize  about  natural  di- 
sasters. The  theory  is  that  global  processes  concentrate,  store,  and 
catastrophically  release  energy  according  to  a  common  pattern, thus 
creating  natural  disasters. 

Energy  Quality  and  Embodied  Energy 

Energy  quality  theory  as  developed  by  H.  T.  Odum  (1971)  and  re- 
cently amplified  (Odum,  1978)  was  built  on  the  laws  of  energy.  The  law 
of  energy  conservation  provides  the  basic  underpinnings  for  all  energy 
circuit  diagrams.  In  a  first  law  energy  diagram,  the  average  heat 
equivalent  stored  or  flowing  per  unit  of  time  is  used  to  evaluate  the 
model.  When  evaluated  in  heat  equivalents,  all  systems  of  man  and 
nature  conserve  energy,  that  is,  the  sum  of  all  the  flows  into  a  sys- 
tem less  the  flows  of  energy  out  must  equal  the  net  changes  in  the 
energy  storages.  The  system  models  in  this  dissertation  meet  the  first 
law  requirement.  Unfortunately,  heat  equivalent  measures  provide  no 
information  relating  to  the  potential  value  of  energy  flow  or  storage 
for  performing  work  or  controlling  a  system. 

Working  with  the  law  of  energy  degradation  provides  insight  into 
assigning  a  value  to  a  particular  energy  flow  or  storage.  The  law  of 
energy  degradation  states  that  in  all  useful  processes  some  energy 
must  be  degraded,  thus  losing  its  ability  to  do  work.  The  concept  is 
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incorporated  into  energy  circuit  language  by  requiring  heat  sinks  or 
energy  degradation  pathways  on  all  energy  flows,  processes,  and  stor- 
ages. The  quantity  of  energy  dissipated  by  each  energy  concentration 
process  is  also  evaluated  in  heat  equivalents  in  energy  circuit  dia- 
grams, even  though  these  energy  pathways  can  do  no  useful  work  on  this 
earth. 

The  concept  of  energy  quality  places  an  embodied  energy  measuring 
cost  in  energy  units  on  the  concentrated  flows  and  storages  which  are 
a  resultant  of  the  energy  convergence  property  of  systems.  Energy 
quality  factors  for  systems  are  defined  to  be  the  ratio  of  the  heat 
equivalent  energy  produced  or  upgraded  by  a  system  to  the  quantity  of 
energy  required  to  power  the  system. 

The  maximum  power  principle  makes  it  possible  to  relate  value  cal- 
culated on  a  basis  of  energy  used  to  value  based  on  effect  of  energy 
flow.  The  maximum  power  principle  proposes  that  systems  that  maximize 
their  flows  of  energy  survive  in  competition  (Lotka,  1922).  The  theo- 
retical relationship  between  cost  quality  and  effect  quality  is  grounded 
in  the  idea  that  the  surviving  system  must  effectively  use  its  energy 
flows  and  storages  to  survive  in  competition.  Hence,  in  the  long  run, 
energy  cost  must  equal  energy  effect. 

Solar  energy  is  used  as  the  basis  for  all  energy  quality  calcula- 
tions in  this  dissertation.  By  constructing  a  global  energy  web  model, 
which  was  evaluated  in  heat  equivalent,  the  energy  quality  of  the  major 
energy  convergence  and  control  processes  of  the  world  were  calculated 
based  on  the  quantity  of  solar  energy  required  to  power  the  system  as 
suggested  by  Odum  (1978). 
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Energy  Embodied  in  Radiation 

A  method  is  proposed  for  calculating  the  energy  embodied  in  electro- 
magnetic radiation.  For  calculations  in  the  biosphere  web  energy  quality 
factors  are  calculated  based  on  the  quantity  of  energy  required  to  pro- 
duce a  certain  energy  in  open  systems  flow.  Embodied  energy  is  then  cal- 
culated by  multiplying  measured  energy  flows  in  heat  calories  by  the 
energy  quality  factor  (Odum,  1978). 

In  the  proposed  method  for  the  energy  embodied  in  electromagnetic 
radiation,  energy  quality  factors  are  based  on  the  absolute  temperature 
of  the  radiation  energy  source.  The  concept  is  based  on  the  fact  that 
more  energy  is  required  per  degree  Kelvin  to  raise  the  temperature  of 
the  population  of  molecules.  The  mean  energy  of  the  photon  emitted  from 
radiation  peak  of  the  energy  source  is  proportional  to  the  temperature 
of  the  source.  Wien's  displacement  law  states  that  the  photon  wave- 
length maximum  emission  A^x  is  related  to  absolute  temperature  of  the 
radiating  body  by 

Amax  =  T  /Ko\  (1) 

where  R  =  2.89  x  10"3  m  °K  (Rosenberg,  1974).  Given  the  wavelength  of 
electromagnetic  radiation,  the  energy  per  photon  e  may  be  calculated 
from  Planck's  constant. 

emax  =  h  |  (2) 

max 

where  Planck's  constant  h  =  1.58  x  10   Cal  sec  and  c  is  the  velocity 
of  light  (c  =  3  x  108  m  sec"1)- 

Substituting  equation  2  into  equation  1  produces  a  direct  rela- 
tionship between  the  energy  of  the  photons  emitted  from  the  peak  of 
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the  electromagnetic  spectrum  and  the  absolute  temperature  of  the  radi- 
ating body. 

emac  =  k2  T  (°K)  Cal  (2a) 

-20        - 1 

where  k2  =  1.6  x  10    Cal  °K  .  Thus  the  average  energy  of  the  photons 
emitted  from  the  peak  of  the  radiation  spectrum  of  a  blackbody  is  in 
proportion  to  the  temperature  and  the  constant  of  proportionality  is 
known . 

It  is  postulated  that  in  order  to  build  a  system  that  can  support 
and  sustain  high  temperature,  some  form  of  molecular  energy  convergence 
chain  must  be  formed.  The  total  energy  required  to  support  the  system 
increases  in  proportion  to  the  temperature  obtained. 

Since  the  energy  per  photon  is  proportional  to  the  temperature  of 
the  source,  it  may  be  possible  to  relate  the  energy  quality  of  certain 
energy  photons  to  the  energy  in  the  population  of  molecules  required 
to  produce  the  photons. 

The  blackbody  temperature  of  the  sun,  6,000°K  (9.8  x  10"17  Cal 
photons)  formed  a  reference  for  energy  quality  calculations  per  photon . 
For  example  the  earth's  mean  surface  temperature  is  about  300°K;  thus, 
at  its  spectral  peak  it  gives  off  4.92  x  10"   Cal  photons,  which  have 
an  energy  quality  factor  0.05.  A  60,000°K  sunspot  would  emit  photons 
with  a  relative  energy  quality  factor  of  10. 

Divergence  of  radiation  is  not  directly  addressed  by  the  proposed 
method  of  calculating  the  embodied  energy  of  radiation.  The  question 
of  the  relationship  between  photon  flux  density  and  the  quantity  of 
energy  required  to  produce  the  radiation  is  not  answered. 
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Previous  Work  on  the  Theory  of  Order,  Disorder,  and  Pulsing 

Previous  study  of  order-disorder  theory  was  developed  with  energy 
circuit  language,  a  systems  diagram  and  a  mathematical  framework, 
by     H.  T.  Odum  (1971).  A  general  model  of  order  and  disorder  was 
recognized  as  Michael is-Menton  Kinetics.  Odum  and  fellow  researchers 
have  used  energy  circuit  language  to  study  and  model  specific  catastro- 
phic events  previously.  Jeff  Richey  (1970)  disordered  microcosms;  Odum 
(1970c)  measured  order-disorder  relations  with  the  irradiation  of  a  trop- 
ical rain  forest  and  considered  order-disorder  model  of  war  and  the  ef- 
fects of  herbicides  in  South  Vietnam  (1974).  Energy  in  food  chains  and 
webs  were  found  to  maximize  power  for  survival  by  conserving  and  storing 
energy  in  high  quality  consumers  and  releasing  energy  in  pulses. 

A  massive  report  on  the  Rain  Forest  Research  Project  was  prepared 
under  H.  T.  Odum's  direction  by  the  multidisciplinary  scientific  field 
team.  Two  parts  of  the  report  relate  directly  to  the  concept  of  order 
and  disorder  as  applied  to  a  natural  ecosystem.  Dienes  (1970)  developed 
a  kinetic  model  of  the  radiation  response  of  the  forest  ecosystem  to 
field  irradiation  based  on  a  chemical  reaction  type  system  of  equilib- 
rium equations  as  in  equation  3: 

Ki        K2 
Xl   -,   X2  ►  X3  (3) 

K-I 

where  Xi  =  normal  state;  X?  =  partially  damaged  state;  X3  =  lethal  state; 
Ki  and  !<2  =  rate  constant  corresponding  to  radiation  induced  chance;  and 
K_1  =  radiation  independent  rate  constant  describing  recovery  of  the  sys- 
tem. Dienes  used  the  model  to  calculate  survival  versus  distance  from 
radiation  source  for  various  sets  of  hypothetical  coefficients.  In  H.  T. 
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Odum's  (1970c)  ecological  analysis  of  the  system  at  El  Verde,  the  clas- 
sic washout  curve  of  the  chemostat  was  used  to  explain  the  effect  of 
the  radiation  stress,  yielding  similar  results  to  Dienes'  model  and  the 
observed  effect  on  natural  system.  Probably  more  significantly,  however, 
from  the  standpoint  of  this  research,  the  ratio  of  disordering  energy  P: 
to  potential  energy  flux  for  reordering  the  system  was  calculated  by 
Odum  who  found  that  in  a  zone  where  destructive  action  balanced  repair 
rates  the  gamma  ray  energy  input  was  only  about  1/10,000  of  that  due  to 
photosynthesis  (Odum,  1970c). 

A  recent  innovation  in  mathematics  is  a  catastrophe  theory  by  Rene 
Thorn  (1975).  In  Thorn's  theory,  models  for  elementary  catastrophes  are 
constructed  from  higher  order  differential  equations.  Surfaces  are 
generated  from  all  the  points  of  the  functions  where  the  first  deriva- 
tive is  equal  to  zero.  Thorn  produced   fold,  cusp,  sallowtail,  butter- 
fly, hyperbolic,  elliptic,  and  parabolic  surfacesin  this  manner.  Pos- 
sibly the  interesting  point  of  Thorn's  work  is  the  use  of  higher  order 
terms  in  the  differential  equations  describing  the  surface  that  are 
similar  to  the  result  of  interacting  several  energy  flows  in  a  web 
model . 

A  theory  related  to  the  theory  of  order  and  disorder  developed  here- 
in is  the  theory  of  risk  as  it  is  used  in  economics.  Hirshlerfer  and 
Shapiro  (1970)  analyzed  the  appropriate  discount  rate  for  capital  invest- 
ments. In  this  paper  risk  of  catastrophic  loss  is  treated  as  a  deprecia- 
tion flow  similar  to  technological  obsolescence  and  physical  depreciation. 
Risk  probabilities  are  based  on  historical  or  perceived  probabilities  of 
disaster,  such  as  depression,  war,  or  drought.  In  mathematical  form 
Hirshleifer  and  Shapiro  incorporate  risk  into  the  standard  calculation 


of  the  present  value  of  one  dollar  of  future  money,  as  given  in  equation 
4,  by  increasing  the  standard  discount  rate  by  risk  factor. 


P=     1 


.1  +  i  +  r)  (4) 

Where  i  =  rate  of  interest  and  r  is  the  rate  of  risk. 

Thurow  and  Taylor  (1966)  proposed  a  model  for  generation  of  the 
capital  stock  K  as  given  in  equation  5: 

K  =  I  -  dK  -  aK  (5) 

where  I  is  investment,  d  is  physical  depreciation,  and  a   is  technologi- 
cal obsolescence.  Incorporating  risk  into  Thurow  and  Taylor's  model 
yields  equation  6: 

K  =  I  -  dK  -  aK  -  rK  (6) 

where  r  is  the  rate  of  risk. 

Figure  5  is  a  translation  of  the  Thurow  and  Taylor  model  into 
energy  circuit  language  drawn  in  the  format  of  an  order-disorder  model 
of  Figure  1.  Note  that  depreciation  and  obsolescence  are  driven  by 
internally  built  accumulation  of  capital  assets  while  risk  is  an  ex- 
ternally powered  process. 

Previous  Work  on  Natural  Disasters 

Relatively  few  studies  provide  a  basic  theory  of  natural  disasters. 
Most  research  has  dealt  with  descriptions  and  analyses  of  specific 
events,  such  as  reports  on  the  Alaska  earthquake  and  tsunami  on  Good 
Friday,  March  27,  1964  (Rogers,  1970),  Hurricane  Agnes  of  June  1972 
(National  Advisory  Committee  on  Oceans  and  Atmospheres,  1972),  and  the 
Lubbock,  Texas,  storm  of  May  11,  1970  (Thompson,  1970).  Another  common 


Figure  5.  Energy  circuit  translation  of  economic  model  of  risk. 
The  diagram  is  based  on  Thurow  and  Taylor  (1966) 
capital  stock  generation  model. 
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K=   I  -  (d+ajK  -rK 
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disaster  research  procedure  is  to  look  at  one   component   in  several 
types  of  disasters.  Notable  examples  would  be  Dacy  and  Kunreuther's 
(1969)  economic  analysis  of  natural  disaster  economic  disruption  and 
Cochrane 's  (1975)  analysis  of  the  disrupting  effects  of  natural  hazards. 
Dacy  and  Kunreuther  developed  an  economic  framework  for  analyzing  the 
economic  problems  resulting  from  disasters,  analyzed  short-term  recovery 
across  several  case  catastrophes,  looked  at  empirical  evidence  of  long- 
term  disaster  recovery,  and  tried  to  describe  the  report  of  federal 
policy  on  disaster  mitigation.  Cochrane  focused  on  what  populations 
were  destroyed  and  how  much  of  the  burden  of  loss  and  reconstruction  was 
borne  by  whom. 

Approaches  that  are  not  disaster  theory  per  se  do  merit  recognition 
including  Ericksen's  (1975)  scenario  methodology  for  natural  hazard  re- 
search and  Friedman's  (1975)  computer  simulation  of  natural  hazard 
assessments.  Ericksen  synthesizes  observed  results  from  cross-sectional 
analysis  of  disasters  at  progressive  longitudinal  steps  to  qualitatively 
describe  what  might  happen  to  specific  communities  in  future  disasters. 
Ericksen's  methodology  was  based  on  the  work  of  Kahn  and  Wiener  (1967). 
Friedman  (1975)  had  a  similar  but  more  quantitative  approach  in  using 
computer  mapping  to  simulate  the  possible  risk  of  destruction  of  future 
earthquakes,  hurricanes,  and  inland  flooding  for  various  geographic 
areas  of  the  United  States.  Friedman's  interest  was  stimulated  by  his 
employment  by  the  Traveler's  Insurance  Company  and  their  concern  over 
the  spatial  distribution  of  loss  of  life  and  property  in  the  future. 
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Description  of  Study  Situations  and  Previous  Work 

In  this  section  previous  work  on  energy  cycles  of  the  earth  are 
reviewed  and  related  to  natural  disasters.  The  sites  and  systems  of 
special  study  and  the  disordering  events  of  interest  are  described 
including  the  tropical  rain  forest  at  El  Verde,  Puerto  Rico,  and  the 
disordering  effects  of  radiation;  the  setting  for  the  City  of  Johnstown, 
Pennsylvania,  1889,  and  the  disastrous  effects  of  the  failure  of  the 
South  Fork  Dam;  and  the  town  of  El  Progreso,  Guatemala,  and  the  catas- 
trophic effect  of  the  earthquake  of  February  4,  1976. 

Global  Models  of  Atmospheric,  Oceanic,  Biological,  Geological,  and  Urban 
Systems 

General  circulation  models  of  the  atmosphere  are  normally  developed 
using  a  grid  technique  integrating  the  equations  over  large  regions, 
often  a  hemisphere  or  the  whole  globe,  using  a  finite  differencing 
scheme  with  as  small  as  practical  integration  step.  The  grid  might  be 
200  km  square  horizontally  and  divided  into  as  many  as  18  layers  verti- 
cally (Frisken,  1973).  Significant  examples  of  uses  of  this  technique 
are  found  in  Smagorinsky  (1969)  and  Oliger  et  al_. ,  (1970).  General  cir- 
culation models  attempt  to  simulate  the  daily  dynamic  process  of  the 
earth's  atmosphere,  while  cl imatological  models  generally  tend  to  con- 
centrate more  heavily  on  modeling  annual  means. 

Houghton  (1954)  developed  a  schematic  representation  of  the  flows 
of  energy  in  the  atmosphere.  Houghton  estimated  that  the  transport  of 
energy  toward  the  poles  by  the  atmosphere  and  by  the  sea  would  reach  a 
maximum  of  1  x  1017  Cal  d"1  across  the  40°N  latitude  circle.  Budyko 
(1969)  and  Sellers  (1969)  expanded  Houghton's  schematic  model  to  a  math- 
ematical mean  annual  zonally  averaged  energy  balance  model.  The  effect 
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of  solar  radiation  on  variations  on  the  climate  of  the  earth  is  a 
central  feature  of  the  model.  Further  development  of  the  theory  of 
energy  balance  climate  models  by  North  (1975)  has  shown  that  a  small  varia- 
tion of  the  solar  mean  can  trigger  an  ice  age  in  the  model.  Somewhat 
parallel  to  the  developments  of  climate  models,  considerable  attention 
has  been  placed  on  the  development  of  global  biological  productivity 
estimates. 

Estimates  of  global  production  by  Liebig  (1862),  which  were  extra- 
polated from  a  single  community,  were  refined  by  Schroeder  (1919)  by 
distinguishing  four  land  community  types.  Helmut  Lieth  (1975)  developed 
a  spatial  model  of  the  primary  productivity  of  the  biological  systems  of 
the  earth-based  world  climate  data  that  included  20  vegetation  types. 
These  were  not  dynamic;  that  is,  there  is  no  feedback  from  the  vegetation 
cover  of  the  earth  to  the  global  climate.  The  patterns  of  vegetation  of 
the  earth  that  were  analyzed  and  tabulated  by  Whittaker  and  Likens  (1975) 
provided  a  useful  summary  of  the  production  and  stocks  in  the  global 
system. 

Another  approach  to  global  biological  system  modeling  has  been  in 
the  carbon  balance  of  the  world.  Bolin  (1970)  diagramed  a  carbon  model 
of  the  biosphere,  finding  the  vast  majority  of  carbon  trapped  in  sedi- 
ments amounting  to  39  metric  tons  per  square  meter  of  the  earth's  sur- 
face. However,  only  a  few  tenths  of  a  percent  of  the  immense  mass  of 
carbon  is  on  cr  near  the  earth's  surface.  By  estimating  detritus  by 
ecosystem  type  for  the  earth,  Schlesinger  (1977)  found  2.8  kg  m  ~  of 
carbon  for  the  world's  surface  in  detritus. 

Mathematical  models  of  the  sea  are  usually  restricted  to  individual 
components  of  marine  systems  running  on  wery   large  data  bases  with 
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little  effort  toward  development  of  a  comprehensive  model  (Walsh,  1972). 
Ekman  (1905)  showed  wind  stress,  acting  at  the  sea's  surface,  put  momen- 
tum into  the  ocean  in  a  thin  boundary  layer.  Ekman's  theory  is  normally 
included  in  general  circulation  models  of  the  sea  (Holland,  1977). 
Sverdrup  (1947)  developed  a  vorticity  balance  theory  to  which  Stommel 
(1948)  added  bottom  friction  to  explain  closed  basis  circulation.  Holland 
and  Sverdrup 's  work  was  the  basis  for  general  circulation  models  of  the 
sea.  Models  of  the  ocean's  tide  were  based  on  the  incorporation  of 
Newton's  second  law  by  Laplace  in  1776  to  explain  the  action  of  gravity 
on  the  sea  (Hendershott,  1977).  Seasonal  models  of  the  thermoclines 
attempt  to  explain  the  heat  exchange  between  the  atmosphere  and  the 
ocean  that  produces  large  horizontal -scale  changes  at  the  temperature 
structure  in  surface  layers  of  the  ocean  (Niler,  1977).  Charney  (1955), 
who  developed  a  theory  to  explain  geostrophic  adjustment,  used  the 
theory  coastlines  known  as  the  costal  jets  such  as  the  Gulf  Stream. 

Modeling  of  the  sea  is  often  fragmented  among  circulation,  waves, 
nutrient,  productivity,  zooplankton,  and  phytoplankton  models  as  if 
they  were  independent  processes.  Odum  (1975)  considered  ways  of  com- 
bining physical,  chemical,  and  biological  units  in  a  single  model  of 
the  sea. 

Sand  transport  on  beaches  and  shoreline  evolution  geological  sim- 
ulation techniques  have  been  documented  by  Harbaugh  and  Bonham-Carter 
(1970).  The  processes  of  erosion,  transport,  deposition,  and  consoli- 
dation form  the  basis  for  models  of  marine  sediment  models. 

Dynamic  modeling  of  the  earth's  crust,  supported  by  the  Vine- 
Matthews  Hypothesis,  proposes  that  anomalies  in  the  earth's  magnetic 
field  are  due  to  cycles  in  the  earth's  magnetic  field.  The  cycles  sec 
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the  magnetic  properties  of  the  basaltic  rocks  of  the  earth's  ocean  crust 
as  they  formed  on  the  ocean  ridges.  Magnetic  patterns  in  basaltic  rocks 
in  the  crust  underlying  the  sea  floor  and  the  age  of  sediments  on  the 
sea  floor  provide  conclusive  evidence  of  the  sea  floor's  spreading  from 
the  ocean  ledges  at  approximately  2  cm  yr~l   (Clark,  1971).  Hayakawa 
(1976)  included  continental  drift,  midocean  ridge,  volcanoes,  orogenesis, 
and  other  geological  phenomena  in  a  pictoral  energy  balance  model  of  the 
earth's  crust.  Hayakawa' s  model  was  evaluated  by  estimating  energy  flows 
from  geophysical  data. 

A  dynamic  global  web  model  of  human  systems  and  resources  was  made 
by  J.  Forrester  (1970)  with  the  support  and  encouragement  of  the  Club 
of  Rome.  Forrester's  work  with  world  modeling  was  followed  by  an  ex- 
tensively documented  work  by  Meadows  et  al_.  (1974).  Eoth  models  were 
diagramed  in  systems  dynamics  language  (Forrester,  1970)  and  simulated 
with  DYNAMO  (Pugh,  1970). 

Hubbert  (1971)  evaluated  a  schematic  global  energy  web  highlighting 
flows  in  the  atmospheric  system  with  some  attention  to  aggregate  biolog- 
ical, geological,  and  oceanic  energy  flows.  In  Hubbert' s  work,  concen- 
tration and  storages  of  energy  are  central,  especially  those  in  fossil 
fuels.  However,  the  only  catastrophic  releases  of  energy  considered  are 
geologic  in  nature.  H.  T.  Odum  has  developed  several  global  energy 
models  (Odum,  1972;  1973).  A  global  minimodel  was  developed,  evaluated, 
and  simulated  by  Alexander  et_  aj_.  (1976)  to  estimate  long-term  United 
States  fossil  fuel  availability.  A  recent  energy  web  model  of  the  bio- 
sphere (Odum,  1978)  treats  storms  and  volcanoes  as  high  quality  control 
mechanisms  which  consume  energy  concentrated  by  the  atmospheric  and 
continental  sedimentary  cycles,  respectively. 
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Tne  Tropical  Rain  Forest  at  El  Verde,  Puerto  Rico  and  the  Effect  of 
Radiation  ~~         ~" 

As  part  of  a  major  postwar  research  effort  by  the  United  States  to 
quantify  the  effects  of  irradiation  on  the  environment,  a  field  experi- 
ment was  set  up  at  El  Verde,  Puerto  Rico,  by  the  Atomic  Energy  Commis- 
sion in  1963.  A  goal  of  the  research  was  to  perform  a  detailed  study 
of  the  disordering  effect  of  gamma  radiation  on  a  tropical  rain  forest 
ecosystem. 

Radiation  studies  with  similar  design  were  done  at  Brookhaven 
National  Laboratory  and  at  Rhine! ander,  Wisconsin.  At  Brookhaven,  the 
chronic  effects  of  irradiation  on  terrestrial  ecosystems  were  studied 
(Woodwell,  1963).  Woodwell  and  Rebuck  (1967)  analyzed  the  effects  of 
chronic  gamma  irradiation  on  oak  and  pine  forests.  Canoy  (1972), 
writing  a  dissertation  under  the  direction  of  H.T.  Odum,  measured  the 
deoxyribonucleic  acid  (DNA)  of  radiated  vegetation  at  Brookhaven  and 
at  El  Verde  by  relating  molecular  order  and  information  to  the  pattern 
of  stress.  Highest  DNA  levels  were  found  at  boundary  where  ordering 
and  disordering  rates  were  similar. 

The  radiation  source  used  in  Puerto  Rico  under  Dr.  Odum's  direc- 
tion was  later  moved  to  Enterprise  Radiation  Forest  near  Rhinelander, 
Wisconsin,  where  a  similar  exposure  of  a  major  northern  forest  type 
(aspen)  was  conducted  (Rudolph,  1974;  Zavitkovski,  1977).  In  the  pre- 
irradiation  study,  detailed  analyses  of  the  ecological  community-- 
including  solar  radiation,  lichens,  flora,  plant  communities,  temporal 
and  spatial  patterns,  leaf-litter  production,  and  vertebrates  and 
small -mammal  populations—were  performed.  The  effects  of  gamma  irradia- 
tion on  northern  forest  communities  was  predicted  by  Zavitkovski  et  al_. 
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(1974).  In  the  recently  completed  radioecological  studies,  extensive 
investigations  were  made  into  the  ecological  effects  of  gamma  radia- 
tion. Of  special  interest  to  this  research  was  the  radiation  effect 
on  biomass  production  of  ground  vegetation  studied  by  Zavitkovski  and 
Salmonson  (1977)  and  the  response  of  the  forest  ecotone  to  ionizing 
radiation  studied  by  Murphy  et  al_.  (1977).  A  result  of  their  research, 
which  relates  directly  to  the  theory  of  order  and  disorder,  was  the 
large  increase  in  the  number  of  seedlings  they  found  in  the  irradia- 
tion recycled  area. 

Puerto  Rico  is  a  Caribbean  island  located  between  18°N  and  19°N 
latitude  and  65°W  longitude.  El  Verde  station  was  constructed  in 
the  tropical  rain  forest  region  of  the  Luquillo  Mountains  of  eastern 
Puerto  Rico  (see  Figure  6a).  The  steady  flow  of  warm,  moist  air  from 
the  1.2  m  sec-1  prevailing  easterly  trade  wind  provides  a  uniformly 
wet  climate  with  a  mean  temperature  of  22.6eC,  a  relative  humidity  of 
91  percent,  an  absolute  humidity  of  18.7  g  m"3,  an  insulation  of  3,830 
Cal  m"2  d"1,  and  an  annual  rainfall  of  3,750  mm.  Climate  fluctuations 
were  found  to  be  small,  for  example  temperature  had  only  a  21°  to  25°C 
diurnal  range  and  a  2°  seasonal  mean  fluctuation  (Odum  et  al_. ,  1970). 

The  vegetation  of  the  area  is  the  lush  tropical  forest  traditionally 
characteristic  of  rain  forest,  with  high  diversity  of  plants  and  insects, 
buttress  roots,  broad  thin  leaves,  bromeliads,  lianas,  trunk  bark  heavy 
with  growth,  and  an  open  ground  story  in  deep  shade.  The  area  is  known 
as  the  Tabanuco  Forest  after  the  principal  tree,  the  tabanuco.  With  a 
rainfall  of  greater  than  102  mm  every  month,  droughts  were  net  experienced 


Figure  6.  Location  maps  and  diagram  of  Puerto  Rico  and  El 
Verde  experiment. 

(a)  Map  of  Puerto  Rico  showing  east  location  of 
El  Verde  (Odum,  1970a). 

(b)  Main  base  map  for  studies  at  the  El  Verde 
site  (Odum,  1970a). 

(c)  El  Verde  station  and  vicinity  showing  some 
principal  locations  of  activity  (Odum, 
1970a). 
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by  the  forest  ecosystems.  For  example,  the  tropical  rain  forest  in  the 
El  Verde  study  area  had  never  experienced  fire  (Odum,  1970a). 

The  rain  forest  project  study  was  concentrated  around  the  radia- 
tion source  and  two  control  locations.  Figure  6b  is  a  plan  view  of  the 
irradiation  and  the  south  control  area  while  Figure  6c  diagrams  the 
ecological  and  irradiation  experiment.  The  irradiation  was  provided  by 
a  10,000  curie  cesium  source  for  93  days  in  early  1965.  Another  6  years 
were  required  for  analyzing  the  results  of  the  irradiation  stress  on 
the  ecosystem,  including  comparison  studies  with  two  control  sites. 
The  south  control  site  was  undisturbed  while  the  north  control  was 
denuded  of  vegetation  to  resemble  the  irradiation  site.  Substantial 
stress  was  imposed  on  the  ecosystem  by  the  high  field  gamma  radiation. 
Johnstown,  Pennsylvania,  and  the  Failure  of  the  South  Fork  Dam 

Johnstown,  Pennsylvania,  is  a  heavily  industrialized  city  located 
121  km  east  of  Pittsburg,  Pennsylvania,  in  a  valley  of  the  Allegheny 
Mountains  (see  Figure  7a).  The  town  was  founded  at  the  confluence  of 
two  medium-sized  mountain  streams,  Stony  Creek  and  the  Little  Conemaugh 
River,  which  join  at  central  Johnstown  to  become  the  Conemaugh  River. 
Above  Johnstown,  the  two  streams  and  their  tributaries  drain  650  sq  mi 
of  valley  and  mountainside  (Gilbert,  1977).  The  point  of  confluence 
provided  water  transport  and  the  level  floodplain  encouraged  development. 
These  coupled  with  the  area's  rich  coal  deposits  made  Johnstown  a  natural 
place  for  eighteenth-century  European  settlers  to  start  a  town  (Figure 
7b). 

Johnstown  received  national  attention  when  on  May  31,  1889,  6  to  8 
in  of  rain  sent  Stony  Creek  and  the  Little  Conemaugh  into  the  streets 
of  the  town  and  breached  what  was  then  the  world's  largest  earthen  dam, 


Figure  7.  Location  map  of  Johnstown  and  the  South  Fork  Dam. 

(a)  Map  of  Pennsylvania  showing  location  of 
Johnstown  in  relationship  to  other  major 
Pennsylvania  cities. 

(b)  General  map  of  the  region  of  the  Johnstown 
disaster  and  the  watershed  of  the  South 
Forks.  This  map  is  the  first  correct  one 
published,  the  decided  differences  between 
it  and  those  which  have  yet  appeared  else- 
where being  due  to  inaccuracies  in  the 
latter.  It  reproduces  by  kindness  of  G. 
W.  Colton  &  Co.  from  an  old  and  scarce 
county  map  of  1867,  on  a  scale  of  1  1/4 

in  per  mi.  The  topography  on  this  map  is 

made  accurately  to  scale  from  sketches 

made  on  the  ground  (Engineering  News,  1889a). 
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14  mi  up  the  Little  Conemaugh  at  South  Fork,  releasing  a  420-acre  lake, 
which  was  65  ft  deep  at  places  (Engineering  Society  of  Western  Pennsyl- 
vania, 1889).  The  water  in  the  lake  overflowed  the  dam  in  a  low  spot 
in  the  center  and  eroded  the  entire  900-ft-long,  70-ft-high  earthen 
structure  on  the  down  river  side,  allowing  the  earthen  wall  to  fail 
under  the  extreme  pressure,  and  sending  a  wave  descending  into  Johnstown. 
The  entire  540  million  cu  ft  of  water  was  released  behind  a  40-ft  wave 
that  took  approximately  one  hour  to  reach  Johnstown,  \lery   large  quan- 
tities of  energy  were  released  as  the  wave  roared  down  the  valley  carrying 
buildings,  trees,  houses,  factories,  and  even  locomotives.  It  hit  the 
town  with  virtually  no  warning,  scouring  most  of  the  town  to  bedrock, 
drowning  2,200  people,  and  creating  $17  million  property  damage.  Figure 
6b  shows  the  path  of  the  floodwave. 

Floods  are  nothing  new  for  Pennsylvania,  for  of  2,800  communities 
in  the  state,  2,468  are  located  in  proven  flood  zones  (Gilbert,  1977). 
However,  Johnstownians,  believing  this  to  be  a  freak  accident,  rebuilt 
on  the  same  floodplain  only  to  be  hit  by  another  torrential  rainstorm 
on  March  17,  1936.  The  amount  of  rain  water  was  compounded  by   melting 
snow  from  the  mountains.  In  the  1936  flood  no  dams  broke,  but  the 
flood  crest  took  25  lives  and  did  $41  million  in  property  damage  to  a 
city  which  had  more  than  doubled  in  size.  With  two  floods  in  50  years 
the  U.S.  Army  Corps  of  Engineers  spent  $7  million  on  a  flood  control 
project  to  make  Johnstown  the  "flood  free  city"  (Gilbert,  1977). 

On  July  19,  1977,  Johnstown  experienced  a  rainstorm  which  was 
expected  to  move  away  by  midnight  but  stopped,  releasing  an  incredible 
11.87  in  of  water  between  July  19  at  9  pm  and  4  am  on  July  20,  creating 
a  severe  flash  flood  and  breaching  five  small  earthen  dams.  The  most 
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serious  was  the  Laurel  Run  Dam,  which  was  42  ft  high  and  impounded  22 
acres,  generating  a  15-ft  wave  directly  into  to  the  Comenaugh  River. 
This  third  Johnstown  flood  took  approximately  100  lives  and  was  respon- 
sible for  $350  million  in  property  damage  (Gilbert,  1977). 

Excellent  field  documentation  of  the  1889  failure  of  the  South 
Fork  Dam  make  the  test  of  constructing  and  evaluating  a  dynamic  model 
for  the  event  feasible.  The  Engineering  Record  (1891)  and  Engineering 
News  (1889b)  provided  detailed  engineering  descriptions  of  the  failure 
of  the  South  Fork  Dam,  the  path  of  the  flood  wave,  and  the  destruction 
of  Johnstown.  In  addition,  several  authoritative  historical  accounts 
of  the  flood  have  been  published.  The  most  notable  are  a  Ph.D.  disser- 
tation by  Shappee  (1940)  and  a  book  by  McCullough  (1968). 

Reports  on  the  1936  flood  are  much  more  sketchy  even  though 
property  damage  as  measured  in  dollars  was  greater.  However,  some 
water  flow  data  were  compiled  by  the  United  States  Geological  Survey. 
The  1977  flood  has  brought  Johnstown  new  fame  and  several  rather  exten- 
sive scientific  fact  collections  are  now  available  from  the  National 
Oceanic  and  Atmospheric  Administration  (1977)  and  United  States  Geologi- 
cal Survey  (Ambruster,  1978;  Brua  and  Humphries,  1978). 

Gilliland  (1975)  modeled  the  failure  of  a  slime  pond  dam  and  the 
disordering  effect  of  the  phosphorus  spill  on  the  Peace  River  in  energy 
circuit  language. 
Guatemala  and  the  Earthquake  cf  February  4,  1976 

The  Guatemala  earthquake  of  February  4,  1976,  and  the  catastrophic 
destruction  of  £1  Progreso  was  used  as  the  field  test  of  the  theory  of 
energy  concentration,  storage,  and  catastrophic  release  triggering  sec- 
ondary catastrophic  events.  Much  of  the  supporting  fieldwork  in  this 
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research  was  done  in  the  Department  of  El  Progreso  with  a  land  area  of 
19,922  km2  which  is  mapped  in  Figures  6,  8,  and  9.  Detailed  analysis 
was  performed  on  the  capital  of  the  department,  also  named  El  Progreso, 
which  has  a  municipal  land  area  of  220  km2  and  a  central  urban  area  of 
approximately  1.5  km2. 

The  source  of  energy  in  the  Guatemalan  earthquake  was  the  2  cm  per 
year  annual  westerly  motion  of  the  North  American  Plate  of  the  earth 
crust  with  respect  to  the  Caribbean  Plate.  The  plate  movement  had 
stored  great  quantities  of  potential  energy  in  stress  built  up  in  the 
Motagua  Fault  over  the  past  two  centuries.  Figure  8c  is  a  map  of  the 
junction  of  the  North  American  Plate  and  Caribbean  Plate  of  the  tec- 
tonic system  of  the  earth  (Tarr  and  King,  1976). 

El  Progreso  experienced  almost  total  destruction  by  the  earthquake 
of  February  4,  1976,  due  to  its  relatively  close  proximity  to  the 
Motagua  Fault  rupture  and  earthquake  epicenter  and  to  the  softness  of 
the  alluvial  soil  on  which  the  adobe  structures  of  the  town  were 
constructed. 

Previous  research  on  the  Guatemala  earthquake  has  been  primarily 
concerned  with  scientific  documentation  of  the  geological  characteris- 
tics of  the  event,  its  size,  magnitude,  and  the  geographic  distribution 
of  the  destruction  generated  by  the  earthquake.  The  most  useful  scien- 
tific summary  studies  were  carried  out  by  the  Comite  de  Emergencia, 
Guatemala  (1976)  and  the  United  States  Geological  Survey  (Espinosa  et 
al . ,  1976).  The  maps  in  Figure  9  summarize  the  data  on  the  damage  to 
adobe  structure  and  the  estimated  earthquake  intensity  distribution  in 
Guatemala  (Espinosa  et  al_.,  1976) .  In  addition  to  development  of  a  re- 
construction plan  for  El  Progreso  (Alexander  and  Sipe,  1977),  many  other 


Figure  3.  Location  maps  of  Guatemala,  El  Progreso,  and  the 
Motagua  Fault. 

(a)  Map  of  Guatemala  showing  the  location  of  the 
department  of  El  Progreso  (Alexander  and 
Sipe,  1977). 

(b)  Map  of  the  Department  of  El  Progreso  showing 
the  location  of  the  municipality  in  relation 
to  the  other  seven  municipalities  (Alexander 
and  Sipe,  1977). 

(c)  Motagua  fault  zone  junction  of  Caribbean  and 
North  American  Plate.  Jordan  (1975)  modi- 
fied by  Tarr  and  King  (1976). 
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Figure  9.  Maps  of  damage  to  adobe  structures  and  earthquake 
intensity  in  Guatemala. 

(a)  Contour  map  showing  damage  to  adobe-type  struc- 
tures caused  by  the  February  4  earthquake 
(Espinosa  et  al_. ,  1976). 

(b)  Isoseismal  map  of  intensity  distribution  from 
spread  of  the  main  event.  Intensity  is  mea- 
sured by  the  Modified  Mercian i  Scale 
(Espinosa  et  al . ,  1976) . 
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foreign  governments  are  sponsoring  reconstruction  projects  throughout 
Guatemala.  The  only  other  formal  study  of  the  effect  of  the  Guatemala 
earthquake  on  El  Progreso  was  a  cadastre  mapping  project  sponsored  by 
the  Guatemalan  National  Geographic  Institute  (Giron,  1976).  The  des- 
truction was  so  great  that  identification  of  landmarks  to  identify 
property  ownership  was  a  substantial  problem.  El  Progreso  was  studied 
in  detail  by  the  author  (Alexander  and  Sipe,  1977). 


METHODS 

Description  of  Modeling  Language  and  Symbols 

The  symbols  used  in  the  systems  diagrams  were  established  by  Howard 
T.  Odum  (1971)  and  are  part  of  the  energy  circuit  language.  The  language 
combines  several  powerful  approaches  which  show  energetics   and  provide 
insight  into  the  mathematical  description  of  a  system,  while  also  illus- 
trating the  process  under  study  in  a  holistic  manner.  Energy  circuit 
language  contains  a  hierarchy  of  symbols,  which  allow  the  diagraming  of 
several  levels  of  complexity  in  one  model. 

Aggregate  symbols  are  summarized  in  Figure  10.  They  consist  cf  the 
cycling  receptor  module,  the  self-maintaining  module,  the  production 
module,  and  the  logic  module.  The  modules  are  constructed  by  combining 
members  of  the  set  of  mathematical  symbols  presented  in  Figure  11. 

Energy  circuit  models  are  normally  constructed  by  arranging  math- 
ematical symbols  within  group  symbols  to  show  systems  and  to  facilitate 
comprehension  of  the  concepts  embedded  in  the  model.  Standard  conven- 
tion also  calls  for  arranging  the  model  so  that  subsystems,  which  con- 
centrate dilute  energy,  are  at  the  left  while  energy  flows  increase  in 
quality  toward  the  right  of  the  diagram.  Heat  sinks  are  always  used 
to  show  loss  of  degraded  energy  from  the  system.  This  creative  use  of 
symbols  and  organization  provides  several  levels  of  hierarchy.  The 
highest  level  is  formed  by  arranging  the  diagram  by  energy  concentra- 
tion or  quality,  the  intermediate  by  use  of  group  symbols,  while  the 
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Figure  10.  Energy  circuit  language  group  symbols  (H.  T.  0dum,1971 ) 
Each  group  symbol  is  made  up  of  a  set  of  mathemat- 
ical symbols. 

(a)  The  cycling  receptor  has  the  condition  of  con- 
servation of  matter  that  is  ordered  by  input 
energy  flow  and  is  degraded  and  recycled  with 
output  of  energy. 

(b)  The  self-maintaining  consumer  module  processes 
energy  and  feeds  back  control  to  other  parts 
of  the  system. 

(c)  The  production  module  fixes  low  quality  energy 
and  the  control  of  a  consumer  module. 

(d)  The  logic  module  performs  discrete  operations 
based  on  control  inputs. 
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Figure  11.  Energy  circuit  language  mathematical  symbols  used 

to  construct  modules  in  Figure  8  (H.  T.  Odum,  1971). 

(a)  Outside  source  of  energy  supply  to  the  sys- 
tem controlled  by  external  system  and  not 
affected  by  energy  demand  on  the  internal 
system;  a  forcing  function  E. 

(b)  Constant  flow  source  of  energy  from  outside 
the  system. 

(c)  A  pathway  whose  flow  is  proportional  to  the 
quantity  in  the  upstream  storage;  J  =  KjE. 

(d)  Storage  of  some  quantity  in  the  system.  The 
rate  of  change  of  the  storage Q  with  time 
equals  the  inflows  minus  outflows;  Q  =  J  -  KQ. 

(e)  Interaction  of  two  flows  to  produce  an  output, 
which  is  a  function  of  the  two  flows  or  stor- 
ages driving  the  flows.  For  multiplicative 
interaction;  f(X,Y)  =  KXY. 

(f)  Transactor  symbol  for  which  money  flows  in 
one  direction  and  energy  flows  in  the  other 
with  price  (P)  adjusting  one  flow  Jx  in  pro- 
portion to  the  other;  Jj  =  PH  . 

(g)  Sensor  of  the  magnitude  of  flow,  J  symbol 
also  used  to  sense  storages  Q;  J2  =  RJ  . 

(h)  Constant  gain  amplifier  multiplies  the  in- 
flow Jx  by  a  constant  A.  Power  is  extracted 
from  energy  source  E;  J2  -  AN. 

(j)  Logic  switch  turns  flew  on  or  off  based  on 
logic  such  that  the  output  equals  the  input 
when  logic  is  one  and  zero  when  logic  is 
off;  on  J2  =  J],  off  J2  =  0. 

(k)  Logic  cemparitor  compares  A  with  B  to  provide 
a  logic  signal  C;  if  A  >  B  C  =  1;  if  A  <  B 
C  =  0. 
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the  dynamic  and  kinetic  detail  is  constructed  with  the  mathematical 
symbols. 

Development  of  Models 

The  models  in  this  dissertation  were  developed  in  a  sequential 
order  starting  with  the  basic  order-disorder  model  in  Figure  1.  A 
family  of  this  type  of  models  was  constructed  and  analyzed  mathemat- 
ically by  setting  up  each  model  in  differential  equation  form,  setting 
the  rate  of  change  equal  to  zero,  and  solving  for  the  steady-state 
solutions.  Dynamic  characteristics  were  next  tested  by  simulating  each 
model  in  the  family  of  order-disorder  models  as  depicted  by  Figure  1. 
A  preferred  order-disorder  model  was  tested  by  modeling  the  previous 
studies  on  the  effect  of  gamma  radiation  stress  of  a  tropical  rain 
forest  at  El  Verde,  Puerto  Rico. 

A  catastrophic  consumer  pathway  was  added  to  the  basic  order- 
disorder  as  shown  in  Figure  2.  The  family  of  similar  models  was  ana- 
lyzed and  simulated  to  test  the  models'  dynamic  behavior.  The  model  of 
catastrophic  energy  release  was  tested  by  simulating  the  failure  of  the 
South  Fork  Dam  (which  generated  the  Johnstown  flood  of  1889)  and  the 
Guatemala  earthquake  of  1976. 

Next,  two  order-disorder  models  were  connected  in  cascade  and 
tested  by  simulation  of  the  Guatemala  earthquake  and  the  secondary  des- 
truction of  El  Progreso  and  the  Johnstown  flood  and  the  secondary  des- 
truction of  Johnstown.  Figure  3  describes  the  basic  configuration.  By 
connecting  units  in  a  web,  a  more  general  model  was  developed  as  a  meth- 
odological guide  for  calculation  of  a  power  spectrum  for  global  disasters 
and  for  estimation  of  quantity  of  solar  energy  embodied  in  generic  types 


49 


of  natural  hazards  and  their  disruptive  pulses  of  energy.  The  model 
summarized  in  Figure  4  is  a  global  energy  budget  model  that  relates 
the  atmospheric,  oceanic,  biological,  geological,  and  urban  processes 
responsible  for  concentration  and  storage  of  solar  energy  to  power 
storms,  floods,  earthquakes,  fires,  and  volcanoes. 

Limiting  Factors  and  Differential  Equations 

The  energy  circuit  models  developed  in  this  research  incorporate 
internal  and  external  limiting  factors.  Internal  limits  result  from 
the  shortage  of  disordered  matter  in  the  cycling  modules.  External 
limits  are  provided  by  energy  flow  sources  such  as  solar  radiation. 
Both  limits  to  growth  are  incorporated  in  diagrammatic  and  differential 
equation  form  of  the  models.  Figure  12  is  a  typical  example,  illus- 
trating the  use  of  internal  and  external  limits  to  growth  in  the  same 
order-disorder  model.  The  differential  equation  for  the  energy  storage 
term  Ql   is  given  below: 

Qi  =  Jr  Kj  Q2  -  K2  Qj  (7) 

The  energy  flow  J  is  split  between  J0  and  J  as  in  equation  8: 


J  =  J0  +  Jr-  (3) 

The  energy  consumption  of  the  system  is  given  by  equation  9: 

Jo  =  K0  Jr  Q2-  (9) 

Substituting  equation  9  into  equation  8  and  solving  for  Jp  yields  equa- 
tion 10: 


Figure  12.  Order-disorder  model  in  energy  circuit  language, 
illustrating  the  incorporation  of  internal  and 
external  limiting  factors.  The  internal  limit 
is  provided  by  conservation  of  matter  in  the 
order-disorder  loop.  The  external  limits  to 
growth  are  a  result  of  the  use  of  a  constant 
flow  energy  source. 
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Substituting  equation  10  into  equation  1  yields  the  desired  results  in 
equation  11 : 

Qi  =  k0  o2  +  i    -K^i-  <") 

Rearranging  equation  9  defines  Kq 

J. 


Ko 


J7QT  (12) 


Requiring  conservation  of  mass  such  that  the  total  mass  in  the 
storages  Qi  and  Q2  to  constant  C  yield  equation  13: 

C  =  K3  Qj  +  Q2,  (13) 

where  K  is  the  mass  per  unit  energy  for  the  system  under  study. 
Solving  for  Q  yields  equation  14: 

Q2  =  C  -  K3  Q1.  (14) 

Substituting  equation  14  into  equation  11  to  remove  Q2  yields  equa- 
tion 15: 

KiJC  -  KiJKgQx 


1  +  K0C  -  KgKsQ: 


K2Qi.         (15) 


Equation  15  is  a  differential  equation  describing  the  behavior  of 
the  model  in  Figure  12  and  incorporating  the  internal  limit  imposed  by 
conservation  of  matter  and  the  external  limit  of  a  constant  flow  energy 
source. 

Data  Assembly  and  Evaluation  of  Models 

The  energy  circuit  models  in  Figures  1  to  4  were  used  as  a  guide 
for  collection  and  organization  of  data.  Data  were  assembled  from 
literature  for  evaluation  of  global  models.  A  table  was  constructed  for 
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each  energy  flow  and  storage  in  the  model.  The  table  contains  the  energy 
value  of  the  flow  or  storage  with  each  calculation  documented  by  a  note 
in  Appendix  I.  Each  data  source  was  referenced  to  the  appropriate  re- 
search and,  when  possible,  several  measures  or  estimates  were  used  for 
each  reference  value. 

Simulation  Procedures 

A  set  of  simultaneous  nonlinear  differential  equations  was  developed 
for  each  of  the  energy  circuits  models  of  the  case  study  disasters.  Ini- 
tial conditions  and  coefficients  for  each  model  were  calculated  from  the 
table  of  energy  flows  and  storages  that  was  constructed  to  document  the 
model.  Mathematical  solutions  were  used  to  develop  the  behavioral  prop- 
erties of  the  basic  order-disorder  model  used  as  a  basis  for  this  re- 
search. However,  this  was  not  feasible  for  more  complex  models. 

Numerical  approximation  techniques  and  solution  by  analogous  system 
were  used  to  study  and  simulate  the  dynamic  behavior  of  the  disaster 
models.  Numerical  approximation  was  the  most  straightforward  simulation 
method.  During  early  stages  of  model  development,  simulations  were  tested 
on  a  Hewlett  Packard  25  programmable  calculator.  However,  more  extensive 
computer  runs  were  performed  with  the  simulation  language  DYNAMO  (Pugh, 
1970).  Working  with  DYNAMO  offered  the  advantages  of  a  large,  dynamic 
range  of  modern  digital  computers.  Thus,  the  models  could  be  simulated 
without  going  through  a  scaling  procedure  necessary  with  analog  simula- 
tion. DYNAMO  models  were  developed  directly  from  the  energy  circuit  dia- 
gram so  that  one  level  equation  was  written  for  each  state  variable 
(storage  tank),  while  one  rate  equation  was  written  for  each  energy  flow 
pathway.  Digital  computer  programs  such  as  DYNAMO,  simulate  nonlinear 
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differential  equation  solutions  by  mathematical  approximation.  To  insure 
reasonable  accuracy  with  Euler  integration  technique,  a  time  step  of  0.1 
percent  of  the  disaster  duration  was  used  as  a  maximum.  Thus,  many  com- 
putations are  required  for  one  simulation  run,  rendering  the  digital 
much  slower  than  analog  computer  methods  that  continuously  solve  the  set 
of  nonlinear  equations.  A  Tektronix  digital  computer  terminal  and  time- 
sharing techniques  were  used  to  carry  out  model  development  and  simula- 
tion. Some  of  the  models  were  scaled  and  wired  on  the  analog  computer 
to  provide  hands-on  immediate  response.  Analog  simulation  was  used  for 
sensitivity  analysis  to  study  relative  response  of  the  system  to  each 
parameter.  Information  gained  from  analog  simulation  was  used  to  plan 
runs  on  the  DYNAMO  version  of  the  model.  The  data  for  storages  and 
flows  that  were  determined  to  be  sensitive  by  analog  or  digital  simula- 
tion were  reverified  and  checked  for  accuracy  with  a  second  or  third 
reference,  and  when  possible  by  several  independent  techniques  of  cal- 
culating the  energy  flows  and  storages  used  in  evaluating  the  models. 

Calculation  of  Ratios 

In  addition  to  evaluation  and  simulation,  the  energy  circuit  models 
in  this  dissertation  were  used  as  systematic  guides  to  energy  analysis 
of  energy  flows.  Several  types  of  ratios  of  energy  flows  were  calcu- 
lated to  add  to  the  quantitative  understanding  of  natural  hazards  and 
their  disruptive  effects.  As  explained  in  this  section,  three  important 
ratios  calculated  are  energy  quality,  embodied  energy,  order-disorder 
ratios,  and  disaster  amplifier  ratios. 

Energy  quality  factors  were  calculated  from  the  evaluated  biosphere 
web  as  diagramed  in  Figure  4  (Odum,  1978).  Under  the  web  concept  all  of 
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the  energy  flows  are  interrelated  in  such  a  way  that  each  is  produced 
as  a  by-product  of  all  the  others.  Taking  this  point  of  view,  the 
total  incident  solar  energy  to  the  world  is  required  to  produce  all 
energy  flows  and  build  all  the  structures.  In  the  example  of  global 
net  primary  production,  rain,  wind,  dry  air,  and  nutrients  of  the  soil 
are  all  known  to  be  necessary  energy  flows  as  direct  radiant  solar 
energy.  Once  the  global  energy  web  is  evaluated,  the  energy  flow  in  a 
specific  process  is  divided  by  the  total  solar  energy  necessary  to 
power  the  earth.  For  this  research  the  mean  solar  incident  radiation 
was  used  as  the  basis  for  all  energy  quality  calculation  except  for 
electromagnetic  radiation. 

Embodied  energy  factors  were  calculated  by  multiplying  an  energy 
flow  in  a  specific  case  by  the  global  quality  factor,  which  was  devel- 
oped from  the  energy  web  model  in  Figure  4. 

The  order-disorder  model  diagramed  in  Figure  13  illustrates  the 
energy  flows  used  in  calculating  the  ratio  of  the  energy  required  to 
produce  order  to  the  energy  required  to  produce  disorder  in  the  system. 
The  ratio  was  calculated  by  dividing  flow  A  by  flow  C  in  Figure  13. 
When  the  ratio  of  the  storage  of  order,  which  is  recycled  to  disorder 
(flow  B),  to  the  flow  of  energy  required  to  pump  the  recycling  process 
is  calculated,  flow  B  is  then  divided  by  flow  C.  For  purposes  of 
clarity,  the  first  is  referred  to  as  the  order-disorder  ratio,  while 
the  second  is  designated  as  the  disaster-amplifier  ratio. 

Calculation  of  Spatial  Distribution  of  Seismic  Energy 

A  method  is  proposed  for  calculation  of  the  spatial  distribution 
of  seismic  energy  from  observed  field  intensities  measured  on  the 


Figure  13.  Order-disorder  model  illustrating  energy  flows  used 

in  calculating  the  order-disorder  ratio  and  the 

disaster-amplifier  ratio. 

Order-disorder  ratio 

Energy  required  to  produce  order    _  A 
Energy  required  to  produce  disorder    C 

Disaster  amplifier  ratio 

Storage  of  order  disordered  externally  _  A 
Energy  required  to  produce  disorder     C 
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Modified  Mercalli  scale.  Intensity  is  a  subjective  numerical  index  de- 
scribing the  effect  of  an  earthquake  on  man,  on  structure,  and  on  the 
surface  of  the  earth.  The  Modified  Mercalli  scale  of  1931  has  intensity 
rankings  from  I  to  XII  (see  legend  on  Figure  30b  for  definition  of  scale; 
Espinosa  et  al_. ,  1976;  Richter,  1958).  Acceleration  A  was  calculated 
from  the  intensities  discussed  above  using  equation  16  (Gutenberg  and 
Richter,  1969). 

log  a  =  1/3  -  1/2  (16) 

If  energy  was  distributed  uniformly,  one  could  find  the  dissipation 
per  unit  of  land  area  by  dividing  the  total  energy  by  the  total  land 
area  affected;  however,  this  is  not  a  realistic  assumption.  For  the 
calculations,  it  was  assumed  that  energy  was  uniformly  distributed  in 
each  isoseismal  zone.  Equation  17  mathematically  gives  this  relation- 
ship: 

»1-ES1/A1-  <17> 

where  ei  is  the  energy  dissipated  per  unit  area  of  zone,  i,  E  .  is  the 
total  seismic  energy  dissipated  in  zone  i,  and  A.  is  the  land  area  of 
zone  i.  Energy  was  related  to  acceleration  by  combining  the  following 
two  equations  (Gutenberg  and  Richter,  1969): 

Esi  =  c  t  (A^/T)2  (18) 

where  c  is  a  transmission  constant,  t  is  the  duration  of  the  wave  group, 
and  T  is  the  period  of  vibration,  all  of  which  were  constant  for  the 
case  under  study.  Ai  is  the  amplitude  of  the  ground  motion  for  zone  i. 
Amplitude  A.,  was  related  to  acceleration  a,  as  shown  in  equation  19 
(Gutenberg  and  Richter,  1969): 
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ai  =  (A./T)  4tt2/T  (19) 

and  then  substitute  equation  19  into  equation  18,  combining  constants 
to  form  a  new  constant  c,  which  yields  equation  20. 

e.  -  c  a.2.  (20) 

The  constant  of  proportionality  c  was  then  evaluated  by  substituting 
equation  17  into  equation  20  to  get  equation  21. 

c  =  E  /("   A.a  2)  (21) 

Substituting  equation  16  into  equation  21  and  rearranging,  yields  the 
desired  relationship  between  energy  distribution  and  earthquake 
intensity: 


E  -  (I   C  A-0'66  I_1).  (22) 

i=l 


Plan  of  Study 

The  theory  of  cycles  of  order  and  disorder  was  first  developed  by 
analyzing  several  possible  configurations  of  order-disorder  models.  The 
analysis  was  carried  out  by  calculation  of  the  steady-state  solution  to 
the  models  in  differential  equation  form.  The  dynamic  properties  of  the 
model  were  tested  using  simulation  techniques.  As  a  refinement  of  the 
order-disorder  concept,  a  consumer  module  was  added  to  the  order-disorder 
model.  The  catastrophic  energy  release  properties  of  the  model  were  ini- 
tially analyzed  for  steady-state  properties  of  the  differential  equation 
form  of  the  family  of  models.  Dynamic  characteristics  were  analyzed 
using  simulation  techniques. 
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An  energy  basis  for  disasters  was  developed  by  connecting  order- 
disorder  models  to  form  a  global  energy  web  model  with  pathways  for 
cascading  of  catastrophic  release  processes.  The  resulting  model  was 
then  used  to  evaluate  the  energy  quality  of  natural  disasters  and  to 
illustrate  the  control  action  of  their  pulsing.  The  first  section  of 
the  results  develops,  evaluates,  and  documents  this  model  by  diagraming 
the  atmospheric,  oceanic,  biological,  geological,  and  urban  systems  of 
the  world,  their  energy  flows,  and  their  interconnection.  Then  energy 
ratios  were  calculated  from  energy  flows  in  the  evaluated  models.  Quan- 
titative comparison  of  order-disorder  ratios  were  made  with  all  of  the 
energy  flows  involved  in  a  common  unit  of  measure  as  calories  of  solar 
equivalent  energy. 

To  verify  the  theory  of  cycles  of  order  and  disorder,  three  disas- 
ters were  studied  with  models,  energy  evaluation,  and  simulation  data 
from  an  irradiation  disordering  experiment  in  a  tropical  rain  forest 
was  used  to  evaluate  a  pulsing  order-disorder  system.  The  model  was 
simulated  and  compared  to  field  observations.  An  order-disorder  model 
containing  a  catastrophic  consumer  module  was  evaluated  for  the  failure 
of  the  dam,  which  created  the  major  disaster  at  Johnstown.  Two  models 
were  connected,  so  their  effect  could  be  cascaded  and  evaluated  for  a 
major  earthquake  and  the  secondary  disruption  of  several  cities.  The 
flood  and  earthquake  were  simulated  and  the  results  were  compared  to 
field  observation.  Finally  ratio  of  ordering  and  disordering  energies 
were  assembled  to  consider  theories  of  energy  flow  in  disasters. 


RESULTS 

Simulation  and  Analysis  of  Order-Disorder  Models 

Results  of  the  study  of  several  possible  order-disorder  models  are 
presented  in  this  section,  starting  with  a  simple  model  and  progressing 
to  more  complex  examples.  Differential  equations  to  describe  the  behav- 
ior of  each  are  also  presented  along  with  time  simulations  of  each  model 
All  of  these  models  have  the  common  characteristic  of  materials  cycling 
as  a  mechanism  to  capture,  concentrate,  and  store  energy  and  are  mathe- 
matically constrained  to  conserve  matter.  Figure  14a  is  the  basic 
order-disorder  model  which  was  first  applied  in  enzyme  reactions 
(Michael ius-Menton,  1913).  The  behavior  of  the  level  of  order  Q1     is 
described  by  differential  equation  22. 

Qi  =  M1Q2  "  Mi  (22) 

Requiring  conservation  of  matter  yields  equation  23. 

Q2  =  QT  -  Ql  (23) 

Substituting  equation  23  in  equation  22  yields  a  single  differential 
equation  which  describes  the  model  behavior  of  equation  24. 

Qj  =  k1E(QT  -  QJ  -  K2Q!  (24) 

To  find  the  steady-state  solution,   Qx  was  set  equal  to  zero  and 
equation  24  was  solved  for  Q,  ,  yielding  equation  25. 
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3!SS  "  (klET^V  «T  (25) 

Assuming  E  constant  and  substituting   kj  for  k:E  in  equation  25, 
produces  equation  26.  Let  QT  =  1 ,  K3  =  0.6,  K2  =  0.2,  and  Ei  =  1 


=  0.75.        (26) 


wiss   k  +  k2  yT   0.6  +  0.2 

This  form  of  the  solution  appears  often  in  problems  involving  dis- 
tribution of  a  resource  between  two  states.  Figure  14a  also  shows  the 
result  of  a  DYNAMO  (Pugh,  1970)  simulation  of  this  order-disorder  model. 
Note  that  the  steady-state  value  determined  both  mathematically  and  by 
simulation  is  the  same. 

Similar  analysis  and  simulation  were  performed  for  each  model  in 
Figures  14b  through  h.  The  analyses  are  summarized  in  Table  1,  while 
the  simulation  results  of  each  model  with  the  same  pathway  coefficients 
and  energy  source  values  were  shown  parallel  to  the  energy  circuit  model 
diagram.  The  documentation  of  each  simulation  model  is  contained  in 
Appendix  II. 

The  series  of  models  presented  in  this  section  all  used  only  one 
differential  equation  in  the  simulation  model.  The  second  equation 
dropped  out  when  the  constraint  of  holding  total  matter  constant  was 
applied. 

The  oraer-disorder  model  in  Figure  14b  incorporates  the  addition 
of  a  feedback  of  some  of  the  order,  Ql   to  pump  more  energy.  The  con- 
servation of  matter  imposed  gives  this  model  logistic  characteristics 
as  seen  in  the  simulation  results.  This  model  has  two  steady-state 
solutions.  Figure  14c  incorporates  an  external  energy  source,  E2  to 
pump  the  recycle  process.  The  simulation  results  are  the  same  as 


Figure  14.  A  comparison  of  several  order-disorder  models  and 
their  simulation  results.  In  all  simulations  the 
following  values  were  used: 

Ei  =  1,  E2  =  1,  Ki  =  0.6,  K2  =  0.2. 
Also  conservation  of  matter  was  imposed  on  each 
model . 

(a)  Simple  order-disorder  model  and  simulation 
results  with  one  energy  source  Ei. 

(b)  Order-disorder  model  and  simulation  results 
with  autocatalytic  feedback  and  one  energy 
source  E2. 

(c)  Order-disorder  model  and  simulation  results 
with  an  ordering  energy  source  Ei  and  a 
disordering  energy  source  E2. 
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DME 


Figure  14.  Continued 

(d)  Autocatalytic  order-disorder  model  and  simula- 
tion results  ordering  energy  source  Ei  and 
disordering  energy  source  Ez- 

(e)  Autocatalytic  order-disorder  model  and  simula- 
tion results  with  cooperation  feedback  in  order 
loop  and  cubic  drain  on  disorder  loop. 

(f)  Autocatalytic  order-disorder  model  and  simula- 
tion results  with  square  drain  in  disorder 
loop  and  an  disordering  energy  source. 
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Figure  14.  Continued 

(g)  Dual  autocatalytic  order-disorder  model  and 
simulation  results  which  behaves  like  a 
flip-flop. 

(h)  Autocatalytic  order-disorder  model  and  sim- 
ulation results  with  constant  flow  energy 
source. 

(i)  Simple  order-disorder  model  simulation  re- 
sults with  two  differential  equations. 
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TABLE  1 .  STEADY-STATE  SOLUTIONS  AND  STEADY-STATE  VALUES  FOR  THE  ORDER- 
DISORDER  MINIMODELS  IN  FIGURE  14. 

Model  No.     Steady-state  Solution  Steady-state  Value 


Qiss  ~  kxE  +  k2 


yiss   kiEi  +  koE- 


■l11!  T  N2C2 


0.75 


ki  EiQT  "  K? 
Qiss=     k^   2°r0  0'67or0 


0.75 


_  Ki  Ei^t  "  k2E2or  0  0 . 67  or  0 

^iss      kjEi 


n      Ki  Ei  qt 


EXQT 

or  0  0 .75  or  0 


'iss   Ex  +  k2E2 


Qiss  =  *  "  .  %2or  0  1  to  0 

-1  ~  ^l 
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Figure  14a,  thus  this  is  another  form  of  the  Michael ius-Menton  kinetics. 
Figure  14d  is  the  logistic  model  with  the  addition  of  an  external  energy 
source  to  pump  the  recycle  process.  Figure  14e  illustrates  the  effect 
of  using  the  square  of  the  order  0^  to  pump  order  and  disorder  with  two 
external  energy  sources.  Figure  14f  is  another  configuration  of  the 
order-disorder  model  which  has  logistic  behavior  and  also  uses  dual 
energy  sources  for  order  and  disorder  processes.  Figure  14g  is  an  un- 
usual order-disorder  model,  which  behaves  as  a  flip-flop  circuit.  This 
may  be  unrealistic  for  there  is  no  depreciation  on  the  storages  except 
through  the  dual  production  functions  of  order  and  disorder.  Figure  14h 
incorporates  internal  and  external  limits  to  form  a  logistic  behavior. 
The  external  limit  is  provided  by  a  constant  flow  energy  source  and  the 
internal  limit  provided  by  material  cycling. 

To  determine  that  this  mathematical  manipulation  would  not  affect 
the  simulation  result,  the  model  in  Figure  14i  was  simulated.  The 
model  diagramed  in  Figure  14i  is  the  same  as  Figure  14a.  This  simula- 
tion was  performed  using  differential  equations  27  and  28. 

Qi  =   KjEiQa  -  K2Qi  (27) 

Q2  =  K2Qi  -  K3E!Q2  (28) 

To  insure  stability  in  the  model,  (Kx)  was  set  equal  to  (K3)  as  in  a 
Lotka  Loop.  Figure  14i  also  shows  the  result  of  the  simulation,  and 
it  was  identical  to  the  results  obtained  in  Figure  14q.  Thus  identical 
simulation  results  were  obtained  for  the  order-disorder  models  with  two 
differential  equations  in  Figure  14i  when  the  coefficients  were  balanced 
as  with  one  differential  equation  and  the  conservation  of  matter  con- 
straint in  Figure  14q. 


71 


Switching  Pulse  Models 

A  previous  energy  circuit  model  used  a  threshold  approach  to  modeling 
disasters  or  disorder  in  processes,  such  as  fire  or  epidemic.  A  switch 
controlled  by  a  logic  threshold  was  used  to  model  the  disruption  process 
(Odum,  1974).  Two  logic  configurations  of  the  model  are  given  in  Figure 
15.  The  behavior  of  the  model  is  simple  and  adequate  for  a  first  order 
approximation  to  modeling  catastrophic  events.  When  the  level  in  the 
storage  Qj  exceeds  a  present  threshhold  T,  the  disaster  pathway  K3  is 
logically  switched  on  and  remains  conductive  until  the  energy  storage  Qj 
is  depleted  below  its  self-extinguishing  threshhold,  T2.  Also,  the  model 
may  be  triggered  by  applying  a  logic  pulse  at  P.  The  limitations  of  this 
model  are  apparent  when  the  simulation  is  examined  over  the  period  of  the 
disaster.  However,  this  model  does  not  explain  the  dynamic  behavior  of 
the  diaster.  It  simply  replicates  observed  results'  destruction. 
Exponential  Surge  Model 

To  facilitate  the  simulation  of  actual  disaster  processes  such  as 
fire,  earthquakes,  and  floods,  a  dynamic  model  was  developed  capable 
of  explaining  energy  concentration,  storage,  and  catastrophic  release 
in  the  disaster  process.  As  a  minimum,  the  model  needed  the  long-term 
characters  of  the  logically  switched  one  in  Figure  15.  Previous  work 
by  other  researchers  with  fire  has  shown  the  energy  in  the  fire  to 
increase  exponentially  until  the  source  of  fuel  is  consumed  by  the  fire 
(Byram  et  al . ,  1966).  These  previous  findings  suggested  that  an  exponen- 
tial growth  model  would  be  appropriate  for  modeling  the  fire.  Figure  16q 
is  an  energy  circuit  order-disorder  model  with  the  logic  pathway  replaced 
by  an  autocatalytic  consumer,  fire.  In  this  model,  the  fire  was  triggered 
by  an  energy  pulse  applied  to  the  fire  energy  storage  Q3.  Once  initiated, 


Figure  15.  Energy  circuit  disorder  model  of  disaster  typical  of 
previous  fire  and  epidemic  models  of  other  re- 
searchers. Disorder  process  is  logically  switched. 

(a)  Model  with  external  trigger  required. 

(b)  Model  with  internal  threshold  trigger 
mechanism. 
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Figure  16.  Order-disorder  energy  circuit  model  with  catastrophic 
energy,  pulse  generated  by  consumer  module. 

(a)  Model  with  external  control  pulse  required. 

(b)  Model  with  internal  threshold  for  pulse 
generation. 
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the  fire  grew  exponentially  until  the  fuel  source  was  consumed  to  the 
point  of  self-extinction.  However,  this  model  could  not  self-start. 

To  incorporate  the  self-starting  characteristic  into  the  model, 
the  linear  drain  pathway,  which  represents  depreciation,  was  fed  into 
the  disaster  energy  storage  Q3.  As  long  as  the  drain  or  depreciation 
of  Q3  is  greater  than  the  inflows,  the  model  was  stable,  and  the  order 
or  potential  energy  in  Ql   could  be  increased.  However,  if  0^  became 
large  enough  so  that  the  depreciation  of  Qj  exceeded  the  losses  from 
Q3,  the  autocatalytic  fire  model  reached  a  critical  point  where  it 
could  self-start  and  grow  until  the  energy  source  Q}  was  consumed. 
This  model  was  also  capable  of  periodic  oscillation  as  shown  in  Figure 
16.  The  storage  of  energy  grows  to  a  threshold  and  periodically  dis- 
charges. 

This  model  has  several  advantages  over  the  logic  model  in  addition 
to  replication  of  upper  and  lower  thresholds.  Because  it  actually 
models  the  disaster,  if  preventive  measures  can  be  taken  in  time,  the 
disaster  may  be  extinguished.  For  example,  the  fire  may  be  controlled 
or  extinguished  by  cooling  the  disaster  energy  storage  Q3  as  is  done  in 
fire  fighting  with  water.  It  is  also  possible  to  control  the  disaster 
growth  and  extinguishing  rate  based  on  the  characteristics  of  the 
stored  fuel . 

The  order-disorder  model  in  Figure  16b  had  good  dynamic  character- 
istic for  disaster  model  used  in  this  research.  The  model  conforms  to 
the  convention  of  conservation  of  matter  such  that  the  sum  of  storages 
are  constant  QT  as  in  equation  29. 

Qi  +  Q2  =  nr  (29) 

The  level  of  order  or  stress  in  the  disaster  Qx  was  described  by  equa- 
tion 30. 


Figure  17.  Simulation  of  the  periodic  buildup  and  pulse  genera- 
tion property  of  the  exponential  surge  model  in 
Figure  16b. 
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Qx  =  k!EQ2  -  k2  (Qx  +  QiQg).  (30) 

Catastrophic  energy  release  was  controlled  by  the  autocatalytic  energy 
release  pathway  k2(Qi  +  Q1Q3).  When  Q3  is  small,  this  term  is  dominated 
by  the  linear  depreciation  k2Qi-  The  value  of  energy  storage  in  the 
decay  pathway  Q3  is  described  by  equation  31. 

Q3  =  (k„  -  k3)Q1Q3  -  k5Q3  +  KeQi  (31) 

Equations  30  and  31  may  be  used  to  solve  for  steady-state  values  of  Qj 
and  Q3  by  setting  the  0^  and  Q3  equal  to  zero,  yielding  equations  32 
and  33  . 

_      kiEQT  (32) 

giss   kxE  -  kj  -  kiQ3(ss) 

.   _       Mi(ss)  (33) 

^iss   kaQ^ss)  +  k5  -  k^Q^ss) 

Also,  threshold  for  catastrophic  release  is  given  by  equation  34: 

k5Q3 


Qlt  =  (kg  "  kj  (Q3  +  k6) 


(34) 


that  is,  if  Ql  >   Qjt.  the  model  will  catastrophically  switch  to  the 
disorder  state. 

The  decay  path  of  the  disaster  is  controlled  by  the  value  of  k5, 
which  is  the  time  constant  of  the  catastrophic  decay. 

Evaluation  of  Global  Energy  Web  and  Disasters 

Figure  18  is  an  energy  circuit  model  that  attempts  to  synthesize 
concepts  about  the  earth  using  concepts  from  geology,  meteorology, 
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oceanography,  economics,  and  ecology.  It  shows  the  global  structural 
mechanisms  that  concentrate  energy  in  space  and  time,  store  energy,  and 
eventually  release  the  energy  in  catastrophic  pulses  that  seem  natural 
disasters  to  people. 

Studies  of  energy  flow  and  concentration  processes  in  the  biosphere 
are  usually  constrained  by  the  boundaries  of  the  discipline.  Parts  of 
the  global  energy  model  developed  in  this  section  have  been  studied  by 
meteorologists,  oceanographers,  ecologists,  geologists,  and  economists, 
with  each  group  concentrating  on  kinetics,  energy  flows,  and  catastrophic 
phenomena  in  their  respective  fields.  However,  energy  is  known  to  flow 
across  the  areas  covered  by  disciplines  in  the  real  world.  A  catastrophic 
event  in  one  field  often  triggers  secondary  disruptions  in  many  fields. 

Modeling  the  concentration,  storage,  and  catastrophic  release 
mechanism  produced  fundamental  insight  into  the  behavior  of  disasters. 
Conceptually,  each  natural  disaster  occurrence  is  viewed  as  a  disrup- 
tion of  a  storage  following  a  peak  in  the  normal  global  energy  cycle. 
Earthquakes,  volcanoes,  tidal  waves,  rock  slides,  and  other  disruptive 
geological  hazards  are  viewed  as  momentary  nonlinear  pulses  of  energy 
generated  by  bottlenecks,  which  produce  a  temporary  concentration  of 
potential  energy.  Floods,  lightning,  hurricanes,  and  tornados  result 
from  a  rapid  release  of  excess  potential  energy  stored  in  atmospheric 
systems.  Forest  fire  is  the  result  of  catastrophic  release  of  chemical 
potential  energy  stored  by  biological  systems  in  organic  matter.  Urban 
fire  and  war  can  be  modeled  as  releases  from  storages  of  society,  often 
disordered  by   earthquakes,  floods,  storms,  tsunamies,  and  volcanoes. 

The  world  model  suggests  interrelationships  of  natural  disasters 
and  the  energy  trajectory  required  for  them.  It  allows  energy  flows 
from  various  energy  sources  to  be  compared. 
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Rational  quantitative  comparison  of  sectors,  disaster  phenomena, 
and  ratios  of  flow  can  be  made  if  all  of  the  energy  flows  involved  are 
put  in  a  common  unit  of  measure  such  as  heat  equivalents,  Calories  of 
embodied  solar  energy,  and  energy  quality  factors.  These  were  evaluated 
for  the  world  model  and  are  given  by  sector  next. 
Atmospheric  Circulation  and  Storms 

Figure  19  is  the  atmospheric  sector  of  the  global  natural  hazard 
model  in  Figure  18,  which  was  the  web  used  for  calculation  of  the  solar 
energy  quality  of  rain,  wind,  and  storms.  To  calculate  the  energy  qual- 
ity of  storms,  it  was  first  necessary  to  evaluate  the  model  from  the 
available  global  energy  data.  Figure  19a  presents  the  model  in  concep- 
tual form,  while  19b  is  in  mathematical  form  with  differential  equations. 
Table  2  documents  the  energy  flow  and  storage  evaluation  of  the  world's 
atmospheric  system  with  detailed  calculations  provided  in  Appendix  I. 

Storms,  such  as  cyclones,  hurricanes,  tornados,  and  thunderstorms, 
are  powered  by  warm  moist  air  conducted  toward  polar  latitudes  of  the 
earth  as  excess  heat  from  tropical  latitudes.  Storms  are  formed  by  the 
atmospheric  systems  acting  as  a  mechanism  to  release  excessive  quanti- 
ties of  potential  energy  that  could  not  be  released  in  a  competitive  way 
by  simple  diffusion.  The  release  of  excess  energy  built  up  across  fronts 
produces  turbulence.  The  wind  pumps  more  energy  through  the  atmospheric 
system.  The  turbulence  also  provides  a  positive  feedback  to  the  ocean 
by  producing  waves  which  in  turn  supply  more  water  vapor  energy  to  power 
the  storm  and  to  provide  rain  for  the  land  systems  of  the  earth. 

The  advent  of  the  satellite  made  direct  measurement  of  global  para- 
meters possible.  Measurements  from  space  platforms  essentially  confirm 


Figure  19.  Energy  circuit  diagrams  of  atmospheric  system. 

(a)  Conceptual  drawing  of  energy  circuit  model 
of  atmosphere  used  to  illustrate  energy 
quality  calculation  for  weather  and  wind. 
Energy  flows  are  in  Cal  m"2d_1  while  stor- 
ages are  in  Cal  m~2. 

(b)  Atmospheric  sector  of  global  energy  web 
model  in  Figure  18  with  differential  equa- 
tions. This  energy  circuit  cycling  recep- 
tor module  was  evaluated  for  the  major 
flows  of  energy  in  the  atmosphere.  Table 

1  provides  documentation  for  the  model. 
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the  Smithsonian  mean  value  of  the  solar  radiation  to  be  1.94  g-Cal  cm 

-1  -2  -1 

min   or  27,940  Cal  m  d  (Byers,  1974).  The  world  model  shows  the 
relation.  Solar  flares,  which  have  been  shown  to  have  an  11-year 
cycle  (Eddy,  1976),  release  high  quality  cosmic  radiation  which  is 
responsible  for  a  2  percent  fluctuation  around  this  mean  (Schneider 
and  Mass,  1975),  sometimes  called  the  "solar  constant."  The  solar 
radiation  incident  to  the  earth  was  taken  to  be  7,000  Cal  m"2d~1. 
This  was  calculated  by  dividing  the  solar  mean  by  the  ratio  of  the 
area  of  the  earth's  plane  to  the  area  of  the  earth's  surface. 

Working  with  5  years  of  data  from  satellites  (1962-1966),  Vonder 
Haar  and  Suomi  (1971)  found  the  global  annual  average  solar  radiation 
to  be  7,000  Cal  nf2d_1,  and  the  radiation  absorbed  to  be  4,900  Cal 
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m  d  .  Radiation  hitting  the  earth  was  either  absorbed  or  reflected 
by  the  clouds,  land  surface,  and  plants.  The  difference  between  the 
energy  absorbed  and  the  incident  energy,  known  as  albedo,  was  found  to 
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be  2,100  Cal  m  d  .  The  global  average  infrared  radiation  was  also 
measured  by  Vonder  Haar  and  Suomi  (1971)  from  satellite  data  and  found 
to  be  4,900  Cal  m~2d_1.  Houghton  (1954)  found  23  percent  of  the  solar 
energy  incident  to  the  globe  was  directly  consumed  in  powering  the  hydro- 
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logical  cycle  to  1,610  Cal  m  d  ,  while  the  majority  of  the  remaining 
energy,  3,290  Cal  m"  d"  ,  was  used  in  direct  heating  land  and  sea  and 
powering  the  photosynthetic  cycle.  Indirect  forms  of  solar  energy, 
however,  are  equally  important  to  the  operation  of  the  earth's  energy 
web.  The  energy  flows  evaluated  above  were  global  averages;  however, 
the  solar  radiation  was  shown  by  Houghton  (1954)  to  be  distributed  in 
such  a  way  that  the  incident  radiation  per  square  meter  of  surface  area 
drops  as  latitude  increases.  Figure  20a  is  an  early  plot  of  this 


Figure  20.  Global  distribution  of  solar  energy  with  latitude. 

(a)  Distribution  with  latitude  of  absorbed  and 
outgoing  radiation  as  computed  by  Houghton; 
solar  radiation  absorbed  by  earth  and  atmo- 
spheres; long-wave  radiation  leaving  the 
atmosphere.  The  abscissa  is  scaled  by  the 
cosine  of  latitude  (Houghton,  1954). 

(b)  Mean  meridional  profiles  (averages  within 
latitude  zones)  of  components  of  the  earth's 
radiation  budget  measured  during  the  period 
1962-66.  The  abscissa  is  scaled  by  the 
cosine  of  latitude  (Vonder  Haar  and  Soumi , 
1971). 

(c)  Total  poleward  energy  transport  required  by 
the  radiation  budget  (RT)  and  its  major  com- 
ponents.  The  free  energy  transport  by  ocean 
currents  is  AF,  LaC  is  the  transport  of  energy 
in  latent  form  of  water  vapor  and  AC,  the 
sensible  heat  plus  potential  energy  transported 
by  the  atmosphere  (Vonder  Haar  and  Suomi ,  1971). 
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phenomenon  developed  by  Houghton  (1954).  Note  that  there  was  a  net  flux 
of  energy  toward  the  poles  necessary  to  offset  the  cooling  processes  at 
the  polar  region.  Houghton's  findings  indicate  that  below  40°  latitude 
heating  exceeded  cooling,  and  above  40°  the  reverse  was  true.  A  more 
accurate  estimation  of  the  relationship  was  available  from  satellite 
data  and  was  plotted  in  Figure  20b  (Vonder  Haar  and  Suomi ,  1971).  The 
planetary  albedo  was  also  plotted  on  Figure  20b.  The  net  result  of  the 
imbalance  of  heating  and  cooling  was  a  poleward  energy  transport,  R^, 
required  to  balance  the  budget  and  is  shown  in  Figure  20c.  Vonder  Haar 
and  Suomi  divided  the  total  energy  transport  among  three  major  compo- 
nents; the  free  energy  transport  of  the  ocean  with  a  mean  of  50  Cal 
m~^d~  ,  the  transport  of  energy  in  the  latent  heat  of  water  vapor  with 
a  mean  of  25  Cal  m~2d_1,  and  the  sensible  heat  plus  potential  energy 
transported  by  the  atmosphere  with  a  mean  of  75  Cal  m"2d"1.  Thus  ap- 
proximately 100  Cal  m"2d"1  is  atmospherically  transported  towards  the 
poles. 

The  buildup  of  potential  gradients  of  energy  was  released  in  the 
model  through  the  mean  circulation  flux  of  meridional  motion  and  eddy 
flux.  The  100  Cal  m~2d_1  transported  to  higher  latitudes  was  used  as 
the  energy  source  for  the  consumptive  processes.  Byers  (1974)  sug- 
gested that  the  split  was  approximately  equal  between  the  circulation 
flux  of  the  atmosphere  and  the  eddy  flux  between  the  equator  and  30° 
latitude,  the  net  direct  effect  was  from  the  circulation  flux  while 
at  higher  latitudes,  the  eddy  flux  term  was  dominant.  Thus  50  Cal 
m~2d_1  were  pumped  from  the  order  to  disorder  storages  in  the  model 
to  represent  circulation  flux  while  the  other  50  Cal  m  2d  l   was  used 
to  power  the  eddy  flux  of  storms  in  the  disorder  pathway. 
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Oceanic  Circulation 

Figure  21a  presents  the  oceanic  systems  of  the  world  in  conceptual 
form  with  the  primary  energy  flows  evaluated.  Sluggish  general  circula- 
tion over  much  of  the  earth  converges  into  fast  currents  and  associated 
vortices  in  a  few  zones,  especially  in  the  Gulf  Stream  and  its  counter- 
parts in  other  oceans.  Figure  21b  is  the  oceanic  model  in  mathematical 
form  with  differential  equations.  Table  3  documents  the  evaluation  of 
energy  flows  and  storages  in  the  model  supported  by  calculations  and 
references  in  Appendix  I. 

The  oceanic  system  of  the  world  was  modeled  as  a  reservoir  of 
water  and  chemicals  in  continuous  circulation.  Solar  energy  in  the 
form  of  radiation,  wind,  and  rain,  plus  the  tidal  pull  of  the  sun  and 
moon  are  the  principal  energy  sources  responsible  for  the  cycling  of 
matter  in  the  sea.  Excess  heat,  built  up  at  lower  latitudes  in  the 
ocean  by  solar  heating,  is  transported  towards  the  poles'  tide  driven 
by  deep  currents  and  wind-driven  waves  and  surface  currents. 

Solar  energy  provided  the  major  energy  input  to  power  the  oceanic 
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system  by  delivering  3,280  Cal  m  d   (Weyl ,  1970).  The  sea  reflected 
back  an  average  of  350  Cal  m~  d~  to  the  sky  (Weyl,  1970),  consumed 
1,570  Cal  m~  d"  to  power  the  evaporation  process,  fixed  2  Cal  m  ^d 
in  photosynthesis  (Whittaker  and  Likens,  1975)  and  transported  50  Cal 
m'2d_1  toward  the  pole  (Vonder  Haar  and  Suomi ,  1971).  This  left  1,318 
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Cal  m  d  ,  which  was  lost  to  the  atmosphere  as  heat.  On  the  average, 
3.4  x  105  Cal  m~2  was  stored  in  the  sea  in  the  form  cf  heat  (Weyl,  1970) 

—  2  —  1 

The  evaporation  process  consumes  2.9  kg  m  d   of  water,  which  was  re- 
placed by  2.6  kg  in~2d_1  of  precipitation  and  0.3  kg  m"2d~1  of  runoff 
from  the  land  masses  of  the  earth. 


Figure  21.  Energy  circuit  models  of  the  oceanic  systems  of  the  earth. 

(a)  Conceptual  diagram  of  energy  concentration  of  sea 
that  illustrates  energy  flows  necessary  for  quality 
calculations.  Energy  flows  are  in  Cal  m-2d_1  and 
storages  in  Cal-2. 

(b)  Energy  circuit  model  of  oceanic  system  of  global 
model  in  Figure  18  with  differential  equations. 
Solar,  wind,  and  tidal  energy  produce  global 
currents  to  transport  excess  heat  at  lower  latitudes 
thus  aiding  in  equalizing  global  energy  distribution 
and  powering  the  cycles  of  the  sea. 
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Wind  and  tide  provided  the  necessary  energy  sources  to  transport 
the  excess  heating  at  the  equator  toward  the  poles  by  delivering  15 
Cal  m  d   of  wind  energy  of  which  0.44  Cal  m~2  was  stored  in  the  poten- 
tial kinetic  energy  transfer  of  waves,  while  0.13  Cal  m"2d"1  was  dissi- 
pated on  the  shores  of  the  continental  land  masses.  The  tide  driven  by 
the  gravitational  pull  of  the  sun  and  moon  contributed  1.5  Cal  m"2d_1  to 
the  ocean  system  by  powering  the  deep  currents  of  the  sea.  The  net  ef- 
fect of  the  wind  and  tide  was  a  transfer  of  50  Cal  m"2d_1  poleward  and 
a  substantial  mixing  of  the  sea. 

Just  as  the  global  imbalance  of  temperature  generates  large  scale 
patterns  of  circulation  in  the  atmospheric  system,  excess  heating  of 
the  sea  produces  gigantic  anticyclonic  gyers.  The  strongest  ocean 
currents  are  found  at  the  western  peripheries  of  gyres,  such  as  the 
Gulf  Stream  and  the  Kuroshio.  The  gyres  of  the  hemispheres  of  the 
earth  converge  by  colliding  equatorial  currents.  The  mean  radius  of 
gyres  is  approximately  2,500  km  and  the  average  velocity  of  the  current 
at  the  periphery  is  about  0.1  m  sec"1,  thus  the  period  of  circulation 
is  approximately  5  years  (Monin  et  al_. ,  1977).  The  mean  current  velo- 
city for  oceans  is  0.035  m  sec'1  thus  0.55  Cal  m'2  is  stored  in  the 
kinetic  energy  of  sea  circulation  (Monin  et  al_. ,  1977). 

The  glaciers  and  polar  ice  cap  store  2.9  x  107  km  (Skinner  and 
Turekian,  1973),  which  is  equivalent  to  5.7  x  10?  g  m"2  of  ice.  The 
turnover  time  for  the  ice  of  the  earth  was  estimated  at  4  x  10k  yr'1 
(Skinner  and  Turekian,  1973).  Therefore,  the  rate  of  new  ice  formation 
is  4  g  m  -<j  ,  which  is  offset  by  melting  of  an  equal  amount  if  the 
global  ice  storage  is  in  equilibrium.  This  melting  of  the  polar  ice 
cap  and  glaciers  requires  0.4  Cal  m"2d_I  of  heat  from  the  sea. 
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Photosynthetic  Production,  Animals,  and  Fire 

Figure  22  is  a  conceptual  model  of  the  order-disorder  model  of  the 
biosphere  with  animal  life  and  fire  in  the  consumer  role.  Figure  22b 
is  the  model  in  mathematical  form  with  the  storages  and  flows  of  energy 
evaluated  and  documented  in  Table  4.  Note  the  strong  similarity  between 
the  application  of  the  theory  of  order,  disorder,  and  pulsing  to  biolog- 
ical, atmospheric,  and  oceanic  systems  of  the  world  energy  web  in 
Figure  18. 

Many  fundamental  similarities  exist  between  the  biological  process 
of  material  cycles  that  form  ordered  structure,  the  structures  continual 
disordering  process,  and  the  atmospheric  and  oceanic  systems  of  the 
earth  previously  diagramed. 

The  photosynthetic  system  of  the  biosphere  forms  the  solar  energy 
concentration  system,  while  animals,  microorganisms,  and  forest  fire 
were  in  continual  competition  to  consume  the  excess  storages  of  energy 
in  the  plant  community,  to  aid  in  recycling  matter,  and  to  perform  many 
necessary  control  actions  by  feeding  back  part  of  their  concentrated 
energy. 

Sun,  heat,  wind,  rain,  and  geological  minerals  were  the  primary 
energy  sources  for  the  production  of  organic  matter  by  plant  life.  The 
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global  average  of  4,900  Cal  m  d   of  solar  energy  was  the  dominant 
energy  source  when  measured  in  heat  Calories.  Wind  and  rain  delivered 
approximately  15  Cal  m  d  .  Organic  and  geologic  minerals  contributed 
a  flow  of  0.25  g  m"2d_1,  or  1  Cal  m"2d_1. 

Lieth  and  Whittaker  (1975)  have  made  extensive  quantitative  studies 
of  the  primary  production  of  the  photosynthetic  processes  of  the  bio- 
sphere and  compiled  extensive  data  on  global  net  primary  production. 


Figure  22.  Energy  circuit  models  of  biological  energy  concentration  pro- 
cess of  the  world. 

(a)  Conceptual  version  of  order-disorder  process  in  biosphere 
with  energy  flows  in  Cal  nr2d-1  and  energy  storages  in 
Cal  m-2. 

(b)  Biological  section  of  global  energy  web  model  in  mathe- 
matical form.  This  cycling  receptor  module  was  evaluated 
for  energy  flows  in  biological  systems.  The  differential 
equation  form  of  the  model  is  also  included. 
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Based  on  annual  statistics  of  the  U.S.  Fire  Service,  it  was  estimated 
that  0.3  Cal  m  d   of  biological  material  was  consumed  by  fire. 

A  third  use  of  organic  matter  is  through  the  sedimentary  cycle. 
Schlesinger  (1977)  estimated  the  input  of  carbon  into  the  world 
detritus  pool  at  3.75  x  1013  kg  C  yr'1 ,   which  is  an  energy  input  of 
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2.1  Cal  m  d  .  The  world  stock  of  carbon  in  detritus  was  estimated 
1.4  x  1015  kg  C  which  is  equivalent  to  2.4  x  10k   Cal  m"2.  This  gives 
a  turnover  time  for  the  carbon  stock  of  about  40  years.  Hubbert  (1971) 
estimated  a  world  fossil  fuel  net  production  rate  of  0.1  Cal  m"2d"1. 
It  is  likely  that  some  detritus  is  degraded  in  the  fossil  fuel  produc- 
tion process.  A  minimal  consumption  of  detritus  by  the  earth's  fossil 
fuel  production  system  is  about  0.3  Cal  m~2d_1  or  approximately  10  per- 
cent of  the  total  detritus  produced.  They  found  the  net  primary  produc- 
tion of  the  biosphere  to  be  0.92  g  m~2d-1  dry  matter,  or  3.9  Cal  m"2d_1, 
of  chemical  energy  photosynthetically  fixed  by  the  biosphere.  This 
estimate  was  found  to  be  in  reasonable  agreement  with  the  findings  of 
other  researchers  (Golley,  1972;  Clson,  1970;  Ryther,  1969;  SCEP,  1970; 
Whittaker,  1970;  and  Whittaker  and  Likens,  1973).  Leith  (1975)  sug- 
gested that  these  estimates  may  not  need  major  revision  unless  human 
interference  changes  the  global  primary  production  drastically.  The 
total  estimated  biomass  for  plant  life  was  also  found  by  Leith  and 
Whittaker  (1975)  to  be  3,600  g  m"2  of  dry  matter  or  15,300  Cal  m"2  of 
chemical  energy. 

Consumer  systems  of  the  world  model,  which  were  powered  by  organic 
matter,  were  also  responsible  for  recycling  minerals  and  controlling 
the  function  of  the  photosynthetic  producer  systems.  Whittaker  and 
Likens  (1973)  found  consumption  of  organic  matter  by  herbivorous  animals 
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to  be  0.62  Cal  m  d  ,  and  the  consumption  of  man  to  be  about  0.27  Cal 

-2  -1  -2   1 

m  d  ,  for  a  total  consumption  of  0.9  Cal  m  d  .  They  also  estimated 
the  world  animal  biomass  to  be  4.5  g  m~  (dry  matter)  or  19  Cal  m~2 
which  was  only  .13  percent  of  the  standing  stock  of  plants. 

Forest  fire  was  of  considerable  importance  to  the  control  of  both 
plant  and  animal  populations.  Of  all  plant  biomass,  99.8  percent  was 
found  in  the  terrestrial  systems  of  the  world,  leaving  virtually  all 
of  the  chemical  potential  energy  stored  in  the  standing  stock  of  plants 
exposed  to  fire  (Whittaker  and  Likens,  1973). 
Plate  Tectonics,  Volcanoes,  and  Earthquakes 

Figure  23a  presents  the  global  geological  submodel  in  conceptual 
form  while  Figure  23B  is  the  mathematical  form  with  differential  equa- 
tions. Table  5  documents  the  model  and  summarizes  the  calculations  of 
this  section. 

The  geologic  processes  of  the  earth  are  of  considerable  importance 
in  the  cycles  of  the  biosphere  and  in  the  causes  of  disasters.  As  the 
tectonic  plates  of  the  earth  collide,  stress  energy  is  stored  in  the 
vicinity  of  the  plate  junctions.  Earthquakes  were  modeled  as  a  catas- 
trophic release  of  energy,  which  then  feeds  back  to  pump  more  energy  to 
drive  the  global  geologic  process.  As  the  plates  of  the  earth's  crust 
collide  the  thinner  oceanic  plates  may  be  driven  under  the  thicker  con- 
tinental plates,  folding  under  the  eroded  organic  matter  from  the  land 
and  the  sediments  deposited  on  the  continental  shelf. 

The  extreme  pressures  and  heat  transform  the  organic  matter  to 
fossil  fuel,  which  was  modeled  as  an  energy  source  for  the  earth's 
volcanic  system,  the  continental  uplift  process,  and  the  cities  of  the 


Figure  23.  Energy  circuit  model  of  earth's  geological  system. 

(a)  Conceptual  energy  circuit  model  of  geological  system, 
which  forms  the  basis  for  energy  quality  of  geologi- 
cal systems  of  the  earth.  Energy  flows  are  in  Cal 
nr2d_1  and  storages  in  Cal  m~2. 

(b)  Energy  circuit  model  of  the  plate  tectonic  system  of 
concentrates  energy  of  lissosphere  to  power  volcanoes 
and  earthquakes.  The  interaction  of  rain  to  produce 
floods  is  also  shown. 
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world.  As  the  ocean  plates  are  consumed  under  the  continental  plates, 
new  plates  are  formed  in  the  deep  trenches  of  the  ocean  from  uplift 
of  molten  mantle  material.  In  some  cases  the  formation  of  ocean  crust 
produces  sea  earthquakes  and  volcanoes.  These  catastrophic  events 
dissipated  some  of  their  energy  through  tsunami  or  seismic  sea  waves, 
which  accelerated  continental  erosion. 

Conceptually,  the  material  cycling  process  of  the  geologic  system 
of  the  earth  was  viewed  as  two  interlocked  order-disorder  loops.  In 
the  first  loop,  the  ocean  crust  was  continuously  recycled  into  the 
ocean  crust  mantel  of  the  earth  by  folding  the  ocean  crust  under  the 
continental  crust  at  the  plate  junctions  by  production  of  new  crust  at 
the  ocean  trenches  by  the  mantle.  Scientific  evidence  exists  to  prove 
that  the  ocean  plates  are  in  constant  motion.  The  ocean  plates  slide 
away  from  ocean  ridges  at  a  rate  of  approximately  2  cm  yr  ' ,  so  that 
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approximately  0.35  g  m  d   of  new  ocean  plate  material  was  formed  from 
hot  material  flowing  up  from  the  mantle  of  the  ocean  ridge  (Clark,  1971) 
As  the  ocean  plates  were  created  by  the  ridge,  the  plates  become  larger, 
applying  great  pressure  to  the  continental  plates,  which  were  about 
five  times  thicker  and  several  hundred  times  older  than  the  older  ocean 

-2 

plate.  Sediments  from  the  continental  plates  also  contribute  0.2  g  m 
d"1  to  the  mass  of  the  ocean  plates.  The  rate  of  creation  of  ocean 
crust  was  balanced  by  a  loss  of  0.58  g  m~  d"1  as  the  ocean  plates 
folded  under  the  continental  plates. 

The  fact  that  the  continental  crust  was  several  hundred  times 
older  than  any  ocean  crust  also  supports  the  concept  that  plates  of 
the  earth  were  part  of  a  separate  but  related  order-disorder  cycle 
in  which  the  erosion  process  of  the  atmosphere  wore  down  continental 
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plates  as  they  were  rebuilt  by  geostatic  forces.  The  continental  crust 
was  worn  down  over  time  by  climatic  forces  such  as  rain  and  wind  and 
transported  to  the  sea  by  the  world's  streamflow.  An  average  of  2.6  kg 
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m  d   of  precipitation  falls  on  the  ocean  annually,  while  1.8  kg  m  d 

-2  - 1 

falls  on  land.     Using  the  1,450  Cal   m     d     ,  which  was  consumed  by  evapora- 

- 1  _  -2-1 

tion  systems,  yields  0.6  Cal  g   of  rain.  That  is,  1,000  Cal  m  d   was 
used  to  drive  the  erosion  of  land  processes.  The  annual  runoff  to  the 

-2-1  -2-1 

sea  amounts  to  19.3  g  m  d  ,  which  carries  0.15  g  m  d  of  material 
from  the  continental  plate  to  form  sediment  on  the  ocean  crust.  Note 
that  the  ocean  crust  was  2.4  times  larger  than  the  continental  crust; 

-2-1  -2-1 

therefore,  0.51  g  m  d   of  continental  erosion  equals  0.21  g  m  d 

-2  -1 

of  ocean  sedimentation  and  the  0.51  g  m  d   loss  of  the  continental 
crust  was  offset  by  geostatic  and  volcanic  uplift  of  an  equal  amount. 
The  energy  required  to  lift  this  material  to  an  average  height  of  20 
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km  from  an  average  height  of  3.5  km  was  computed  to  be  0.05  Cal  m  d 

-2  - 1 

Hayakawa  (1976)  found  geologic  uplift  to  consume  0.025  Cal  m  "d  , 
which  was  in  reasonable  agreement  with  the  agreement  with  the  calcula- 
tion in  equation  13. 

Because  of  the  relatively  rigid  nature  of  the  plates  of  the  earth 
by  comparison  to  the  fluid  ocean  or  atmosphere, it  was  believed  that  the 
geological  process  of  floods,  glaciers,  volcanoes,  and  earthquakes  per- 
form an  important  role  in  plate  recycling.  Floods  were  modeled  as  an 
interaction  of  the  atmosphere  process  of  rain,  ocean  waves,  and  the 
continental  plates.  The  role  of  floods  in  the  global  web  was  modeled 
as  a  prime  transporter  of  material  from  the  continental  crust  to  the 
ocean  crust,  thus,  accelerating  the  cycling  of  the  geologic  system  of 
the  earth. 
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Flood  waters  have  very  high  soil-carrying  capacity,  so  that  turbu- 
lence built  up  in  floods  moves  yery   large  quantities  of  land  to  lower 
elevations.  Assuming  that  1  percent  of  all  rain  was  involved  in  flooding 
(yields  1.8  g  nf'd"1  global  average  of  flood  water  elevated  an  average  of 
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1,000  m)  results  in  a  gravitational  potential  energy  of  0.04  Cal  m  d 
for  floods. 

Volcanoes  provide  an  offsetting  force  in  the  contribution  to  conti- 
nental uplift  process.  The  average  global  energy  flow  was  estimated  by 
Hayakawa  (1976)  to  be  5  x  10-5  Cal  m"  d"  .  The  volcanoes  in  the  model 
were  powered  by  organic  sediments,  solar  heating,  and  the  interaction  of 
deep  heat  flowing  out  of  the  continental  surface.  An  equally  important 
geological  process  was  the  seismic  energy  release  of  earthquakes. 
Richter  (1958)  estimated  the  seismic  energy  global  average  to  be  1  x  10 
Cal  m  d  .  Earthquakes  were  modeled  as  the  catastrophic  release  of 
strain  energy  built  up  by  the  interaction  of  the  ocean  and  the  continen- 
tal plates  driven  by  the  creation  of  new  ocean  plates  and  the  conven- 
tion of  the  mantle.  Earthquakes  were  fed  back  to  accelerate  the  erosion 
process  by  fracturing  continental  plate  material. 
Urban  Structure,  Fire,  Disasters,  and  War 

Figure  24  is  an  energy  circuit  submodel  of  the  urban  sector  of  the 
world  system  showing  the  internal  disordering  effect  of  depreciation, 
obsolescence,  fire,  and  war  and  the  external  disordering  effect  of  floods, 
storms,  earthquakes,  volcanoes,  and  forest  fires.  Table  6  documents  the 
energy  flows  and  storages  in  the  model. 

The  cities  of  the  earth  formed  the  fifth  order-disorder  system  in- 
cluded in  the  global  model.  The  structure  or  capital  assets  of  the  ur- 
ban systems,  including  all  structures,  machinery,  and  human  assets,  were 


Figure  24.  Urban  order-disorder  model  with  depreciation,  fire,  war,  and 
natural  disasters  in  disorder  process.  Table  5  documents 
the  evaluation  of  the  flows  and  storages  in  the  model. 

(a)  Conceptual  model  with  energy  flows  Cal  nr2d-1  and 
storages  in  Cal  m-2. 

(b)  Mathematical  form  of  model  with  differential  equations. 

Q^  =  K6lQ2lQl9(Q7+Q9)QmQl6-K65Q22Ql9-K67Ql9Q23-K6^Ql9 

-K65Q19(Q23+Q18+Q10+Q6+Q20) 

Q22=^5-K53Ql'+ 

Q2  3=K6  7Q2  3"^6  8Q2  3"^6  5Ql9Q2  3- 
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modeled  as  order.  The  disordering  forces  modeled  were  depreciation,  ob- 
solescence and  urban  renewal,  urban  fire,  disasters,  and  war.  The  com- 
bined action  of  these  destructive  forces  was  used  to  cycle  ordered 
matter  in  the  form  of  city  structure  to  disorder. 

Cities  were  modeled  as  a  consumer  of  fossil  fuel  and  biological 
material.  The  urban  fossil  fuel  consumption  of  the  world  in  1937  was 
only  0.056  Cal  m"2d-1  (Oxford  Economic  Atlas  of  the  World,  1965).  By 
1974,  the  world  fossil  fuel  consumption  had  grown  to  0.23  Cal  m~2d_1 
(United  Nations  Statistical  Year  Book,  1975,  1976). 

The  world  consumption  of  biological  material  by  humans  was  estimated 
by  Whittaker  and  Likens  (1975)  to  be  0.27  Cal  m"2d_1,  yielding  a  total 
chemical  energy  consumption  for  man  of  0.45  Cal  m~2d_1.  The  energy 
storage  in  urban  structure  of  the  world  was  calculated  from  depreciation 
rates  of  assets.  Solow  (1962)  showed  that  incorporating  technical  ob- 
solescence and  physical  depreciation  gave  an  average  annual  loss  of  ap- 
proximately 8  percent  of  capital  stocks.  Based  on  an  energy  consumption 
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of  0.45  Cal  m  d  ,  a  depreciation  rate  of  the  capital  stock  of  the 
world  cities  was  calculated  to  be  2,000  Cal  m2.  By  comparison  this  was 
13  percent  of  the  biological  standing  stock  of  the  globe.  The  biolog- 
ical material  and  fossil  fuel  consumed  by  the  urban  sector  were  primarily 
carbon  and  thus  vented  in  the  form  of  gas  back  to  the  atmosphere.  Urban 
renewal  was  modeled  as  an  internal  recycle  process,  which  used  energy 
from  the  production  function  of  the  urban  sector  to  recycle  structure. 

Urban  fire  was  the  largest  destructive  force  with  an  estimated  mag- 
nitude of  four  times  the  total  of  all  other  catastrophic  losses.  Eased 
on  fire  statistics  for  the  United  States,  it  was  estimated  that  urban 
fire  consumes  at  least  4  percent  of  the  total  productivity  and  1  percent 
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of  the  urban  assets.  This  amounts  to  a  loss  of  nearly  0.06  Cal  m~^d_1 
for  urban  systems  due  to  fire. 

The  second  major  internal  distructive  force  was  revolution  and 
war.  War  was  modeled  as  the  consumption  of  excess  storages  in  urban 
systems,  which  was  then  fed  back  to  produce  disordering  of  the  urban 
assets.  The  energy  consumption  of  war  and  revolution  was  calculated 
by  taking  the  ratio  of  the  world  military  expenditures  to  the  world 
gross  national  product.  This  ratio  was  multiplied  by  the  world  energy 

-2 

consumption  to  get  war  consumption.  It  was  found  that  0.05  Cal  m 
was  consumed  by  war. 

The  destruction  of  urban  assets  by  natural  disasters,  such  as 
storms,  floods,  and  earthquakes,  was  modeled  as  a  stress  on  the  cities 
of  the  world.  The  total  loss  due  to  these  natural  disasters  was  found 
to  be  0.006  Cal  m~2d_1  for  urban  systems.  Global  losses  from  each 
disaster  were  estimated  and  used  as  a  basis  for  comparison  of  disaster 
severity.  Order-disorder  ratios  for  each  disaster  were  calculated  and 
are  presented  in  Table  7.  Global  average  energy  flows  in  each  disaster 
are  listed  along  with  order-disorder  ratios  in  Calories,  heat,  equiva- 
lent, and  Calories  solar  equivalent. 
The  Global  Energy  Web  and  Energy  Quality  of  Disasters 

The  energy  circuit  made  of  the  world  in  Figure  25  was  evaluated 
to  serve  as  a  guide  for  energy  quality  calculation  in  this  research. 

-2  -i 

The  model  was  powered  by  7,000  Cal  m  d   of  solar  energy  entering  at 
the  left  and  feeding  through  the  web  to  the  right.  In  each  process, 
considerable  energy  was  dissipated  in  the  form  of  waste  heat,  which 

-2-1 

when  collected  produced  7,000  Cal  m  d   of  infrared  radiation  outbound 
from  the  earth. 
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Energy  quality  factors  for  the  major  energy  flow  processes  of  the 
earth  are  summarized  in  Table  8.  The  energy  flow  in  the  process  is 

-2-1 

presented  in  Cal  m  d   and  the  mean  energy  factor  for  each  process  in 
Calories  of  solar  energy  per  Calorie  of  the  process.  Note  that  Table 
8  is  arranged  in  descending  order  of  energy  flow  and  thus  ascending 
order  with  respect  to  energy  quality.  The  energy  quality  of  electro- 
magnetic radiation  was  calculated  by  a  different  method.  For  infrared 
radiation  the  ratio  of  the  average  energy  per  photon  was  used.  This 
was  also  equal  to  the  ratio  of  the  average  temperature  of  the  energy 
source.  In  all  energy  quality  calculations  solar  energy  was  used  as 
the  reference. 

In  the  sections  that  follow  the  quality  factors  developed  based 
on  the  energy  web  of  the  biosphere  were  applied  to  specific  case  study 
order-disorder  events. 


Evaluation  and  Simulation  of  the  Radiation  Stress  on  a  Tropical 
Rain  Forest  at  El  Verde,  Puerto  Rico 


This  section  contains  the  results  of  evaluation  and  simulation  of 
an  order-disorder  model  of  the  effect  of  gamma  ray  irradiation  on  a 
tropical  rain  forest.  Figure  26  is  a  much  aggragated  conceptual  model 
of  the  tropical  rain  forest  at  El  Verde,  Puerto  Rico.  In  the  model, 
solar  energy  provides  the  necessary  power  to  develop  order  in  biomass 
and  diversity  from  soil.  Recycling  is  provided  by  respiration  of  plants, 
animals,  and  microbes.  The  structure  built  by  the  forest  ecosystem  was 
fed  back  to  provide  canopy  mechanisms  for  the  collection  of  solar  energy. 
The  model  was  evaluated  at  20  m  from  the  irradiation  source.  At  this 
distance  from  the  center  of  radiation  it  was  observed  that  the  structural 
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TABLE  8.  ENERGY  QUALITY  FACTORS  FOR  GLOBAL  PROCESSES  AND  DISASTERS 


Global  Process  or  Disaster      Global  Mean   Global  Energy    Note* 

Energy  Flow   Quality  Factor 
Cal  irr2  d'1     Cal  Cal"1 


Solar  energy 

to  earth 

7,000 

1.0 

2-1 

Solar  energy 

absor 

bed 

by  earth 

4,900 

1.4 

2-1 

Surface  radiation  i 

on  earth 

3,290 

2.1 

2-4 

Evaporation 

1,600 

4.4 

3-6 

Weather 

100 

70 

2-15 

Wind  and  rain 

i 

15 

4  x  102 

2-16 

Primary  production 

of 

plants 

7.8 

9  x  102 

4-6 

Hurricanes 

4 

1.7  x 

103 

2-18 

Tide 

1.54 

4.5  x 

103 

3-16 

Animal  respiration 

0.65 

1.1  X 

lO1* 

4-9 

Cities 

0.45 

1.5  x 

10u 

6-2 

Ice 

0.39 

2.0  x 

10^ 

3-23 

Forest  Fire 

0.3 

2.3  x 

10^ 

4-10 

Waves 

0.13 

5.4  x 

10^ 

3-13 

Fossil  fuel 

0.11 

6.4  x 

10^ 

5-19 

Urban  fire 

0.065 

1.1  x 

105 

6-8 

War 

0.05 

1.4  x 

105 

6-9 

Uplift 

0.05 

1.4  x 

105 

5-14 

Floods 

0.04 

1.7  x 

105 

5-17 

Earthquake 

0.001 

4.7  x 

106 

5-16 

Volcanoes 

0.0005 

1.4  x 

107 

5-22 

The  calculations  that  document  this  table  are  shown  in  Appendix  I. 


Figure  26.  Order-disorder  model  of  gamma  ray  stress  on  tropical 
rain  forest  at  El  Verde,  Puerto  Rico.  Solar  energy, 
nutrients  and  biomass  interact  to  produce  plant 
sugar,  which  is  transformed  into  biomass  and  struc- 
ture. Biomass  is  recycled  through  litter  to  micro- 
bial and  animal  consumers.  With  the  supplementary 
stress  generated  by  the  gamma  radiation,  additional 
biomass  is  recycled,  the  rate  of  which  is  a  func- 
tion of  the  radiation  intensity. 
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loss  due  to  radiation  equaled  the  structure  generated  by  primary  produc- 
tion (Odum,  1970c). 

Figure  27a  is  the  energy  circuit  model  used  for  the  simulation  of 
the  order  and  disorder  cycling  process  in  the  rain  forest  ecosystem. 
The  storages  and  pathway  are  labeled  with  their  mathematical  coeffi- 
cients in  differential  equation  form.  Figure  27b  is  the  model  with 
the  storages  and  flows  evaluated. 
Evaluation  of  the  Rain  Forest  Model 

Table  9  summarizesthe  energy  calculations  necessary  to  evaluate 
the  model  in  Figure  24,  while  detailed  calculations  are  presented  in 
Appendix  I.  All  calculations  were  made  in  Calories  to  simplify  sim- 
ulation of  the  model.  The  solar  insulation  J  was  measured  at  the 
site  to  be  3,830  Cal  m"2d_1  (Odum  et  al_. ,  1970). 

The  rain  forest  in  El  Verde,  Puerto  Rico,  had  a  gross  primary 

-2  - 1 

production  of  32  g  m  d   organic  matter,  which  v/as  equivalent  to 
145  Cal  m  2d  1   (Odum,  1970a).  The  production  of  organic  matter  by 
the  photosynthetic  process  of  the  rain  forest  ecosystem  was  balanced 
by  respiration  of  the  plants,  animals,  and  microbes.  The  loss  due  to 
respiration  was  equal  to  gross  production  at  20  meters  (Odum,  1970c). 
The  total  standing  stock  of  the  biomass  Q1   was  found  to  be  42  kg  m"2 
(Odum,  1970c). 

The  10,000  curie  cesium  137  source  R  was  found  to  deliver  744  Cal 
d   in  gamma  rays.  To  calculate  the  radiation  dose  rate  received  by 
the  forest  ecosystem,  the  source  was  assumed  to  be  a  point  source.  By 
considering  inverse  square  loss  and  radiation  absorption  by  forest 
vegetation  and  atmosphere,  the  exposure  rate  at  20  m  was  calculated  to 


Figure  27.  Evaluated  order-disorder  model  of  gamma  ray  radiation 
stress  on  tropical  rain  forest  in  El  Verde,  Puerto 
Rico.  Table  9  documents  the  numerical  evaluation  of 
the  model . 

(a)  Energy  circuit  model  with  mathematical  coeffi- 
cients and  differential  equation  used  in 
simulation. 

(b)  Evaluated  energy  circuit  model. 
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Figure  28.  Gamma  radiation  field  and  effects  on  tropical  rain 
forest  at  El  Verde. 

(a)  Summary  of  radiation  field  at  close  range  as 
determined  by  dosimeters  (Odum  and  Drewry, 
1970). 

(b)  Graphs  of  optical -density  characteristics 
for  the  South  Control  Center  and  the  Radia- 
tion Center  (Desmarais  and  Helmuth,  1970). 
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Figure  29.  Results  of  simulation  of  order-disorder  gamma  ray 
model  in  Figure  27. 

(a)  Plot  of  order  (biomass)  Qi  and  disorder 
(litter)  Q2  in  kilograms  per  meter  squared 
versus  time  after  radiation  source  is  con- 
tinuously activated  in  days  for  distance 

of  1,  2.5,  5,  10,  20,  30,  and  50  m  from  the 
source. 

(b)  Plot  of  order  (biomass)  in  kilograms  per 
meter  squared  versus  distance  from  the 
source  after  93  days. 
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be  1.28  x  10   Cal  cm  d  ,  which  was  also  verified  by  comparison  to 
field  measurement.  See  Figure  28  from  Odum  and  Drewry  (1970).  The 
absorbed  dose  at  20  m  was  calculated  to  be  2.5  x  107  Cal  gm'V1  based 
on  a  measured  forest  biomass  of  42,450  g  m   (Odum,  1970c). 
Simulation  of  the  Order-Disorder  Process  in  the  Rain  Forest  Ecosytem 

Figure  29  demonstrates  the  result  of  the  tropical  rain  forest  order- 
disorder  model  simulation  covering  the  93-day  radiation  experiment  at 
radius  of  1,  5,  10,  20,  30,  and  50  m  from  the  source.  At  the  end  of 
the  period  of  irradiation,  the  radiation  source  was  turned  off  to  let 
the  system  rebuild.  Note  that  in  the  simulation  results  at  20  m,  the 
growth  of  the  forest  was  just  offset  by  the  radiation  stress. 

Figure  29b  presents  the  simulated  radiation  dose  and  the  dis- 
ordering effect  of  the  irradiation  with  distance.  Biomass,  as  measured 
by  optical  density  of  the  rain  forest  standing  stock  Qi,is  plotted  in 
Figure  28b  and  compared  with  measurements  of  radiation  dose  (Figure  28a). 
Note  the  similarity  between  the  measured  and  simulated  biomass  standing 
stock. 
Calculation  of  Order-Disorder  Energy  Ratios  for  the  Rain  Forest 

A  previous  calculation  of  the  order-disorder  ratio  for  the  tropical 
rain  forest  system  at  El  Verde  was  performed  by  Odum  (1970c)  comparing 
the  dose  rate  to  photosynthetic  energy.  Odum  found  1  Cal  of  gamma  ener- 
gy was  required  to  disorder  what  10,000  Cal  of  photosynthetic  energy 
could  order.  The   order-disorder  ratio  for  the  tropical  rain  forest  was 
also  calculated  from  the  evaluation  of  the  model  in  Figure  27.  The  gross 
primary  production  of  137  Cal  m~2d~  was  divided  by  the  absorbed  dose  of 
1.1  x  10"2  Cal  :i"2d~1.  This  yields  an  order-disorder  ratio  of  13,000 
Cal  of  production  per  Cal  of  radiation. 
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The  disaster-amplifier  ratio  was  calculated  by  dividing  the 
total  loss  of  biomass  due  to  the  irradiation  by  the  total  absorbed  dose 
over  the  93-day  exposure.  At  20  m  from  the  source  10  g  m~^d  *was  dis- 
ordered by  the  radiation.  This  amounts  to  a  total  loss  of  3.9  x  103 
Cal  m"2  over  the  period  of  exposure.  The  total  absorbed  dose  was  1 
Cal  m"2.  This  yields  a  disaster-amplifier  ratio  of  3.9  x  103  Cal  of 
biomass  per  Cal  of  radiation. 

For  compairson  with  other  systems,  ratios  were  also  calculated  in 
solar  equivalents.  The  embodied  solar  energy  in  biomass  was  calculated 
from  the  quality  factor  for  gross  primary  production  in  Table  8.  Two 
methods  were  used  to  estimate  a  range  of  quality  factors  for  gamma 
radiation. 

The  first  method  is  based  on  the  ratio  of  blackbody  equivalent 
temperature  of  the  cesium  gamma  source  to  the  blackbody  temperature 
of  the  sun.  From  Weins  law  the  blackbody  temperature  of  6.6  x  10 5  ev 
gamma  rays  is  1.5  x  109°K.  The  blackbody  temperature  of  the  sun  is 
6,000°K  (Gates,  1962).  This  gives  a  low  estimate  of  the  energy  quality 
of  the  gamma  radiation  of  2.5  x  105  Cal  solar  per  Cal  radiation. 

The  second  estimation  is  based  on  the  ratio  of  the  energy  required 
to  produce  the  10,000  Ci  cesium  source  to  the  744  Cal  d" 1  of  radiation 
produced  by  the  source.  The  cesium  source  cost  $25,000  in  1965  (Cdum, 
1970a).  Using  the  ratio  of  2.5  x  10u  Cal  coal  equivalent  per  dollar 
(Odum  and  Odum,  1976)  and  a  quality  factor  for  Cal  of  6.4  x  1014  Cal 
solar  per  Cal  coal  from  Table  3,  yields  a  total  of  4  x  1013  Cal  solar 
energy  required  to  produce  the  10, COO  Ci  radiation  source.  The  half 
life  of  cesium  137  is  30  yr   (Gloyna  and  Ledbetter,  1969).  Starting  from 
744  Cal  d"1  a  total  of  3.1  x  106  Cal  of  radiation  emitted  over  the  life  of 
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the  source.  This  results  in  a  high  estimate  of  the  energy  quality  of 
the  gamma  radiation  of  5  x  105  Cal  solar  per  Cal  gamma  radiation. 

With  a  low  estimate  of  2.5  x  105  Cal  per  Cal  and  a  high  estimate 
of  5  x  106  Cal  per  Cal  an  intermediate  quality  factor  of  1  x  106  Cal 
solar  per  Cal  gamma  radiation  was  used.  Applying  the  intermediate 
quality  factor  to  the  absorbed  dose  at  20  m  gives  an  absorbed  dose 
ra 
93  d  exposure.  Solar  energy  embodied  in  the  gross  primary  production 


te  of  1  x  10"2  Cal  solar  m~2d_1  and  a  total  dose  of  9.3  x  105  for  the 


_2  - 1 

was  calculated  by  multiplying  the  137  Cal  m  d   production  by  the 
primary  production  quality  factor  in  Table  8,  yielding  1.23  x  105  solar 
Cal  m~2d_1.  This  gives  an  order-disorder  ratio  of  12  solar  Cal  produc- 
tion per  solar  Cal  of  radiation.  The  disaster-amplifier  ratio  was  also 
calculated  in  solar  equivalents.  The  total  loss  of  biomass  was  3.5  x 
105  solar  Cal  m~2  due  to  a  total  absorbed  dose  of  9.3  x  105  solar  Cal. 
The  disaster-amplifier  ratio  was  thus  4  solar  Cal  biomass  per  solar 
equivalent  Cal  of  radiation. 

Evaluation  and  Simulation  of  the  Johnstown  Flood  of  1389 

Figure  30  is  an  energy  circuit  model  of  the  failure  of  the  South 
Fork  Dam  and  the  destruction  of  Johnstown.  Pennsylvania,  in  1889.  The 
flow  of  water  into  the  reservoir  was  modeled  as  the  interaction  of  the 
land  in  the  South  Fork  basin  and  the  rain.  The  model  was  designed  so 
that  when  the  water  in  the  spillway  exceeded  capacity  the  overflow  was 
fed  back  to  rupture  the  dam.  The  flood  wave  grew  exponentially  until 
the  source  of  energy,  the  reservoir,  was  exhausted,  generating  an 
energy  pulse.  The  pulse  was  transmitted  down  the  Conemaugh  Valley  in 


Figure  30.  Energy  circuit  model  of  Johnstown,  Pennsylvania,  1889, 
and  the  disruption  generated  by  the  catastrophic 
failure  of  the  South  Fork  Dam.  Table  11  documents  the 
evaluation  of  the  Johnstown. 

Qi  -  Ki  J  Q2  -  K3  Qi  -  K„  Qx  -  K3  Ql   Q3 

Q3  =  K3  Ql   +  K5  Qi  Q3  -  K6  Q3  -  K7  Q3  Q„ 

Q2  =  QT  -  Ky  Qi 

Qu  =  (Km  -  K13)  Qk   Q5  Q6  Q7  -  K10  Q6 

Qs  =  C  -  Kn  Q4 

Q6  =  K3  R  -  K9  Q„  Q5  Q6  Q7  -  K10  Q5 

Qy  =  -K12  Q^  Qs  Qs  Q? 
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the  form  of  a  flood  wave.  The  destruction  of  Johnstown  was  performed 
by  using  the  pulse  energy  to  disorder  the  town. 
Evaluation  of  the  Model 

Table  10  documents  the  evaluation  of  the  model  Johnstown  flood 
of  1889.  To  analyze  and  model  the  disruption  and  disorder  created  by 
the  great  flood,  it  was  necessary  to  estimate  the  energy  released. 
The  flood  was  modeled  as  the  interaction  of  the  geologic  processes  of 
the  earth  and  the  hydrological  system  of  the  atmosphere  that  produces 
storms. 

The  South  Fork  reservoir  Q  drained  1.5  x  108m2.  Not  all  of  the 
4  x  10~5m3sec_1  rainfall  J  was  observed  flowing  into  the  reservoir  in 
the  form  of  runoff.  Based  on  an  observed  rate  of  rise  in  the  reservoir 
of  7  x  10~5m  sec"1  (Engineering  Society  of  Western  Pennsylvania,  1889) 
and  a  reservoir  surface  area  1.69  x  105m2,  a  net  increase  of  118  m3 
sec   of  water  was  stored  in  the  South  Fork  reservoir.  In  addition, 
the  spillway  of  the  South  Fork  Dam  K6Q3  was  designed  to  release  99  m3 
sec"1  (Engineering  Society  of  Western  Pennsylvania,  1898).  However, 
Shappee  (1940)  estimated  the  peak  capacity  to  be  169  m3sec-1.  Com- 
bining the  net  increase  in  storage  with  the  spillway  loss  yields  a 

3    ~  1 

total  inflow  of  368  m  sec  ,  or  approximately  58  percent  of  the  rain- 
fall. 

The  potential  energy  stored  in  a  cubic  meter  of  water  in  the  South 
Fork  reservoir  calculated  by  equation  35. 

E  =  2.3  x  10"1*  mgh  Cal  m"3,  (35) 

where  m  =  mass  of  one  m3 

h  =  elevation  of  South  Fork  reservoir  with  respect  to  Johnstown 
g  =  acceleration  of  gravity 
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the  gravitation  potential  energy  was  found  to  be  6.94  x  103  Cal  m"3. 

The  total  energy  consumption  for  Johnstown  was  estimated  to  be 
1.4  x  109  Cal  d  ,  of  which  71  percent  was  coal  and  the  remaining  29 
percent  was  mainly  lumber,  firewood,  and  food.  The  net  production  of 
the  city  was  1.2  x  109  Cal  d"1,  which  was  used  to  offset  depreciation 
and  consumption  of  1.8  x  108  Cal  d"1.  The  remaining  4  x  107  Cal  d"1 
was  consumed  by  the  growth  of  the  city.  This  would  be  an  annual  growth 
rate  of  6  percent  which  was  supported  by  the  rapid  industrial  expansion 
of  Johnstown  between  1880  and  1889  reported  in  Engineering  News  (1889a). 
The  tangible  assets  for  Johnstown  Qh   were  estimates  based  on  the  His- 
torical Statistics  of  the  United  States  (1975)  and  found  to  have  an 
energy  storage  of  1 .8  x  1012  Cal.  This  was  approximately  four  times 
the  annual  energy  consumption. 

The  failure  of  the  South  Fork  Dam  produced  a  wave  that  took  56 
min  to  travel  the  5.57  x  lOSn  to  Johnstown  at  an  average  speed  of  7.66 
m  sec  .  As  the  wall  of  water  moved  down  the  Conemaugh  Valley,  the 
water  from  the  torrential  rain  lubricated  the  wave;  however,  the  top 
of  the  water  moved  faster  because  of  lower  friction,  creating  a  ball 
of  water  and  debris  which  rolled  down  the  valley.  When  the  wave  hit 
Johnstown,  the  height  was  observed  to  be  10  m.  However,  it  was  as 
high  as  30  m  in  a  narrow  gorge.  The  kinetic  energy  of  the  10  m  wave 
which  hit  Johnstown  traveling  at  a  velocity  of  7.66  m  sec"1  was  made 
up  of  two  components.  First,  the  linear  motion  term  described  by 
equation  36. 

eL  =  |  mV2  =  123  Cal  m"3,  (36) 
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which  yields  123  Cal  m   of  kinetic  energy.  However,  the  turbulence 
generated  by  the  rolling  action  of  the  wave  also  contained  considerable 
kinetic  energy  as  given  in  equation  37  (Kinsman,  1965). 

Et  =  0.3  H2  Cal  m"2  (37) 

where  H  is  the  height  of  the  wave  in  meters  based  on  the  observed  10  m 
wave  height  the  turbulence  was  32.3  Cal  m~2. 

To  calculate  the  kinetic  energy  of  the  flood  wave  the  linear 

_  3 

energy  of  123  Cal  m   was  multiplied  by  the  volume  of  water  in  the 

flood,  1.36  x  107m3  (Engineering  News,  1889b),  yielding  a  total  linear 

kinetic  energy  release  of  1.67  x  109  Cal.  The  total  turbulent  kinetic 

energy  release  was  calculated  by  multiplying  the  area  of  flood  wave 

1.36  x  106m2  by  the  turbulent  energy  32.3  Cal  m"2,  which  gave  4.4  x 

107  Cal.  Thus  a  total  kinetic  energy  of  1.7  x  109  Cal  was  released  in 

the  flood.  The  kinetic  energy  K7Q3QI+  of  the  flood  was,  therefore,  1.8 

percent  of  the  potential  energy  stored  in  the  reservoir.  Application 

of  the  energy  quality  factor  for  floods  of  1.7  x  105  solar  Cal  per  Cal 

yields  3  x  101!+  Cal  of  solar  energy  embodied  in  the  flood. 

Simulation  of  the  South  Fork  Dam  Failure  and  the  Destruction  of 
Johnstown 

Figure  31a  is  the  results  of  the  simulation  of  the  dam  failure  with 
time.  The  total  storage  in  the  South  Fork  reservoir,  the  kinetic  ener- 
gy in  the  form  of  turbulence,  and  the  rate  of  overflow  over  the  dam  or 
through  the  dam  are  plotted  as  a  function  of  time.  Note  that  once  a 
critical  overflow  was  reached,  the  growth  of  the  kinetic  energy  in  the 
dam  failure  wave  was  only  limited  by  the  quantity  of  water  remaining 
in  the  reservoir  to  power  the  wave. 


Figure  31.  Simulation  results  of  Johnstown  flood  of  1889  in 
Figure  30. 

(a)  Results  of  dam  failure. 

(b)  Results  of  distruction  of  Johnstown. 
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Figure  31b  presents  the  results  of  simulation  of  the  catastrophic 
destruction  of  Johnstown,  Pennsylvania.  The  results  of  the  simulation 
approximate  qualitative  descriptions  of  the  event.  Coefficients  were 
calculated  from  the  flow  and  storage  evaluation  in  Table  10. 
Disaster-Amplifier  Ratios  for  Johnstown  and  the  Flood 

The  ratio  of  the  energy  loss  in  the  disaster  to  the  energy  in  the 
flood  was  calculated  in  Calories  heat  equivalent  and  calories  solar 
equivalent.  The  energy  delivered  to  Johnstown  by  the  flood  was  found 
to  be  1.7  x  109  Cal .  The  total  catastrophic  depreciation  undammage 
Ki+6Qt+Q3  was  T-S5  x  1012  Cal.  Dividing  the  total  destruction  by  the 
energy  required  to  produce  the  disorder  gave  a  disaster-amplifier  ratio 
at  970  Cal  of  destruction  to  1  Cal  of  flood  energy. 

For  comparison  purposes  the  disaster-amplifier  factor  was  recalcu- 
lated with  both  energy  flows  in  solar  equivalents.  Energy  quality 
factors  for  floods  and  cities  were  taken  from  Table  3.  The  embodied 
energy  in  the  flood  was  2.9  x  1011+  solar  Cal  while  the  energy  embodied 
in  the  loss  of  city  structure  was  2.4  x  1016  solar  Cal.  The  disaster- 
amplifier  ratio  was  thus  85  solar  Cal  of  city  per  solar  Cal  of  flood 
energy. 

Evaluation  and  Simulation  of  the  Guatemala  Earthquake  of  1976 

Figure  32  is  an  energy  diagram  of  the  Guatemala  earthquake  and  the 
destruction  of  El  Progreso.  The  interaction  of  the  North  American  Plate 
with  the  Caribbean  Plate  converges  stress  energy  into  the  Motagua  Fault 
system.  When  a  threshhold  is  reached,  excess  energy  is  released  in  the 
form  of  an  earthquake  pulse.  The  seismic  pulse  is  transmitted  through 
the  crust  of  the  earth  and  coupled  to  secondary  systems.  The  destruction 


Figure  32.     Earthquake  order-disorder  model   for  Guatemala  earthquake. 
See  Table  11   for  energy  flow  and  storage  explanations. 

*i  =  MTTT  "  M1Q3-  K2Q1 

Q2  =  Ce  -  K15Q2 

Q3  ■  ^QiQs  +  K7Q1-  K5Q3-  K0Q3Q4 

Q4  =  K9E  Q^Qs  -  K10Q3Qt+  -  KuQi, 

Qs  =  Ct  -  Kl3Qk 
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of  El  Progreso  is  performed  by  disordering  the  structure  of  the  town. 

The  earthquake  energy  pulse  provides  the  energy  for  the  recycling 

process. 

Evaluation  of  the  Model 

Results  of  evaluating  the  model  are  given  in  Table  11.  The  stress 
was  a  buildup  of  the  Motagua  Fault  system  as  the  North  American  Plate  Ej 
moved  west  at  a  rate  from  2.1  to  2.2  cm yr-1  relative  to  the  stationary 
Caribbean  Plate  (Jordan,  1975;  Plafker,  1976).  The  total  change  in 
strain  energy  during  the  earthquake  was  estimated  by  multiplying  the 
strain  release  per  unit  volume  by  the  volume  of  fault  relaxed.  The 
United  States  Geological  Survey  found  this  to  be  a  low  stress  drop 
event  with  between  3  and  18  bars  of  strain  relieved  by  the  rupture 
around  the  fault.  The  fault  break  was  240  km  long,  7  km  deep,  and  a 
slip  of  3  m  wide  (Espinosa  ejt  aj_. ,  1976).  This  yields  an  estimated 
active  fault  volume  Q2  of  4.83  x  1016  cm3.  Releasing  3  bars  of  strain 
over  the  whole  volume  produced  a  spontaneous  release  of  3.1  x  1012  Cal 
during  the  earthquake.  The  result  of  this  calculation  was  in  good 
agreement  with  the  2.62  x  1012  Cal  reported  (Espinosa  et  al . ,  1976). 

The  change  in  strain  energy,  Qls  was  accounted  for  either  by 
released  seismic  energy,  E<-,  and  by  work  done,  W,  by  sliding  on  the 
fault  as  shown  in  equation  37  (Dieterich,  1974): 

Qi  =  Es  +  W  (37) 

Many  authors  report  methods  for  relating  the  seismic  energy  release, 
in  ergs,  to  the  magnitude  of  the  quake  on  the  Richter  scale  M.  The 
recommended  form  is  given  in  equation  38  (Gutenberg  and  Richter,  1969): 

log  E$  =  9.4  +  2.14M  -  9.954M2  (38) 
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This  gives  a  seismic  energy  release  of  2.28  x  1022  ergs  (5.44  x  TO11 
Cal )  for  the  M  =  7.5  earthquake  under  study  (Espinosa  et  aj_. ,  1976). 
From  equation  37,  it  was  found  that  75  percent  of  the  energy  released 
in  the  Guatemala  earthquake  was  dissipated  in  moving  the  fault,  while 
the  other  25  percent  was  released  as  seismic  energy. 

The  2.2  cm  yr~     westward  velocity  of  the  North  American  Plate 
relative  to  the  Caribbean  Plate  for  203  years  yields  a  total  offset 
of  446.6  cm.  However,  a  maximum  offset  of  340  cm  was  observed  with 
an  average  displacement  of  108  cm  (Plafker,  1976).  Comparison  of  the 
average  stress  buildup  motion  to  the  average  displacement  gives  an  ef- 
ficiency of  24.1  percent.  The  loss  of  the  stress  not  accounted  for  in 
the  earthquake  was  accounted  for  by  many  minor  slips  over  the  past 
several  centuries.  Harlow  (1976)  was  able  to  show  that  over  100  of  the 
mini-earthquakes  occurring  the  year  before  the  major  event  resulted 
from  minor  slips  in  the  Motagua  Fault  (K2Qi). 

Assuming  that  the  mean  of  the  mini-earthquake  released  5.37  x  108 
Cal  of  strain  energy  explained  the  difference  between  the  calculated 
446  cm  and  the  observed  108  cm.  The  average  mini -quake  magnitude  would 
be  a  little  greater  than  four  on  the  Richter  scale.  This  magnitude  is 
in  agreement  with  reported  shocks  (Harlow,  1976).  To  build  up  the 
necessary  storage  of  2.62  x  1012  Cal  to  supply  the  earthquake  in  207 
years,  1.26  x  1010  Cal  of  energy  is  required  to  offset  the  depreciation 
due  to  losses  in  the  fault  resulting  from  minor  slips.  Therefore,  a 
total  of  6.6  :<  1C10  Cal  yr~l   was  delivered  by  the  2.1  cm  yr~     plate 
motion. 

Table  12  summarizes  the  numeric  calculation  necessary  to  relate 
the  known  earthquake  intensity  to  the  distribution  of  siesmic  energy 
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dissipated  in  Cal  m   for  each  isoseismal  zcne  in  Guatemala.  Figure  33 
is  the  results  of  a  Log-log  plot  of  the  distribution  of  seismic  energy 
throughout  Guatemala.  Double  log  coordinates  were  used  so  that  the 
large  range  of  values  could  be  included. 
Simulation  of  the  Earthquake  and  Secondary  Distribution 

Figures  34a  and  34b  are  the  results  of  simulation  of  the  catastro- 
phic release  of  energy  in  the  Guatemala  earthquake  and  the  secondary 
destruction  of  El  Progreso  as  modeled  in  Figure  32.  In  Figure  35a  the 
energy  released  is  plotted  as  a  function  of  time.  Once  the  earthquake 
is  triggered,  the  event  grows  exponentially  until  the  strain  in  the 
fault  is  relieved.  This  generates  a  pulse  of  siesmic  energy,  which  is 
then  transmitted  radially  in  a  control  action.  Figure  34b  is  a  plot 
of  the  destruction  of  El  Progreso  by  the  disordering  action  of  the 
siesmic  pulse. 
Calculation  of  Disaster-Amplifier  Ratios 

With  the  distribution  of  earthquake  energy  known  by  zone  for 
Guatemala  and  the  percentage  of  destruction  calculated  for  each  zone, 
the  relationship  between  earthquake  energy  dissipated  and  local  dis- 
order and  distribution  was  examined. 

The  gravitational  potential  energy  for  a  typical  Guatemalan 
house  was  calculated  from  equation  39  (Spears  and  Zemansky,  1955)  and 
found  to  be  97  Cal . 

Vg  =  mgh,  (39) 

where  V  is  the  gravitational  potential  energy,  m  is  the  mass  of  the 
house,  g  is  the  gravitational  constant,  and  h  is  the  mean  height  cf 
the  mass.  The  foundation  area  of  the  house  was  6  m2;  therefore,  the 
gravitational  potential  energy  is  15  Calm"2. 


Figure  33.  Energy  distribution  in  Guatemala  from  the  major  event, 
February  4,  1976  earthquake.  Log-log  plot  of  land 
area  affected  versus  earthquake  energy-dissipated. 
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Figure  34.  Simulation  results  of  earthquake  model  in  Figure  32. 
Seismic  energy  release  and  destruction  of  El  Progreso, 
Guatemala,  versus  time  after  earthquake. 
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Figure  35.  Structural  damage  of  towns  in  Guatemala  as  a  function 
of  earthquake  energy  dissipated  by  the  adobe  houses. 
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When  the  energy  of  the  earthquake  was  dissipated  over  the  land  not 
all  of  the  energy  was  absorbed  in  the  man-made  structures.  The  wave 
moved  through  a  crust  7  km  thick,  so  one  might  picture  a  house  sitting 
on  a  column  of  earth  7  km  high.  A  1  m2  column  of  earth  would  have  a 
mass  of  1.54  x  107  kg  while  a  square  meter  of  house  has  a  mass  of  512 
kg.  Since  the  houses  are  made  of  earth,  linear  distribution  of  energy 
was  assumed  throughout  the  entire  column  of  earth  crust  and  the  house. 
Equation  40  gives  the  amount  of  energy  absorbed  in  the  house. 

eh  =  ei  ^mh/mc^  (40) 

where  e.  is  the  energy  absorbed  by  the  house,  e.  is  the  energy  dissi- 
pated per  square  meter  in  zone  i,  m,  is  the  mass  of  the  house,  and  m 
is  the  mass  of  the  square  meter  column  of  earth  crust.  Column  9  of 
Table  12  shows  energy  absorbed  per  square  meter  of  house  based  on  the 
above  method. 

From  the  energy  per  unit  area  absorbed  in  the  house  and  the  en- 
ergy per  unit  area  stored  in  the  structure,  the  effect  of  the  two  for 
towns  throughout  various  intensity  zones  in  Guatemala  was  compared. 
Table  13  summarizes  the  findings.  The  structural  energy  released  by 
the  earthquake  was  calculated  by  taking  the  gravitational  potential 
energy  stored  in  the  structures  of  the  town  and  converting  it  to 
kinetic  energy  based  on  the  percentage  damaged.  The  disaster-amplifier 
factor  was  then  calculated  for  each  city  by  dividing  the  gravitational 
kinetic  energy  released  by  the  earthquake  seismic  energy  absorbed  by 
the  structure.  Figure  35  shows  the  results  of  relating  the  simulated 
destruction  to  earthquake  seismic  energy  dissipated  by  adobe  structures 
in  Guatemala.  It  was  found  that  on  the  average  1  Cal  of  earthquake 
wave  energy  released  7,716  Cal  of  gravitational  potential  energy  by 
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toppling  adobe  houses.  The  magnitude  of  the  earthquake  energy  dissipated 
per  square  meter  was  calculated  from  reported  field  intensities  collected 
by  the  United  States  Geological  Survey  (Espinosa  et  aK  ,  1976)  and  mapped 
in  Figure  35.  The  damage  data  were  reported  by  the  Guatemalan  National 
Reconstruction  Committee  (Comite  de  Emergencia,  1976). 

The  disaster-amplifier  ratio  was  also  calculated  using  Calories  of 
embodied  solar  energy.  The  mean  seismic  energy  absorbed  per  square  meter 
of  structure  and  the  quantity  of  urban  structure  disorder  were  converted 
to  solar  Calories.  The  energy  quality  factors  presented  in  Table  8  were 
used  for  this  calculation.  It  was  found  that  25  solar  Cal  of  structure 
were  disordered  by  1  solar  Cal  of  seismic  energy. 


DISCUSSION 

From  models  and  measurements  of  disasters  a  general  pattern  of  pulsing 
was  found  for  many  aspects  of  nature.  Interdependence  of  pulses  was  re- 
lated through  a  world  web  model  as  a  general  theory  of  disasters. 

The  Dynamics  of  Pulsing  Systems 

The  exponential  surge  model  diagramed  in  Figure  2  was  found  to  apply 
to  all  cases  examined.  It  has  an  order-disorder  model  with  a  catastro- 
phic consumer  in  the  recycle  loop.  This  model  uses  pulsing  to  improve 
material  cycling  between  the  storages  of  order  and  disorder  in  Figure  1, 
and  accounts  for  the  catastrophic  properties  of  disasters.  The  trans- 
formation from  a  stable  system  to  a  pulsing  action  may  be  either  self- 
activated  or  externally  triggered. 
The  Cycling  of  Matter 

Conservation  of  matter  was  incorporated  in  the  models  by  requiring 
that  the  sum  of  the  ordered  and  disordered  matter  equal  a  constant. 
Structure  was  built  from  disorder  by  ordering  matter.  Likewise  disorder 
was  created  from  order  by  recycling  matter.  In  the  basic  exponential 
surge  model,  two  recycle  pathways  exist.  In  the  linear  pathway  the 
recycle  process  is  continuous  and  proportional  to  the  quantity  of  matter 
stored.  When  sufficient  stocks  are  accumulated  for  catastrophic  re- 
cycling to  be  supported,  a  catastrophic  consumer  is  generated  that  is 
capable  of  ^ery   rapid  recycle  and  energy  pulse  generation. 
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The  ability  to  cycle  matter  may  be  a  fundamental  design  principle 
of  self-sustaining  systems.  When  systems  build  up  large  storages  of 
matter  they  may  become  sluggish  and  suffer  from  stagnation.  Many 
forces,  both  internal  and  external,  continually  work  to  recycle  the 
stored  matter  and  prevent  stagnation.  The  exponential  surge  model  was 
used  to  diagram  the  systems  of  the  global  energy  web,  the  Johnstown 
flood,  and  the  Guatemala  earthquake.  The  simulation  results  presented 
illustrate  the  material  cycling  character  of  this  flexible  configura- 
tion of  the  order-disorder  model. 
Energy  Convergence  and  Storage 

The  exponential  surge  configuration  of  the  order-disorder  model 
presented  in  Figure  2  and  used  throughout  this  research,  uses  the 
ordered  structure  and  the  recycled  matter  in  the  form  of  disorder  to 
converge  energy  and  grow  order  from  disorder.  The  ability  to  concen- 
trate and  store  energy  by  transforming  disordered  matter  to  order  is 
believed  to  be  a  property  of  self-sustaining  systems.  This  model  is 
capable  of  self-design  in  the  sense  that  it  can  switch  energy  from 
developing  order  to  control  and  recycle  as  either  one  of  the  other 
processes  becomes  limiting. 
Exponential  Surges  and  Energy  Pulses 

When  the  quantity  of  energy  converged  in  the  form  of  ordered 
structure  is  great  enough  to  trigger  catastrophic  consumption,  the 
exponential  surge  model  catastrophically  recycles  the  matter  stored 
and  releases  its  energy  storage  in  a  control  action.  The  ability  to 
generate  control  pulses  is  believed  to  be  a  common  design  principle 
of  self-sustaining  systems.  The  magnitude  of  the  energy  pulse  re- 
leased is  proportional  to  the  energy  stored  in  the  form  of  order. 
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The  pulse  width  and  height  is  set  by  the  dynamic  properties  of  the 
consumer  module  which  catastrophically  recycles  the  ordered  matter. 
System  Control  through  Cascading 

The  catastrophic  energy  pulses  generated  by  one  system  often 
control  other  systems.  The  diagram  in  Figure  3  illustrates  this 
concept.  Several  possible  control  mechanisms  exist.  In  some  catas- 
trophic events  the  disaster  pulse  energy  provides  the  power  to  recycle 
the  ordered  structure  by  the  interaction  of  the  stored  energy  and  the 
disaster  pulses.  In  other  cases  the  disaster  pulse  triggers  a  sepa- 
rate secondary  disaster  such  as  a  dam  failure  or  a  fire.  A  fourth 
design  principle  suggested  for  self-sustaining  systems  is  the  ability 
to  send  and  receive  information  in  the  form  of  energy  pulses. 
Dynamic  Behavior  of  Energy  Cycling  Models 

The  application  of  the  order-disorder  model  in  Figure  18  to  the 
disordering  effect  of  radiation  stress  on  a  rain  forest  ecosystem  ap- 
proximated the  dynamic  nature  of  the  actual  system  and  was  stable. 
When  the  model  was  applied  quantitatively  to  the  Johnstown  flood  and 
the  Guatemala  earthquake,  it  was  found  to  have  stable  dynamic  charac- 
teristics and  to  provide  graphs  replicating  the  observed  events. 

Global  Energy  Web 

The  energy  circuit  model  of  the  world  conceptually  presented  in 
Figure  18  and  evaluated  by  sector  in  Figures  19,  21,  22,  23,  and  24 
connects  many  order-disorder  models  in  a  web.  The  global  web  allows 
energy  flows  including  disasters  to  be  put  in  relative  prospective 
of  importance  and  impact. 
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Energy  Quality,  A  Universal  Energy  Conversion  Factor 

The  energy  quality  factors  computed  were  based  on  the  weblike 
nature  of  energy  flows  of  the  earth's  systems.  The  quality  factors 
are  the  ratio  of  the  total  solar  energy  required  to  power  the  energy 
convergence  process  of  the  earth  to  the  quantity  of  concentrated  energy 
flowing  in  the  system  of  interest.  Through  the  evolutionary  processes 
of  the  earth,  the  systems  best  adapted  to  the  use  of  all  pulsing  energy 
sources  present  may  have  prevailed. 

The  necessity  of  each  process  to  evolve  until  all  available  sources 
of  energy  are  harnessed  to  maximize  power  may  be  the  reason  the  atmosphe- 
ric, oceanic,  biological,  geological,  and  urban  sectors  of  the  world's 
system  have  developed  into  a  web.  Once  the  web  is  formed  and  evolved 
to  maximize  power  for  the  total  system,  no  sector  can  compete  by  oper- 
ating separetely  because  limitations  in  cycles  would  develop.  Leibin's 
(1862)  law  of  the  minimum  describes  the  behavior  of  processes  without 
good  recycle. 

As  solar  energy  is  converged  through  the  web  of  the  earth,  the 
energy  flowing  in  a  given  process  becomes  less  even  though  large  quan- 
tities of  dilute  solar  energy  were  necessary  to  power  the  concentra- 
tion process.  Thus  the  equality  of  the  energy  in  areas  of  the  flow 
decreases.  Lotka's  (1922)  principle  of  maximum  power  suggests  that 
systems  which  further  converge  high  quality  energy  use  it  effectively 
to  survive.  Energy  wasting  systems  are  eventually  replaced  by  systems 
that  maximize  power.  For  systems  which  have  stood  the  evolutionary 
test  of  time,  it  is  likely  that  the  quantity  of  energy  converged  to 
produce  the  energy  flow  will  be  equal  to  the  effect  of  the  energy 
flow  as  a  feedback  amplifier.  Thus  energy  quality  is  a  universal  way 
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of  comparing  different  forms  of  energy  flowing  in  the  weblike  process 

of  the  earth  showing  which  are  valuable  and  have  large  requirements  and 

impact. 

Embodied  Solar  Energy,  A  Universal  Measure 

The  solar  energy  embodied  in  a  flow  is  calculated  by  multiplying  the 
energy  flow  by  the  appropriate  global  quality  factor.  This  calculation 
estimates  the  quantity  of  solar  energy  converged  by  the  biosphere  to 
produce  a  flow  of  interest.  This  transformation  of  the  energy  flows 
converts  them  to  common  units  of  measure,  embodied  solar  Calories. 
This  makes  possible  the  comparison  of  two  processes.  Also  the  relative 
value  of  two  interacting  flows  can  be  determined.  Table  8  ranks  energy 
flows  according  to  their  value. 
Disater-Amplifier  Ratios,  A  Measure  of  Control 

When  disaster-amplifier  ratios  were  calculated  for  the  evaluated 
models  of  the  tropical  rain  forest,  Johnstown  flood,  and  Guatemala 
earthquake  in  Figures  27,  30,  and  32,  the  control  or  disordering  effect 
was  greater  than  the  quantity  of  energy  delivered  by  the  catastrophic 
pulsing  action  of  the  controlling  system.  Calculation  of  the  disaster- 
amplifier  ratio  in  heat  equivalents  ranged  from  1,000  to  10,000.  Lower 
amplifier  factors  ranging  from  10  to  100  were  found  when  calculations 
were  made  with  the  disaster  energy  pulse  and  the  quantity  of  disorder 
generated  both  measured  in  embodied  solar  Calories. 

The  appropriate  disaster-amplifier  ratios  multiplied  by  the  quan- 
tity of  energy  embodied  in  a  local  hazard  may  be  useful  in  estimating 
total  risk  from  all  hazards  in  specific  locations.  Risk  maps  con- 
structed from  plots  of  the  aggregate  risk  potential  may  be  of  con- 
siderable utility  for  guiding  local  land  use  decisions. 
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The  Frequency  of  Pulsing  and  Control 

The  ability  of  systems  of  man  and  nature  to  adapt  to  the  pulsing 
environment  is  probably  a  function  of  the  frequency  of  the  pulsing 
action.  Examples  of  adaptation  are  common.  Structures  built  and  eco- 
systems living  in  areas  of  daily  tidal  exchange  with  the  sea  show  many 
characteristics  of  adaptation  to  the  tidal  pulsing;  however,  adaptation 
to  less  predictable  pulses  such  as  earthquakes  is  less  pronounced.  The 
magnitude  and  frequency  of  the  pulsing  action  of  the  environment  is 
believed  to  be  of  fundamental  consequence  in  selecting  which  systems 
will  survive  the  test  of  time. 
Control  of  Global  Cycles  by  Solar  Pulses 

Empirical  evidence  links  11 -year  solar  cycles  to  the  cycles  of  the 
global  energy  web.  The  physical  mechanism  for  the  linkage  is  not  well 
known  except  for  the  upper  atmospheric  phenomena.  Sunspct  numbers  are 
correlated  with  the  aurora  boreal  is.  Auroral  displays  are  produced 
when  charged  particles  from  the  sun  interact  with  the  earth's  magnetic 
field,  resulting  in  particle  acceleration  and  collisions  with  air  mole- 
cules in  the  upper  atmosphere  (Eddy,  1976).  Bray  (1967)  found  signifi- 
cant correlation  between  sunspot  and  biological  productivity  through 
radiocarbon  activity  from  527  B.C.  to  1964  A.D. 

The  high  quality  pulses  generated  by  solar  flares  may  couple  to 
the  order-disorder  cycles  of  the  earth,  synchronizing  the  two  processes. 
The  appropriate  disaster-amplifier  factor  for  sunspots  might  be  similar 
to  the  one  calculated  for  the  controlling  effect  of  gamma  radiation  on 
the  tropical  rain  forest  in  Figure  27.  An  energy  quality  factor  for 
sunspots  could  be  obtained  by  estimating  the  quantity  of  solar  energy 
converged  to  produce  a  solar  flare.  The  product  of  the  disaster- 
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amplifier  ratio  and  the  energy  embodied  in  sunspot  energy  incident 
would  then  produce  and  estimate  the  degree  of  control  over  world's 
cycles  by  solar  cycles. 

Figure  36  is  a  conceptual  diagram  suggesting  a  mechanism  for 
coupling  the  cycles  of  the  sun  to  the  cycles  of  the  earth.  Preliminary 
estimates  based  on  disaster-amplifier  ratios  of  gamma  radiation  and 
energy  quality  factors  based  on  the  blackbody  temperature  of  solar 
flares  point  to  sufficient  energy  in  sunspot  activity  to  control 
meteorological,  biological,  geological,  and  human  cycles  of  the 
earth. 

Disasters 

Natural  and  man-made  disasters  have  had  considerable  effect  on 
the  history  of  the  earth.  Many  relatively  smooth  periods  of  slow 
physical  and  cultural  changes  have  been  violently  disrupted  and 
forever  altered  by  the  catastrophic  release  of  energy  in  fires, 
wars,  earthquakes,  floods,  and  storms. 
The  Energy  Cost  and  Benefit  of  Disasters 

The  results  of  the  research  in  this  work  suggest  that  energy 
source-pulse  systems  provide  benefits  considerably  larger  than  harm 
to  the  world  environment.  The  ratio,  of  the  quantity  of  work  done 
to  the  energy  required  to  produce  the  disordering,  was  greater  than 
one  in  each  case.  Comparison  between  the  continuously  disordering 
effect  of  gamma  radiation  on  a  tropical  terrestrial  ecosystem,  the 
two  pulse  disordering  process  of  a  flood  and  an  earthquake  suggest 
that  the  disordering  pulse  is  more  effective  than  an  equivalent 
amount  of  energy  expended  on  3  continuous  basis. 


Figure  36.  Conceptual  model  illustrating  possible  mechanisms 
for  control  of  global  energy  cycles  by  Solar 
Pulsing. 
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When  calculated,  on  the  basis  of  Calories  of  heat  equivalents,  1.1 
x  10"2  Cal  m"2d_1  of  gamma  radiation  disordered  the  work  of  145  Cal  m"2 
d  1   of  photosynthetic  energy  or  one  Calorie   of  disorder  was  offset 
by  10,000  Calories  of  order.  The  same  calculation,  for  the  devastation 
of  Johnstown  by  the  failure  of  the  South  Fork  Dam,  yielded  1.7  x  105 
Cal  m"2  of  devastation  produced  by  1.3  x  103  Cal  m"2  of  flood  energy 
or  approximately  one  Cal  of  disorder  to  offset  1,000  Cal  of  order. 
In  the  third  case  in  El  Progreso,  Guatemala,  3.2  x  10"3  Cal  m"2  of 
energy  from  the  Guatemala  earthquake  of  1976  produced  14.5  Cal  m"2  of 
release  of  gravitational  potential  energy  of  structure  or  approximately 
one  Cal  of  earthquake  energy  absorbed  the  disordering  of  7,700  Cal  of 
gravitational  potential  energy.  In  each  case  the  effect  of  one  Calorie 
of  disordering  energy  was  far  greater  than  the  effect  of  one  Calorie 
of  ordering  energy.  However,  the  above  ratios  may  be  of  little  value 
since  the  calculation  of  the  ratios  were  based  on  the  heat  equivalent 
Calories  rather  than  a  common  denominator  such  as  solar  equivalents. 
The  Distribution  of  Disaster  Pulses  as  a  Control  Device 

The  control  action  of  energy  pulses  generated  by  the  catastrophic 
disordering  of  energy  concentrating  systems  of  the  earth  is  distributed 
conversely  to  the  concentration  gradient.  As  diagramed  in  the  global 
model  in  Figure  18,  energy  is  concentrated  by  all  of  the  ordering  pro- 
cesses of  the  earth,  and  when  catastrophically  released  is  fed  back  to 
control  the  system. 

Throughout  history,  the  occurrence  of  catastrophic  events,  both 
natural  and  man-made,  have  been  used  as  milestones  in  time  and  space. 
Such  terms  as  pre-World  War  I  are  common  place  in  literature  and  eyery 
day  conversation.   Such  references  suggest  that  the  control  action  of 
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disaster  pulses  on  the  events  of  the  world  is  generally  recognized. 
Additional  evidence  for  the  acceptance  of  the  large  effect  of  disas- 
ters on  the  behavior  of  systems  is  implied  by  the  deliberate  hesita- 
tion of  many  analysts  to  model  time  periods  which  contain  catastro- 
phic events.  Also  a  caveat  is  normally  included  with  economic  and 
population  forecasts  which  specifically  excludes  the  event  of  war  or 
natural  disaster  from  the  projections. 

The  global  model,  developed  and  evaluated  in  this  research  and 
presented  in  overview  form  in  Figure  10,  addresses  the  control  func- 
tion of  disasters  on  the  world  energy  web.  Many  important  material 
recycling  processes  are  heavily  dependent  on  energy  pulses.  Further- 
more the  recycling  of  materials  helps  to  maximize  the  power  of  dis- 
ordering systems  which  have  stagnated. 

Specific  examples  of  higher  productivity  in  natural  systems  after 
disordering  have  been  sited  by  various  researchers.  For  example,  Lugo 
and  Snedaker  (1974)  suggested  thatmangroves  organized  their  extent  of 
regrowth  according  to  hurricane  frequency  to  maximize  productivity.  The 
role  of  forest  fire  as  a  regenerative  force  in  nature  has  been  well  docu- 
mented by  Komarek  (1962)  at  the  Tall  Timbers  Research  Station.  The  long- 
term  effect  of  disasters  on  the  productivity  of  urban  systems  has  also 
been  shown  to  be  positive.  After  extensive  economic  analysis  of  the 
long-term  effect  of  disasters  on  communitites,  Dacy  and  Kunrenther 
(1969)  concluded  that  communities   may  benefit  economically  from  di- 
sasters through  rapid  inflow  of  capital  for  rebuilding  and  the  adop- 
tion of  technological  innovations  to  meet  critical  situations.  The 
field  research  in  El  Progreso,  Guatemala  and  Johnstown,  Pennsylvania, 
support  the  positive  benefit  concept.  In  El  Progreso,  in  addition  to 
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large  quantities  of  immediate  relief,  an  entirely  new  structural 
system,  including  water,  sewerage,  roads,  and  drainage,  was  installed 
and  all  of  the  government  buildings  were  replaced  with  more  modern  and 
functional  facilities.  In  the  case  of  Johnstown,  all  three  major 
floods  have  performed  urban  renewal  functions  and  stimulated  new  econo- 
mic investment.  The  rapid  economic  growth  of  Germany  and  Japan  after 
heavy  destruction  from  World  War  II  are  also  examples  of  the  positive 
benefit  of  recycling  obsolete  structure  by  externally  driven  catastro- 
phic forces. 
Common  Aspects  of  Disaster  and  Disorder  in  Different  Systems 

As  previously  stated  one  of  the  objectives  of  this  research  was 
to  develop  an  energy  basis  for  disasters  and  to  show  that  they  were 
an  essential  part  of  the  cycles  of  order  and  disorder  of  the  systems 
of  the  earth.  One  of  the  ways  this  idea  was  addressed  in  this  re- 
search was  to  demonstrate  that  a  similar  process  of  order,  disorder, 
and  pulsing  exist  throughout  the  major  systems  of  the  earth.  All  of 
the  disaster  generating  models  in  this  research  use  a  common  format 
for  the  energy  convergence  and  catastrophic  pulsing.  Note  the  simi- 
larity in  the  energy  convergence  process  for  storms,  floods,  fires, 
war,  earthquakes,  and  volcanoes  as  illustrated  by  Figures  18,  19,  21, 
22,  23,  and  24.  This  is  further  illustrated  by  the  simulations  of 
flood  and  earthquake  destruction  modeled  in  Figures  30  and  32,  re- 
spectively. 

Pulsing  and  Maximum  Power 

The  field  results  in  this  research  suggest  that  primary  and 
lateral  systems   can  be  more  productive  if  pulsed.  To  test  the 
feasibility  of  this  concept,  a  simple  experiment  was  carried  out. 
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Two  order-disorder  models  with  identical  energy  sources,  differential 
equations,  and  coefficients  were  constructed  as  shown  in  Figures  37a 
and  b.  The  model  in  Figure  36a  was  pulsed  disordered  every  25  time 
periods.  The  model  in  Figure  37b  was  allowed  to  grow  undisturbed.  See 
Figures  37c  and  d.  Also  plotted  in  Figure  37d  is  the  result  of  the 
integration  of  the  gross  production  of  both  systems  (pathways  marked  P). 
The  results  are  identical  during  the  initial  growth  period.  At  the  end 
of  25  time  periods,  system  a  was  disordered  and  then  allowed  to  rebuild. 
For  this  case,  the  total  production  was  stimulated  by  the  recycle  gene- 
rated by  the  pulse,  producing  a  greater  total  production. 

The  results  of  this  hypothetical  case  suggest  that  thepulsedis- 
ordering  of  systems  maximizes  power  by  recycling  the  matter  stored  in 
the  system  which  then  may  stimulate  regrowth.  The  counter  intuitive 
nature  of  this  results  raises  many  interesting  questions  about  the  role 
of  pulses  in  systems  of  man  and  nature.  Applications  in  the  world  pos- 
sibly range  from  the  pulse  harvest  of  agricultural  systems  to  possibly 
explaining  the  vitality  of  Germany  and  Japan  post-World  War  II  destruc- 
tion. 

Disaster  Planning  and  Future  Research 

White  and  Haas  (1975)  have  pointed  out  the  critical  need  for  anal- 
itical  tools  to  support  the  development  of  urban  policy  in  hazard  prone 
areas.  Several  possible  applications  of  the  theory  developed  in  this 
research  to  urban  problems  are:  (1)  demonstration  of  similarity  between 
problems  generated  by  natural  and  man  made  disaster  situations;  (2)  cal- 
culation of  the  total  risk  expressed  by  specific  community  to  all  hazards; 


Figure  37.  A  test  of  the  effect  of  pulsing  on  the  principle  of 
maximum  power. 

(a)  Pulsed  order-disorder  model.  The  gross  produc- 
tion of  the  system  is  summed  over  time  by  Q3. 

Qi  ■  K1E1Q1Q2  -  K2(h-SW 

Q2  =  C  -  Qi 
Q3  =  K1E1Q1Q2 

(b)  Order-disorder  model,  not  pulsed,  used  as  a 
control.  The  gross  production  of  the  system 
is  summed  over  time  by  Q7. 

Qs  =  KjEiQsQe  -  Ms 
Qe  =  C  -  Q5 

Q7  =  K1E1Q3Q6 

(c)  Results  of  simulation  of  pulsed  mode!  a. 

(d)  Results  of  simulation  of  model  b  and  plot  of 
integral  of  production  of  models  a  and  b. 
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(3)  simulation  of  effect  of  future  disasters  and  effectiveness  of  dis- 
aster mitigation  techniques;  and  (4)  preparation  of  local  maps  of  total 
risk. 

Effective  public  policy  can  only  be  formed  with  an  adequate  under- 
standing of  the  public.  Too  often  the  many  types  of  disaster  potentials 
that  threaten  the  structure  of  communities  are  treated  as  separate  and 
dissimilar  events  (White  and  Hass,  1975).  An  awareness  of  the  similarity 
of  disaster  mechanism  and  disruptive  effects  is  necessary  to  make  effi- 
cient use  of  public  resources.  In  addition  the  cascading  of  disasters 
renders  planning  for  a  specific  type  of  disruption  unrealistic.  A  well 
trained  and  integrated  civil  defense,  fire,  police,  and  emergency  govern- 
ment agencies  are  essential  to  properly  dealing  with  the  multiple  hazards 
of  urban  areas. 

Disaster-amplifier  ratios  and  the  energy  quality  factors  developed 
may  have  application  in  calculating  the  total  risk  a  particular  community 
faces.  A  quantitative  understanding  of  the  magnitude  of  the  aggregate 
risk  would  provide  a  base  for  public  policy  decisions  involving  commit- 
ment of  resources  to  disaster  mitigation  and  preparedness. 

Simulation  of  the  possible  effects  of  future  disasters  on  a  speci- 
fic community  in  a  manner  similar  to  the  models  developed  for  El  Progreso 
and  Johnstown  forms  a  quantitative  basis  for  disaster  planning.  Models 
incorporating  disaster  mitigation  techniques  could  be  used  to  quantita- 
tively evaluate  various  policy  alternatives  in  disaster  planning. 

Local  risk  maps  could  be  constructed  by  summing  the  result  of  con- 
verting spacial  distribution  of  various  forms  of  pulsing  energy  to  em- 
bodied solar  energy  by  using  quality  factors  from  Table  8  and  disaster 
amplifier  factors  from  appropriate  cases.  The  resulting  multi-hazard 
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risk  map  could  be  used  to  more  adequately  allocate  urban  land  between 
various  uses.  Possibly  land  uses  with  minimum  structure  would  be 
allocated  to  high  hazard  areas. 

It  is  recommended  that  a  specific  locality  be  selected  for  test 
application  and  further  development  of  the  energy  basis  for  disasters. 


APPENDIX  I 
NOTES  TO  TABLES  2-7,  9,  TO  AND  11 


Notes  to  Table  2 

(1)  Total  solar  flux  received  by  Earth  I  per  square  meter  of  surface 
area. 

I  =^x  E 
xo   As  x  L 

where  Ap  is  the  effective  area  of  earth's  plane 

As  is  the  surface  area  of  the  earth  (Van  Nostrands  Scientific 

Encyclopedia,  1976) 

E  is  the  solar  constant  19.4  Cal  m   min   (Byers,  1974) 

.    1.278  x  IP11*  m  x  19.4  Cal  m"2min"1x  1.44  x  103  min  d"1 
lo   "  5.101  x  IP14  m^ 

I  =  7, POP  Cal  m~2d_1 
o 

(2)  Global  average  solar  radiation  absorbed  by  earth  R  (Vonder  Haar 

and  Suomi ,  1971) . 

R  =  3.4  Cal  m"2  min"1  x  1 .44  x  IP3  min  d"1 
n 

Rn  =  4,9CP  Cal  m"2d_1 

(3)  Albedo  radiation  reflected  by  earth 
Albedo  A 

R       .  nnn  r..,     -2.-1 


a  =  l   _H      4,9CP  Cal  m  d  .,n 

A   '  "  I0   "  "  7,PPP  Cal  m"2d-i  "  ,JU 

Ra  =  7,PPP  Cal  m"2d_1  -  4,9PP  Cal  m~2d_1 
a 


Ra  =  2, IPO  Cal  m"2d_1 

a 

(4)  Surface  radiation  heat  loss  due  to  cooling  and  used  by  plants. 
47%  of  solar  flux  (Houghton,  1954] 

7,CPP  Cal  m"2d_1  x  .47  =  3,29C  Cal  m"2d_1 

182 


183 


Notes  to  Table  2  -  Continued 

(5)  Energy  used  in  powering  atmosphere 

23%  of  solar  flux  I 
o 

7,000  Cal  m"2d_1x  0.23  =  1,610  Cal  m'V1 

(6)  Global  average  infrared  radiation  R  outbound,  (Vonder  Haar  and 
Suomi,  1971). 

RL  =  3.4  Cal  m"2  min"1  x  1,400  m"2  d"1 

RL  =  4,900  Cal  m'V1 

(7)  Energy  stored  in  latent  heat  of  world  atmosphere 

Total  atmosphere  water  per  m2       (Skinner  and  Turekian,  1973) 

(1.3  x  10Hm3)  x  (1Q9ni3km3)  x  (1Q6  g  m"3)  x  0.574  Cal  q"1  _ 
5.10  x  10^m2 

1.46  x  102  Cal  m"2 

(8)  Mass  of  atmosphere  based  on  mean  pressure  of  1.013  x  105 

dynes  cm"2  (Goody  and  Walker,  1972) 

im-c-  _  (1.013  x  106  dynes  cm"  )(1  x  10^  cm?  m"2)  (g  cm  sec"2  dynes"1) 
(981  cm  sec-2)        (1  x  103  kg-*) 

Mass  =  1.04  x  10^  kg  m"2 

(9)  Energy  of  global  circulation  and  turbulence  (Monin  et  al . ,  1977). 

(1  x  1021  J)(2.30  x  10""  Cal  J"1)  _  -7n  .  -  -2 
(5.1  x  10^  W) 470  Cal  m 

(10)  Energy  absorbed  by  atmosphere 

23%  of  solar  (Houghton,  1954) 

7,000  Cal  m"2d_1x  0.23  =  1,610  Cal  m"2d-1 

(11)  Evaporation  from  sea  and  land       (Skinner  and  Turekian,  1971) 
3.83  x  1020  g  yr~l 

Evaporation  from  land:  0.63  x  1020  g  yr~l 
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Notes  to  Table  2  -  Continued 

Average  evaporation  per  square  meter  of  earth 

[(3.83  +  .63   )  x  IP20  g  yr  -1]  _  „  dnn  n     -2.-1 
(5.10  x  !W]xp  dyr-i   J  '  2'400  9  m     d 

Evaporation  energy 

2,400  g  m~2d-1  x  .6  Cal  g  =  1,400  Cal  m"2d_1 

(12)  Turbulence  (Byers,  1974) 
50%  of  the  kinetic  energy  of  atmosphere  in  eddy flux 

0.5  x  300  Cal  m'V1  =  150  Cal  m'V1 

(13)  Circulation  (Byers,  1974) 
50%  of  the  kinetic  energy  of  atmosphere  in  circulation 

0.5  x  300  Cal  irfV1  =  150  Cal  nfV1 

(14)  Cooling  due  to  radiation  loss  by  atmosphere 

Heat  loss  from  atmosphere  due  to  cooling  by  first  law  of  thermo- 
dynamics 
2,720  Cal  nfV1 

(15)  Poleward  flow  of  energy  in  weather  system,  see  Figure  20  (Vonder 
Haar  and  Suomi ,  1971) . 

100  Cal  nfV1 

(16)  Total  energy  in  winds,  water,  convection      (Hubbert,  1971) 
15  Cal  m'V1 

(17)  Heat  loss  in  circulation  and  eddy  flux  of  atmosphere  by  first  law 

150  Cal  m~2d_1  +  150  Cal  m"2d_1  -  (100  Cal  m"2d_1  +  15  Cal  nf2d_1) 
=  185  Cal  m"2d_1. 

(18)  Hurricane   and  cyclone   energy  (Dunn  and  Miller,  1964): 

(9.5  x  1015  Cal  d-1)(9  d  hurricane-1)^  hurricane  yr~l)    .   ,   -  ,  -2.-1 


(5.1  x  10il+  W]    (365  d  yr-1) 
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Motes  to  Table  3 


(1)  Solar  energy  47%  of  total  radiation  (Weyl ,  1970) 
7,000  Cal  m'V1  x  .47  =  3,290  Cal  m~2d_1 

(2)  Solar  reflection  5%  of  total  radiation  (Weyl,  1970) 
7,000  Cal  m~2d_1  x  .05  =  350  Cal  m'V1 

(3)  Energy  stored  in  heat  of  sea  (Weyl,  1970) 
3.4  Cal  cm"2  x  1  x  104  cm2  m"2  =  3.4  x  104  Cal  m"2 


(4)  Mass  of  sea  (Skinner  and  Turekian,  1973) 

1.4  x  IP21  kg  _  o  oo  y  ins  kp  m-2 
3.61  x  lQ>m2  "  3-88  X  10  kS  m 

(5)  Energy  stored  in  sea  gyres 

Mean  gyre  velocity  of  3.5  cm  sec-  (Monin  et_  al_. ,  1977) 

e  =  1/2  mv2  =  (3.88  x  106  kg  m2)(.035  m  sec-1)(2.38  x  lO'^Cal  J  T1 

2 

e  =  0.55  Cal  m-2 

(6)  Evaporation-energy  loss  in  latent  heat  (Skinner  and  Turekian,  1973) 

3.83  x  IP20  g  yr"1 _  9  «  v  1n3  „  -2.-1 

(3.61  x  lO1^  m2)  x  (365  d  yr  !)  a 

0.54  Cal  g"1  energy  for  evaporation 

2.9  x  103  g  m"2d_1  x  0.54  Cal  g"1  =  1,570  Cal  m"2d_1 

(7)  Heating,  solar  energy  consumed  in  heating  of  sea  by  conservation 
of  energy,  notes  1,  2,  6. 

(3,290  -  350  -  1,570)  Cal  m"2d-1  -  1,370  Cal  m"2d_1 

(8)  Transport  of  heat  in  sea  water  toward  the  poles  by  tidal  energy 
50  Cal  m"2d_1  (Vonder  Haar  and  Suomi ,  1971) 
Assume  50%  by  tide  and  50%  in  wind  drives  surface  current 

50  Cal  m'V1  x  0.5  =  25  Cal  m"2d-1 
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Notes  to  Table  3  -  Continued 
(  9)  Wind  pumped  surface  current 

-2  -1 

50  Cal  m  d    transport  of  heat  toward  the  poles  (Vonder  Haar 

and  Suomi,  1970). 

Remaining  50%  not  transported  by  tide 

50  Cal  m'V1  x  0.5  =  25  Cal  nfV1 

(10)  Heat  pumped  by  wind  driven  waves  47%  of  cooling 

Excess  heat  =  [3,290  +  1,570  +  50)]  Cal  nfV1  =  1,320  Cal  irfV1 
(1,320  Cal  nfV1  x  .47)  =  620  Cal  nfV1 

(11)  Heat  lost  due  to  surface  cooling  (not  evaporation  or  wind  generated) 
by  law  of  energy  conservation 

1 ,320  Cal  m~  d"  excess  heat 

-620  Cal  m  d   wind  driven  wave  cooling 
700  Cal  m_2d-1  surface  cooling 

(12)  Surface  current  drive  by  wind 
25  Cal  m~2d_1  (see  note  9) 

(13)  Energy  dissipated  by  waves  (Kinsman,  1965) 

(7  x  IP26  ergsyr"1)  x  (2.39  x  IP11  Cal  erg"1)   n  .,  .  .  -2.-1 
(3.61  x  10il+  m2)  x  (365 1  yr-1) =  °-13  Cal  m  d 

(14)  Wind  cooling 

620  Cal  m~2d_1  (see  note  10) 

(15)  Precipitation  (Skinner  and  Turekian,  1973) 

3.47  x  1020g  yr'1 _a  _, 

3.61  x  10^  m2)  (365  d  yr"l)  "  2'633  g  m  d 

(16)  Energy  input  of  tide  dissipated  by  sea  currents  (Weyl ,  1970) 

(2.7  x  1019  ergs  sec"I)(8.64  x  105  sec  d"M(2.39  x  lP_11Cal  erg"1) 

(3.6  x  10^m2) 

1.54  Cal  m~2d_1  (Skinner  and  Turekian,  1973) 
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Notes  to  Table  3  -  Continued 

(17)  Energy  input  of  wind  (Hubbert,  1971) 
15  Cal  m~  d" 

(18)  Energy  input  to  sea  gyres  (Monin  et  aj_. ,  1977) 
Velocity  of  gyer  at  periphery  10  cm  sec" 

Mean  gyer  radius,  2,500  km 

2  (2,500  km)  (1  x  105  cm  km"1)    _  ,  R   ,n3  ri-i 
(10  cm  sec" i)  (8.64  x  W*   sec  d'1)   ' *°  x  ,u  a 

r-  _  0.55  Cal  m   _  -,  „  in-v,i  m-2^-1 
E  =  1.8  x  10yd   ~  3  x  10  Cal  m  d 

(19)  Energy  dissipated  in  gyers  by  first  law  of  thermodynamics 
Dissipation  set  equal  to  consumption 

(20)  Storage  of  energy  in  ice  (Skinner  and  Turekian,  1973) 

Volume  of  water  in  polar  ice  cap  and  glaciers 

2.9  x  107  km3  x  1  x  109  m3  km"3  x  1  x  106  g  m"3  _  ,-  fiq   in?    -2 
5.1  x  10it  m2  °'Dy  x  1U  9  m 

(21)  Mean  residence  time  for  water  in  ice  caps  and  glaciers 
4  x  1Qk   years 

_,   _  storage 5.69  x  107  g  m~2 _ 

ow  "  mean  residence  time   (4  x  10*  yr  i)(365  d  yr-i)  " 

3.9  g  m"2d-1 

(22)  Rate  of  melting  set  equal  to  rate  of  buildup 
3.9  g  m"2d_1 

(23)  Energy  required  to  melt  ice  at  3.9  g  m"2d_1  rate 
Energy  stored  in  ice  relative  to  mean  global  temperature 
Energy  in  AT  300°  -  273°K  =  23°K 

Specific  heat  of  1  Cal  g_1°  K"1 
23°K  x  0.001  Cal0K_1  =  0.023  Cal  g"1 


+  heat  of  fusion  =  0.079  Cal  g"1 
0.1  Cal  g"1 
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Jotes  to  Table  3  -  Continued 


(3.9  g  m"zd_i)(0.1   Cal   g"1)   =  0.39  Cal  m"zd 


-i\  -, 


-2j-l 
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Notes  to  Table  4 

(1)  Mean  plant  biomass  3.6  kg  m        (Whittaker  and  Likens,  1975) 
Energy  storage  3.6  kg  m"2  x  4,250  Cal  kg"1  =  15,300  Cal  m~2 

(2)  Animal  biomass  4.5  g  m"2  x  4.25  Cal  g"1  =  19  Cal  m"2 

(Whittaker  and  Likens,  1975) 

(3)  Energy  storage  in  fire  0.3  Cal  m  d 
Energy  flow  (See  note  10) 

Assume  mean  storage  time  of  1  min  for  forest  fire 

-2  - 1 

0.3  Cal  m  d   x  1  min   0  „  ,n-k   _  ,  -2 
1,400  min  cR 2  x  10   Cal  m 

(4)  Mass  of  soil  (Brady,  1974) 
Mean  bulk  density  approximately  1.3  g  cm-1 

Mean  soil  depth  approximately  50  cm 

Soil  „ass  .  (1xl0''Cm21^00(kg-''":n'"3)(5°Cm)  =  65°  k9  m"2 

(5)  See  note  4  of  Table  1 

(6)  Gross  primary  production  2  x  net  primary  production  (E.  P.  Odum,  1971) 
3.9  Cal  m"2d_1  x  2  =  7.8  Cal  m"2d_1 

(7)  Respiration  =  gross  production  -  net  primary  production  (E.  P.  Odum, 
1971)  7.8  Cal  m~2d_1  -  3.9  Cal  m"2d_1  =  3.9  Cal  m"2d_1 

(8)  Net  primary  production  0.92  g  m~2d-1  (Whittaker  and  Likens,  1975) 
0.92  g  m"2d_1  x  4.25  Cal  g  =  3.9  Cal  m~2d_1 

(9)  Consumption  by  herbivorous  animals   (Whittaker  and  Likens,  1975) 
0.147  x  4.25  Cal  g'1  =  0.62  Cal  m"2d_1 

(10)   Plant  biomass  consumed  by  fire 

Fire  statistics,  U.S.,  1974  --  (Statistical  Abstract  of  the 
United  States,  1975) 
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Notes  to  Table  4  -  Continued 

Area  in  data  sample   5.86  x  1012  m2  for  U.S. 

Number  of  fires      1.45  x  105  in  U.S. 

Area  burned         1.13  x  1010  m2  in  U.S. 

Mean  area  per  fire    7.8  x  10^  m2  for  U.S. 

Plants  consumed  by  fire 

Mean  continental  standing  stock  of  biomass  =  12  kg  m"2 

(Whittaker  and  Likens,  1975) 
Energy  consumed  =  ^||  o"r"^ple  x  Energy  stored  per  unit  area 

(1.3  x  1010m2  yr-1)x(T2.2kg  m-2)x(4.25  x  IP3  Cal  ko"1) 
(5.86  x  10^2  m2)    (365  d  yr1) —  = 

0.3  Cal  m*2d_1 

As  cross  check:  Energy  released  for  forest  fire  3.67  x  109  Cal 

in  U.S.  (Sellers,  1965) 

(7.8  x  10^  m2  per  fire)x(12.2  kg  m"2)x(4.25  x  103  Cal  kg_1)=4  x  109  Cal 

Thus  this  estimates  is  in  good  agreement  with  Sellers,  1965. 

(11)  Transpiration 

338  g  m  2d  l  (Skinner  and  Turekian,  1973) 

(12)  Energy  dissipated  in  fire  (see  note  10) 
0.3  Cal  m" 2d_1 

(13)  Respiration  of  animals 

0.62  Cal  m  d  1   Set  equal  to  consumption  by  first  law  (see  note  9) 

(14)  Animal  biomass  lost  by  fire 

area  burned   -         .  .    .  ,  .  . 

area  sample"  x  Ener9y  stored  in  animal  biomass  per  unit  area 

(1.13  x  1010  m2  yr'1)  ,.  (19  Cal  m'2)   .    ,„-*,.,  -2j-i 
(b.86  x  10-  m*) x  (365  d  yr-M  =  ] '  x  10   Cal  m  ^d 
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Motes  to  Table  4  -  Continued 

(15)  New  soil    (see  Table  5,  note  10) 
0.51  g  n 

(16)  Erosion 


0.51  g  m"2d_1 


0.55  g  m"2d"1  (see  Table  5,  note  10) 

(17)  World  detritus  production  (Schlesinger,  1977): 

(3.75  x  IP10  mt  yr-1)(8.5  Cal  g-1)(l  x  10  g  mt'1).,  „  r  ,  -2.-1 
(5.1  x  1014  mz)  (365  d  yr x) 

(18)  World  detritus  stock  measured  (Schlesinger,  1977)  in  carbon: 

(1.4xl0^mt)(8  5Cal  f)Q   x  106  g  mt"1),^  x  1Qu  Cal  m~2 
5.1  x  1011+  m2 
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Notes  to  Table  5 

(1)  Mass  of  ocean  crust 

Ql2  =  3.61  TlO^  m2     =  K94  x  10?  kg  m"2  (C1ark'   1971) 

(2)  Mass  of  mantle  of  earth 

Ql3  =  5. 11   x°l0l^9m2  =  7'8  x  109  kg  m"2  (Clark'   1971) 

(3)  World  Fossil   Fuel   Reserve  (United  Nations  Statistical  Yearbook, 
1975) 

1974  Known  Resource 

Bituminous  coal     1.07  x  lO*2  mt  x  6  x  10^  Cal  mt"1   =  6.42  x  1013  Cal 
Lignite  coal  3.4  x  10H   mt  x  4  x  106  Cal    mt"1  =  1.36  x  1018  cal 

Petroleum  7.5  x  10*°  "mt  x  1   x  107  Cal  mt-i   =  0.75  x  IQia  Cal 

Gas  5.9  x  1013  m3  x  9  x  103  Cal  m"3=  0.53  x  IQio  Cal 


9.06  x  1013  Cal 


9.06  x  1013  Cal    ,  70   1nu  r  ,  -, 
5.1  x  10^  m-2    =1"78x  10"  Ca1  m  2 

Bituminous  coal  4.30  x  10H  mt  x  6  x  106  Cal  mt"i  =  2.58  x  ]Q18 

Lignite  coal  1.60  x  10H  mt  x  4  x  106  Cal  mt"i  =  0.64  x  1018 

Petroleum  7.3  x  10io  mt  x  1  x  107  Cal  mO  =  0.75  x  10is 

Gas  5.9  x  1013  m3  x  9  x  103  Cal  m"3  =  0.53  x  1013 

Chemical  energy  stored  in  earth's  crust  from  organic  deposits  of 

oil ,  coal ,  gas,  etc. 

Coal  and  lignite  7.78  x  TO18  Cal 

Natural  Gas  0.53  x  1018  Cal 

Petroleum  liquids  0.75  x  TO18  Cal 

Tar  sands  oil  0.44  x  1018  Cal 

Shale  oil  0.27  x  TO18  Cal 
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n    9.77  x  1018  Cal  _  ,  Q  v  lnk   r  ,  -2 
Q  it  =  5j  x  10if  5=2  -  1 -9  x  10"  Cal  m 

(4)  Stress  energy  stored  in  crust  of  earth  generated  by  mean  plate 
motion  of  2  cm  per  year. 

Annual  release  of  seismic  energy  (Richter,  1958) 

2.14  x  1011+  Cal  yr~l 

Assuming  11  yr  average  period  for  earthquakes 

2.15  x  10"  Cal  vr-i  x  11  yr  =  4>6     -a 

5.1  x  1014  mz 

(5)  Mass  continental  crust  of  earth  (Clark,  1971) 

n    1.6  x  1022  kg   ,  n7  „  ln8  \,„     -2 
Ql5  =  1.49  x  10^  m2  =  1  -07  x  108  kg  m 

(6)  Storage  of  energy  in  volcanoes 

9.6  x  1013  Cal  yr~l     Energy  released  by  volcanoes  (Hayakawa,  1976) 
Assume  1  hr  turnover  time 

?cV  ^nfor  V~\*\   ]  /q  itUr^e[  tim6l\  =  2.2  x  10"5  Cal  m"2 
(5.1  x  1014  Cal  nrz)  x  (3.7  x  10*  hr  yr-1) 

(7)  Storage  energy  in  active  earthquake 
Annual  energy  release 

2.15  x  10lt+  Cal  yr"1  (Richter,  1958) 

Assume  average  life  of  earthquake  of  one  minute 

2.15  x  ~\0lk   Cal  yr"1  x  1  min/earthquake  _  R   1n-7  r  ,  -z 
(5.1  x  10^  m2)  '(5.25  x  10^  min  yr~i)     '   8  x  IU   Ulm 

(8)  Mantle  flowing  into  ocean  trench  to  form  new  ocean  crust 

(4.6     x  1013  kg  yr'MdO3  g  kg"1)   =Q  35         -2d"l 
(3.6  x  10^  m2)(365  dyF1)  U-J5  9  m 

Ocean  plates  to  recycle  in  1.5  x  103  yr  (Clark,  1971) 
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(0.7  x  TO22   kg)    (103   g  kg"1) =  n  .,  n     -2,-1 

3.61   x  10^  m2)(1.3  x  108  yr)   x   (365  d  yr~n       U-J0  9  m     a 

(  9)  Ocean  crust  lost  to  mantle 

fi6.?10^'^'"^^9-!'  ■  0-58  9  m*f* 
(3.61  x  TW*   m2)  x  365  d  yr  l  3 

balances  new  crust  formed  in  note  8  and  sediments  in  note  10 

(!^i"  ViSl*  ^V"  Vr1)(1°3  9kq"1)  ■  0-58  9  ^d"1 
(3.61  x  TO*1*  m2)  x  365  d  yr_i  a 

(10)     Global   erosion  driven  by  world  stream  flow  (Clark,  1971] 

3  x  1013  kg  yr-1 =  0  55  a  m"2d_1- 

1.49  x  10^m2x365  d  yr"i       U'"  9  m     a     ' 

however,  sea  is  2.4  times  larger  than  land  therefore  sediments  = 


0.55  g  m"2d"1  _  n  or>  „     -2,1-1 
^p— 2 =  0.23  g  m     d 

(11)  Stress  by  ocean  plates 

Average  plate  velocity  (Clark,  1971) 

2  cm  yr  l 

(12)  Stress  energy  buildup  from  2  cm  average  motion  of  tectonic 
plates  of  the  earth's  crust  (Clark,  1971) 

(13)  Stress  energy  consumed  by  earthquakes  by  conservation  of  energy. 
Stress  energy  consumed  =  uplift   +  earthquakes    +  heat  loss 
Stress  energy  consumed  =  0.05  +  0.005  +  0.1  =  0.155 

(14)  Energy  required  in  uplift  process 

m  =  annual  uplift  7  x  1013  kg  yr-1 

g  =  9.98  x  m  sec"2 

u  _  continental  crust-ccean  crust  _  40  km  -  7  km  _■,  55  x  -jq^  m 
2  2  x  10_:>km  m_l  "  '  * 
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E  =  mgh  =  7  x  1013  kg  yr~l  x  9.88  m  sec"2  x  1.65  x  10u    m 

x  2.39  x  lP-i+  Cal   J"1 
E  =  2.75  x  1015  Cal  yr"1 

c  -  2.75  x  1015  Cal  yr'1  n  nc  r  ,     -2 .-i 

E  =  /t   Ar. inlL     o\7 — 57r-j -.1=  0.0b  Cal   m     d 

(1.49  x  IP14  m2)  x  365  d  yr1 

(15)  Energy  dissipated  in  earthquake  process 

1  x  1015  Cal  yr"1 =  5.4  x  10'3  Cal  m"2d_1 

(5.1  x  IP14  mz)   365  d  yr~^ 

(16)  Global  seismic  energy  released  by  earthquakes     (Richter,  1958) 

2.15  x  IP14  Cal  yr"1 ,  lc   1n-3  -  ,  -2,-1 

5.1  x  10^  m^)  x  365  d  yr-1  =  1J5x  10   Cal  m  d 

(17)  Stress  in  earth's  crust  dissipated  as  heat,  must  exist  by  second 
law  of  thermodynamics  taken  as  twice  uplift  as  a  minimum 
possible  estimate  (see  note  14  for  uplift  calculation) 

0.C5Cal  m"2d_1  x  2  =0.1  Cal  m"2d_1 

(18)  Rate  of  formation  of  new  soil  set  equal  to  erosion  of  soil 

I?  1 10i1»i"q  VCttl  h  ™U  kg~''    ■   c-51  g  ^d"1 

(1.6  x  1014  mz)    (j65  d  yr  l)  s 

(19)  Net  production  of  fossil  fuel 

Working  with  an  estimated  coal  rate  of  3  x  109mtyr_1  (Hubbert, 
1971)  .  (Note:  this  estimate  may  be  low) 

(Sxicgmtyr"1)  (1  x  106  g  mt"1)(6  Cal  g"1) = 

(365  d  yr-i)(5.14x  10^m2)0.88  (coal  %   of  fossil  fuel) 

0.11  Cal  m'V1 
Heat  loss  in  fossil  fuel  production  must  exist  by  second  law  of 
thermodynamics.  Minimum  value  assumed  of  two  times  production. 
0.11  Cal  m'V1  x  2  =0.22  Cal  m"2d_1 
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(20)  Gross  production  of  fossil  fuel 
Gross  =  Net  +  heat  loss 

Gross  =  (0.11  +  0.22)  Cal  m"2d_1  =  0.33  Cal  m"2d-1 

(21)  World  fossil   fuel  consumption 

1937       1.9  x  109  Metric  tons  Coal   equivalent 
(Oxford  Economic  Atlas  of  the  World,   1965) 

(1.91   x  109  T  yr-1)(l   x  106   g  T"1)    (5.5  Cal   g~M       n  n,,  r  ,     -2  .-i 
(5.1   x  101^  m   (365  d  yr"1) =  °'°56  Cal  m     d 

1960       4.31   x  109  Metric  tons  Coal   equivalent 

(Oxford  Economic  Atlas  of  the  World,  1965) 

(4.13  x  109  T  yr-1)(l   x  106   g  T"1)(5.5  Cal   g'1)   _  n  197  r  .     -2,-i 
(5.1   x  10^  m2)    (365  dyH)  °'127  Cal   m     d 

1971       7.14  x  109  Metric  tons  Coal   Equivalent 

(United  Nations  Statistical  Year  Book,  1975,  1976) 

(7J4x  109  T  yr"1)    (1   x  106   g  T~1)(5.5  Cal   g"1)   _  n  „   .  ,     -2  .-i 
(5.1   x  101^1    (365  d  yr-1) "  °'21   Cal   m     d 

1974       7.95  x  109  yr~l    (United  Nations  Statistical   Year  Book,   1975) 

(7.95  x  109  T  yr~l){-{   x  IP6   g  T"x)(5.5  Cal    q"1)       n  „  r  ,      -2,-i 
(5.1   x  10^iP)f365  dyr-1)  =  °'23  Cal   m     d 

Fossil   fuel   consumption   (1970)  (Hubbert,   1971) 

(4  x  1013KWH  yr-1)(859  Cal    KWH"1)      _  n  ln  r  ,      -2  ,-i 
(5.1   x  10^m2)(365  d  yr~i)  "  °'18  Cal  m     d 

(22)  Global   energy  released  annually  by  volcanoes 

(5jVfoffixy365  d  yr-i   =  5J   x  10"U  Cal  ^d"1  (Hayakawa,   1976) 

(23)  Heat  dissipated  in  volcanic  process  by  second  law  must  exist 
estimated  twice  volcanic  release  in  note  20 

(5  x  10_Lf  Cal   m"2d_1)   x  2  =  1   x  10"3  Cal  m"2d_1 
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(24)  Energy  consumed  by  volcanoes  by  law  of  energy  consumption 
(1  x  10"3  +  5  x  10_tt)  Cal  m~2d_1  =  1.5  x  10"3  Cal  m~2d-1 

(25)  Deep  heat 

1.9  x  IP17  Cal  yr"1 =  1   Cal   rrfV1  (Hayakawa,   1976) 

(5.1   x  10il+  mz)   x   (365  d  yr~i) 

(26)  Energy  storage  in  floods 

Assume  1   hr  turnover  time  (see  note  19  for  consumption) 

consumption         m  0.04  jjI/^C   x1   hr  =  1 .7  Cal  nf3 
Wli      turnover  time  24  hr  d-1 

(27)  Energy  consumed  by  floods 

Mass  10%  of  continental   rain   (9.9  x  1015  kg) (0.1)   =  9.9  x  1014  kg 
Elevation  1,000  m 

E  =  mgh  =  9.9  x  1011+  kg  g"1  x  9.9  m  sec"2  x  1,000  m  x  2.3  x 
10_1+  Kg"1  m"2sec     Cal   = 

E  ■  ?;3  *  ^ni^r"1/^  „       -n  "  0-04  Cal  m^d"1 
(1 .49  x  10iH  mz)  x  (365  d  yr  ±) 

(28)  Flow  of  water  in  world's  rivers  to  produce  erosion  of  continental 
crust  to  ocean  crust 

°}l9cX   ]?niu9  JrrrogE  a =tx    =  753  g  m~2d_1  precipitation  on  land 

(3.6  x  10li+  mz)  (365  d  yr  1)  3 

(Skinner  and  Turekian,  1973) 

(29)  Erosion  to  sea  (Skinner  and  Turekian,  1973) 

3  x  1013  kg  yr~l   _n  „  n   m-2,H_1 
(3.6  x  101If  m2)       3 

(30)  Energy  dissipated  in  floods  set  equal  to  energy  consumption  for 
conservation  of  energy  (see  note  19) 

0.04  Cal  m"2d_1 
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(1)  Energy  storage  of  cities 
8%   depreciation  rate 

2%   fire  loss 
1%  war  loss 
0.5%  disaster  loss 
11. 5% total  loss 

-2  -1 

consumption    0.45  Cal  m  d   x  365d   ,  nnn   r,n  -2 
depreciation  =  07T15 =  ]  '400  Cal  m 

(2)  Energy  consumed  by  cities  (Hubbert,  1971) 
Consumption  of  forest  fuel;  1970 

0.18  Cal  rn'V1  (See  note  21,  Table  4) 

Consumption  of  biomass,  1970  (Whittaker  and  Likens,  1975) 

0.27  Cal   nfV1 

-2-1 

Total   consumption  0.45  Cal   m     d 

(3)  Air  pollution  from  cities 

Release  of  1/4.25  g  Cal"1 

0.45  Cal  _  n  ,   n     -2,-1 

4.25  Cal    g-r-  0A    9  m     d 

(4)  Energy  dissipated  by  cities  equal  consumption  by  conservation  of 
energy  less  growth  (see  note  7) 

0.45  Cal  nfV1 

(5)  Operation  of  urban  structure  (Solow,  1954) 
Estimated  depreciation  and  obsolescence  for  capital  assets 

8%  of  stock  per  year 


(1.400C«lm-i)(Q.08),yr-1.0-31  M   ,-td-. 

36o  d  yr  L 
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(  6)  Fire  loss  estimated  from  U.S.  data  (United  States  Statistical 

Abstract,  1975) 

Total  urban  loss  for  1965  due  to  fire  2.5  x  109  yr~l   _ 
National  wealth  14.5  x  10*6 

(1>400C^mH2)(°:?17)yr"1  =  0-065  Cal  m-V1 
365  d  yr  l 

(  7)  Energy  dissipated  by  fire  set  equal  to  consumption  in  note  6 

0.065  Cal  m"2d_1 


(  8)  Energy  stored  in  urban  fire.  Assume  1  hr  turnover  time 

0.065  Cal  m"2d"1  x  1  hr  _  0   7  „  ,n-3  r  ,  -2 

0/1  r — jtt =  2.7  x  10   Cal  m 

24  hr  d  l 

(  9)  Energy  consumed  by  war  (Statistical  Abstract  of  U.S.,  1975) 

World  Military  Expenditures  (1973)  2.7  x  1011  $  yr"1 

World  Per  Capita  GNP  in  1968  dollars   $600  (Meadows  et  aj_. ,  1974) 

World  Population  3.78  x  109  people  (Statistical  Abstract  of  U.S., 

1975) 

2.7  s  IP11  $  yr n  ,, 

(660  $  people-M(3.78  x  10"7  people)  '  u'" 

War  consumption  of  energy  =0.45  Cal  m~2d  l   x  0. 11  =0.05  Cal  m  2d  l 

-2  -1 

(10)  War  dissipation  set  equal  to  consumption  of  0.05  Cal  m  d   from 
note  9. 

(11)  Energy  stored  in  war  equipment  assuming  5-yr  lifetime  of  war 
equipment 

~2,|~1  -3CC  A     w«.~l 

=  90  Cal  m 


0.05  Cal  m"2d"1  365  d  yr"1  _  qn  -,  -2 


0.2 

(12)  Disaster  losses  for  U.S.  in  1964  dollars  (Dacy  and  Kunreuther, 
1969) 
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(1950-1965  averages)       Fraction  of  US  GNP  Fraction  of  wealth 
Storms     2.3  x  108  $  yr                  0.005  0.0016 

Floods     2.1  x  108  $  yr~l  0.0045  0.0015 

Earthquakes  1.9  x  108  $  yr'1  0.014  0.0014 


6.3  x  10  0.014  0.0045 

Where  GNP  =  4.6  x  1010  $  yr'1 
National  wealth  =  14.5  x  1010  $ 
Energy  loss  due  disaster 
(0.45  Cal  m"2d_1  consumption) (0.014  =  6.3  10"3  Cal  m~2d_1) 
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(1)  Calculations  for  storm     loss  based  on  U.S.   data 

1950-1956  Average  annual   U.S.   loss  (Dacy  and  Kunreuther,   1969] 

Storms 

Hurricanes  1.7  x  108  $yr    l 

Tornados 0.5  x  108  $yr~  1 

Strom  total  2.2  x  108  Syr'1 

2.2  x  108  Syr"1  (U.S.   storm  loss)   _   ,  B       in-3 
4.6  x  10B     $yr"MU.S.   GNP)  "  4,°  X   IU 

Energy  embodied  in  loss  due  to  storms  calculated  based  on  energy  con- 
sumption of  cities 

(0.45  Cal   m"2d_1)  x  (4.8  x  10"3)   =  2.1   x  10"3  Cal  m"2d_1 
Storm  energy  to  cities  0.014  Cal   m~  2d_1     (Note  4) 
Disaster-amplifier  ratio  in  Calories  heat  equivalent 

Energy  Disordered  _  2.1   x  10~3  Cal  m  2d   1  .  n  u  r^-,   r3l-i 
Storm  flow 0.014  Cal   nr-cr1  0J5  ~aI   Cal 

Disaster-amplifier  ratio  in  embodied  solar  energy 

D   .  .      .     (2.1   x  IP"3  Cal  m"2d"1city)(1.5  x  10^)   Cal   Cal^city"1 
Kat1°  1n  0.014  Cal   m-2d-*  fl.7  x  1 0 3 )   Cal   Cal"i 

1.4  Solar  Cal   per  Solar  Cal 

(2)  Calculation  for  floods  based  on  U.S.   data 

1950-1965  average  annual   flood  loss   (Dacy  and  Kunreuther,   1959) 

Floods       2.1   x  10s  Syr"1 

2.1   x  1Q8  $yr-1   (Flood  loss)   _  .  ,  „  ln-3 
4.6  x  1010  |yr-l   (GNP)  4-5  x   IU 


Flood  energy  loss   (0.45  Cal  m"2d-1)(4.6  x  10"3) 

=  2.1   x  10"3  Cal   m"2d"1 
Flood  energy  to  cities       1.4  x  10"4   (Note  4) 
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Disaster-amplifier  ratio 

H   X  }nl   r9]  m"^"l  ■  15  Cal  Car1 
1 .4  x  10-4  Cal  m-zd-1 

Disaster-amplifier  ratio  in  embodied  solar  energy 

(2.1  x  10"3  Cal  m"^"1)  1.5  x  ~\Qk     _  -,  ?  ,.  ,   p  ,     r  ■>       r  , 
(1.4  x  10-^  Cal  m-2d-i)(1.7  x  10^)  "  ''3  Solar  Cal  per  So1ar  Cal 

(3)  Earthquake  loss  based  on  U.S.  data  (Dacy  and  Kunreuther,  1969) 

1950-1965  averages 

Earthquake  loss  1.9  x  10  Syr"1  _  a   -,   ln-3 
GNP  4.6  x  10  $yr-i  "    * 

Energy  loss  (0.45  Cal  m~2d_1)  (4.1  x  10"3)  =  1.85  x  10"3  Cal  m"2d-1 

Seismic  energy  to  cities  3.4  x  10~6  Cal  m~2d_1  (Note  4) 

Disaster-amplifier  ratio  _  1 .85  x  10"3  Cal  m  2d  1_   53g  Ca-j  ral_1 

3.4  x  10"6  Cal  m~2d-1 

Disaster-amplifier  ratio  in  embodied  solar  energy 

(1.85  x  10~3  Cal  m'2d~l   )(1.5  x  ]Qk)   u   ,  7  -  ,   -  ,     s  ,   -  , 
(3.4  x  10-*»  Cal  m-2d-l)  (4.7  x  10<>)    K7  So1ar  Cal  per  b0iar  Ul 

(4)  Disaster  energy  available  tocities  and  urban  places 
World  mean  urban  plus  industrial  area  per  person 

500  m2  (Meadows  et  al_. ,  1974) 

Fraction  of  world  surface  area  urban 

500  m2  per  person  x  3.5  x  109  people  _  n  nn,« 
5.1  x  10^  m2  u'uu^ 

Storm  energy  absorbed      (Table  2,  note  18) 

4  Cal  m2  x  0.0034  =  0.014  Cal  m'2d_1 

Flood  energy  absorbed  (Table  5,  note  17) 

0.04  x  0.0034  =  1.4  x  10"u  Cal  m~2d-1 

Seismic  energy  absorbed  (Table  5,  note  16) 

0.0C1  x  0.0034  =  3.4  x  10'6  Cal  m"2d"1 
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Documentation  of  Order-Disorder  Model  of  Gamma  Ray  Radiation  Stress 

(1)  Source  =  1  x  lO^Ci  of  137Cs  which  has  0.662  MeV  Y  97% 

X  =  (1  x  ~\Qk   Ci)  (3.7  x  1010  dissec-^Cf"1) (0.662  x  0.92)  MeV  dis-1 

(1.6  x  10"5  erg  MeV-1) 
X  =  3.6  x  108  erg  sec"1  =  3.11  x  1013  ergs  d_1  =  744.5  Cal  d_1 

(2)  Energy  absorbed  by  forest  air 

biomass  of  rain  forest  =  42,450  g  m~^  -  canopy  height  =  20  m 

(Odum,  1970c) 

42,950  g  m~2         _  n  nn91    _       -3 

■^rr tttc — 5-t      0.00<£     g  cm 

20  m  x  10^cm3m3  3 

density  of  air  is  0.0012  g  cm"3  thus  total   density  is  0.0033  g  cm"3 

with  path  of  25  m  (2,500  cm)   and  absorption  of  0.02  cm     g 

K     =   (2,000  cm)(0.0033  g  cm'3)(0.02  cm2g"i  =  0.132 

(3)  Thus  (0.132)(744.5  Cal  d"1)  =  98.27  Cal  d"1  absorbed 
Therefore  744.5  Cal  -  98.27  Cal  d"1  =  646.2  Cal  d'1 

(4)  Inverse  square  effect: 

4*r2  =  4  x  3.14  x  (2,000  cm2)  =  5.02  x  107  cm2 

(5)  Exposure  per  cm2: 

ffi&lkC  -  I-*  X10-*  Cal  art"' 

(6)  (1.28  x  10"5  Cal  cm"2  d_1)(43,950  g  m"2)(l  cm3g_1)  =  0.552  Cal"2d 

(7)  Absorbed  dose: 

(0.02  cm  g_1)(1.28  x  10"5  Cal   cm"2d_1)   =  2.56  x  10"7  Cal^d"1 

2.56  x  10"7  Cal   g  d"1  =  2.56  x  10"7  Cal   g'1  day  x  4.184  x  1010ergs  Cal-1 

=  1.07  x  1014  ergs"1 


1.071  x  101*  ergs  0.01  rad  erg  ld   1  =   107.1  rad  d 
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(  8)  Total  absorbed  dose/m2: 

(43,950  g  m"2)(2.56  x  10~7Cal   g_1d_1)   =  1.099  x  10"2Cal  m2  d"1 
(   9)   Solar  insulation  daily  mean  (Odum  et  al_. ,   1970) 

3,830  Cal  m'V1 

(10)  Gross  primary  production  of  rain  forest  at  El  Verde  (Odum,  1970c) 
in  g  organic  matter 

32  g  m"2 

(11)  Gross  primary  production  of  a  rain  forest         (Odum,  1970c) 
in  energy  units 

(32  g  m"2d_1)  (4.25  Cal  g"1)  =  145  Cal  nfV1 

(12)  Biomass  of  rain  forest  in  g  organic  matter        (Odum,  1970c) 
42,450  g  m"2 

(13)  Litter  biomass  of  forest  (Odum,  1970c) 
14,630  g  m"2 

(14)  Respiration  assumed,  based  on  reported  system  response  (Odum,  1970c) 
22  g  m"2d_1 

(15)  Radiation  induced  disorder,  assumed  based  on  reported  system 
response  (Odum,  1970c) 
10  g  m  2d  1   at  20  m  from  source 
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Johnstown  Flood 

(1)  Rainfall  into  South  Fork  basin 

Rainfall  6-8  in/12  hr  (Engineering  Society  of  Western  Pennsylvania, 

1889) 

using  7  in/12  hr 

7  in/12  hr  x  0.0254  m  in"1  x  2.77  x  10_lf  hr  sec"1  =  4.1  x  lO^msec"1 

Given  a  total  basin  area  of  60  sq  mi 

60  sq  mi  x  2.59  x  106  m2  mi"2  =  1.554  x  108m2 

Rainfall  into  basin  therefore  is 

Rainfall  x  land  area 

4.1  x  10'5  m  sec"1  x  1.554  x  loV  =  637  m3  sec"1 

3     - 1 

Shappee  (1940)  calculated  591  m  sec 

(2)  Flow  into  reservoir 

rate  in  =  increased  +  rate  out 

3-1  3-1         3-1 

rate  in  =  199  m  sec   +  169.0  m  sec   =  368  m  sec 
second  reference  gives  283  m  sec   (Shappee,  1940) 

(3)  Water  loss  from  reservoir  -  designed  rate 

3,500  ft3  sec"1  (Engineering  Society  of  Western  Pennsylvania,  1889) 

3,500  ft3  sec-1  x  2.83  x  10~2m3  ft"3  =  9,910  m3  sec"1 

Spillway  max,  (Shappee,  1940)  -  6,000  ft3  sec"1  =  169  m3  sec"1 

(4)  Water  over  top  of  dam,  center 
4,000  ft3  sec"1  (Shappee,  1940) 

4,000  ft3  sec"1  x  2.83  x  10"V  ft3  =  133  m3  sec'1 

(5)  Turbulence  from  overflow  (note  4  -  note  6) 
133  m  sec""  -  99  m  sec""  =  34  m  sec"1 
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(Shappee,  1940) 

(Engineering  News,  1889) 

=  1.35  x  107m3 

(Engineering  News,  1889) 

(Engineering  News,  1889a) 

(Engineering  News,  1889a] 


(  6)  Spillway  10  m  long 

99  m3  sec"1 
(  7)  Volume  of  South  Fork  reservoir 

4.8  x  108    ft3 

4.8  x  108         ft3x  2.83  x  10"2  m2    ft 
(  8)     Spillway  storage  capacity 

1,690  m3 
(  9)     Land  area  of  South  Fork  Valley 

1.55  x  103  m2 

(10)  Slope 
Land  covered  by  reservoir 

Calculated  by  dividing  reservoir  capacity  by  depth. 

Reservoir  depth   62  ft 

62  ft  x  0.3049  m  ft_1=  18.9  m 

(11)  Tangible  assets  of  Johnstown,  Pennsylvania  (Historical  Statistics 
of  the  United  States,  1976) 

1889  U.S.  per  capita  tangible  assets 

1.795  x  TO11  (series  F447)  =  $2  365 

7.59  x  107  people  ^'3bb   per  perSOn 

Johnstown  Population  (1889)  10,000  (McCullough,  1968) 

Johnstown  Tangible  Assets  (1889) 

($2,365/person)  (10,000  people)  =  2.4  x  107 

Energy  value  of  1889  dollar  =  Gross  Energy  Consumption  (note  10) 

Gross  Product      (note  10) 

ikM  10?r?a1  J\l)    (365  d  ^  -  1-1  *  1Q5  Cal  S"1 

($4.ol  x  106yr-1) 

Energy  value  of  Johnston   ($2.4  x  107)    (1.1   x  105  Cal   $_1)   = 
2.6  x  1012  Cal 
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(12)  Land  area  of  Johnstown 

1.10  x  106  m2   (USGS  Map,  Johnstown,  Pennsylvania,  1964) 
(McCullough,  1968) 

(13)  Natural  resources 

Estimated  forest  area, 

2  x  108  m2     (USGS  Map,  Johnstown,  Pennsylvania,  1964) 

Mean  biomass  standing  stock  for  deciduous  temperate  forest 

30  kg  m"2      (Whittaker  and  Likens,  1975) 

Total  standing  stock 

(2  x  108m2)(30  kg  nf2)(4.25  x  103  Cal  kg)  =  2.55  x  1013  Cal 

(14)  Coal  reserve,  Johnstown 

Est.   1  x  lOS  1889  annual   consumption 

1   x  109  Cal   d"1  x  365dyr4     x  10 V  =  3  x  1015Cal 

(15)  Net  primary  production  of  natural   system  (Whittaker  and  Likens, 

1973) 

Land  area  2  x  108  m2 

Mean  net  primary  production  for  deciduous  temperate  forest 

(1,200  g  yr~V2)(4.25  Cal  g-1)(2  x  108m2).  9  8   1Q9Cal  d-i 
365  d  yr~l 

(16)  Resource  consumption 

US  consumption  of  timber  products  1900 

Historical  Statistics  of  United  States  (1975) 

(7.14  x  109ft3yr-1)(.0283  m"3  ft3)  (8  x  T(gq  m3)  (4.25  Cal  or1) 
(365  d  yr-   )   7.59  x  HP  people 

2.4  x  TO*  Cal  d"1  person"1 

Food  energy  U.S.  per  capita  1909  (Historical  Statistics  of  U.S., 

1975) 
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3.5  x  104  Cal  d"1  person"1 

Firewood  44%  value  of  forest  products  for  1900 

Estimated  consumption 

1.  x  104  Cal  d"1  person 
Total  renewable  resource  consumption 
Timber  products  2.4  x  104  Cal  d-1  person 
Fuel  wood      1  x  10h   Cal  d"1  person 
Food         0.3  x  104  Cal  d"1  person 

3.8  x  10k   Cal  d-1person 
Total  consumption  estimate  for  Johnstown  1889 

(3.8  x  10**  Cal  d'1)  x  (10,000  people)  =  3.8  x  108  Cal  d"1 

(17)  Respiration  of  natural  system 
Production  -  Consumption 

2.8  x  109  Cal  d"1  -  3.8  x  108  Cal  d"1  =  2.4  x  109  Cal  d"1 

(18)  Land  consumption 

Based  on  50%  land  in  valley  occupied,  low  density  in  Johnstown 

(McCullough,  1968) 

(1.1  x  106  m2)  (.5  occupied)  =  ,   ,Q-7  miCal-i 
1.8  x  10^2  cal         J  x  iu   mui 

(19)  Johnstown  coal  consumption  (Statistical /Abstract  of  the  United  States, 1971 
Based  on  1920  data 

1.5  x  108  BTU/personx  0.252  Cal  BTlf1  =  3.8  x  10  Cal  yr~l   person 

3'8  365°d  y?-!^"1  =  10'000  people  =  1'°  x  1Q9  Cal  d_1 

(20)  Energy  value  of  gross  regional  production  (1889) 
Natural  (note  6)  =  0.38  x  109  Cal  d"1 
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Fossil  Fuel  =   1.0  x  109  Cal  d"1 
1.4  x  10y  Cal  d"1 

Gross  National  Product  _  3.5  x  1010       Mmn  ^  n   \ 
Population 7.59  x  10?       (1929  dollars) 

($461  person-1)O0,000)  =  $4,610,000 

(21)  Net  regional  product       (Historial  Statistics  of  U.S.,  1975) 

Net  National  Product   $3.1  x  IP10   *.nQ      -i 
Population =  7.89  x  109~  =  $408  person 

$408  person  x  1,000  people  =  $4,080,000 

Net  regional  product    _  4,080,000  _n  ao 
Gross  regional  product    4,610,000  "u-ao 

Net  Johnstown  consumption 

(1.4  x  109  Cal  d_1)U88)  =  1.2  x  109  Cal  d"1 

(22)  Depreciation  8%  (Solow,  1965) 
Consumption  10% 

(i-8x3Wd*ffy8,'r)  =  8-9  x  io8  ca' d_l 

(23)  Catastrophic  destruction  by  flood 

Damage  $1.7  x  107  (Gilbert,   1977) 

($1.7  x  107)(1.1   x  105  Calhe     $_1)   =  1.9  x  1012Calhe 

To  convert  to  solar  calories 

Johnstown  energy  consumption  in  solar  equivalents 

Biomass  consumption  (note  15) 

Plant  biomass  quality  factor  (Table  8) 

(3.8  x  103  Calhe     0(1.7  x  103  Cal$e  Cal^"1)   =  7  x  10uCalsed_1 

Fossil   fuel   consumption  (note  6) 

Fossil   fuel   quality  factor  (Table  8) 

(1.   x  109  Cal   d_1)(6.4  x  104  Cal   „  C,    _1)   =  6.4  x  1013Cal     d"1 
\  /  v  se     ne  se 
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Quality  factor  for  Johnstown  =  6.4  x  1013  Cal  d"1  solar 

1.4  x  10y  Cal  d"1 

4.6  x  10"  Cal  Cal"1 

Solar  energy  embodied  in  flood  destruction 

Solar  energy  lost 

(19  x  1012  Cal)  (4.6  x  104  Calse  Cal^"1)  =  8.6  x  1016  Cal  Solar 

(24)  Gravitational  potential  energy  of  water  in  South  Fork 
Reference  Johnstown 

South  Fork  elevation  1,618  ft  (McCullough,  1968) 

Johnstown  elevation  1,120  ft  (McCullough,  1968) 

E  =  mgh  =  (1.36  x  1010  kg)(9.8  m  sec  2)(1 ,618-1 ,120  ft)  x  .304  m  ft"1 

E  =  2.01  x  1013  joules  x  4.186  x  10"3 Cal joules"1'  3.43  x  1010  Cal 

Potential  energy  per  m~3  in  South  Fork  reservoir  with  respect 

to  Johnstown 

E  =  mgh 

E/m3  =  (1,000  kg  nf  3)(9.8  m  sec"2)(l,618  ft  -1,120  ft)  x  0.304  m  ft"1 

E/m3  =  1.65  x  105  joules  m3  x  4.186  x  10"3  =  6.94  x  103  Cal  m"3 

Slope  of  Conemaugh  Valley 

South  Fork  of  Johnstown 

-,     elevation   (1,618  ft-1,120  ft)  0.304  m  ft"1 
5lope   distance  "  16  mi  x  1,390  m  mi"1 

151.4  m 
25744  m 

Slope   5.88  x  10"3m  m"1 

(25)  Wave  energy 

Wave  height  33  f t  x  0.304  m  ft"1  =  10.  m 
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,-i 


,.     ,   ..   16  mi  x  1,609  m  mi     -,  cc       „„  -l 
Wave  velocity  cg  m^n  z   en  fQ,-  m,-»-i  =  7-66  m  sec 
J   56  mm  x  60  sec  mm  x 

Wave  width    500  m  (from  USGS  map) 

Assuming  all  water  in  reservoir  in  rectangular  wave 

■     _   V  1 .3   x   IP7  m3    _   g   r    „    ln3   m 

L    =    T       =    nrCn Tn =    2.6    X    10°    m 

A    500  m  x  10m 
Linear  kinetic  energy  per  m3  of  water 
E  =  1/2  m  V2  =  1/2  (1,000  kg  m"3)(7.66  m  sec"1}2 

E  =  2.933  x  10  ^J!_  m     =   2.933  x  10^  joules  x  4.186  x  10"3 
sec 

Cal  joule"1  m"3  122.81  Cal  m"3 

Turbulence  wave  energy 

125  H2  ergs  cm"2  (Weyl ,  1970) 

where  H  is  in  cm 

H  =  30  ft  x  30.48  cm  ft"1   =   1,040  cm 

Q.    =  125  (1,040)2  ergs  cm"2  x  10,000  cm2  m"2  x  2.39  x  10"11  Cal   erg" 


Q     =  32.3  Cal  m 
w 


f2 

Total  wave  energy  (Engineering  News,  1889b) 

K  (32.3  Cal  m"2)(1.36  x  107  m2)  =  4.4  x  107  Cal 
^26)  Total  wave  energy 
Linear  +  turbulent 
Linear 

(122  Cal  m"3)(1.36  x  107  m3)  =  1 .67  x  109  Cal 
Turbulent 

(32  Cal  m"2)(1.36  x  10s  m2)  =  4.4  x  107  Cal 
Total  flood  energy         =1.7  x  109  Cal 
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Guatemala  Earthquake 

(1)  Strain  in  bars  -  relating  the  230  km  rupture  and  the  108  cm 
average  fault  displacement  (Plafker,  1976)  of  the  Motagua  Fault 
during  the  1976  Guatemala  Earthquake  on  Dieterich  (1974) 
graphical  representation  of  scaling  laws  of  seismic  energy 
yield  three  bar  strain  release.  This  point  also  matches  ob- 
served seismic  energy  release  on  Dieterich' s  monograph. 

ET  =  strain  per  unit  volume  x  volume  (see  note  2) 

Ey  =  3  bar  x  4.83  x  1016  cm3 

ET  =  3  x  105  dynes  car2  x  4.83  x  1015cm3 

ET  =  1.3  x  1023  dynes  cm  =  ergs 

ET  =  2.62  x  1012  Cal 

(2)  Volume  of  material  under  strain  =  V  =  Iwh 

lenght  =  230  km,  width  =  3  m,  height  =  7  km.  Plafker  (1976) 
gives  a  volume  of  4.83  x  1015  cm3. 

(3)  Observed  value  of  stored  energy  at  time  of  rupture  based  on 
back  calculation 

2.62  x  1012  Cal  (Espinosa  et  al_. ,  1976) 

(4)  Strain  energy  source,  westernly  nation  of  North  American  Plate 
2.12  cm  yr"1  x  3.17"8  yr   sec"1  =  6.72  x  10"8  cm  sec'X Jordan,  1975) 

(5)  Strain  energy  source,  Caribbean  Plate 

0  cm  yr"1  =  0  cm  sec"1  (Jordan,  1975) 

(6)  Annual  buildup  of  strain  energy  = 

total  buildup  E 

pp-T +  depreciation  (7) 
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(  7)  Depreciation  on  annual  energy  lost  through  minor  fault  slips; 
100  small  shocks  of  estimated  average  magnitude  of  4  on 
Richter  scale  report  (Harlow,  1976) 

R  =  9.4  +  2.14  or  -  0.054M2   (Gutenberg  and  Richter,  1969) 
Es  =  2.25  x  1019  ergs 

Depreciation  =   100  yr~l   x  2.55  x  1019ercs  =  2.25  x  1021ergs  yr~l 
Depreciation  =  2.25  x  1021  ergs  yr"1  x  2.39  x  1011  Cal  erg"1  = 
3.17  x  10"3  yr   sec"1  x  5.37  x  1010  Cal  yr"1  =  1.76  x  103  Cal  sec"1 
(  8)  Total  seismic  energy  released  in  main  event 
M  =  7.5  (Espinosa  et  a]_. ,  1976) 
where  M  is  the  magnitude  on  Richter  scale 


I2 

\2.Z    _„_,.  -  c  /i/i  „  mil 


E     =  9.4  +  2.14M2  (Gutenberg  and  Richter,   1969) 


Es  =  5.44  x  1011  Cal 


Es  =  2.28  x  10z~  ergs  =  5.44  x  1011  Cal 
Energy  lost  in  plate  motion 
(  9)  ET  =  Es  +  W  (Dieterich,  1974) 

where  Ey  =  2.62  x  1012  Cal 
I11 

(10)  W  =  2.07  x  1011  Cal 

(11)  Material  strained 

4'832xQ^16rcm3  =  2.4  x  1014  cm3  yr'1   x  3.17  x  10"8  yr   sec-1= 

7.61  x  106  cm  sec"1 

(12)  Momentary  storage  M  earthquake  pulse 
1.3  x  106  Cal 

(13)  Calculations  for  typical  house 
3  =  6  m,  w  =  4.25,  h  =  2.5  m,  wall  thickness  =  0.3  m, 
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density  of  adobe  =  1.538  kg  m"3,  long  wall  volume  =  2  (6x2. 5x. 3) 

+  2(4. 25x2. 5x. 3)  =  15.37  m3 
mass  =  15.37  x  1.538  kg  m-3  =  2.4  x  104  kg 

average  h  =  -j- 

PE  =  mgh  =  2.4  x  10^  x  9.8  x  1.25  =  2.9  x  105  J 

PE  walls  =  2.9  x  1012  ergs 

Roof: 

mass  =  40  m2x  0.05  m  thick  x  1,922  kg  m3  =  3,840  kg 

PE  =  mgh  =  3,840  x  9.8  x  3  =  1.13  x  103  J 

PE  house  =  4.07  x  105  J   =  4.07  x  1012  ergs  =  97.2  Cal 

Area  of  contact  =  thickness  x  p  =06  m2 

PE  house/m2  =  97'2Ja1  =  15  Cal  m"2 
6  mz 

Dissipation  of  energy 

Volume  of  earth  under  house  =  6.25  m2  x  7,000  m 

=  4.37  x  TO*  m3 

density  of  plate  =  2,800  kg  m~3 

mass  of  plate  =  1.2250  x  103  kg 

mass  of  house  =  2.42  x  }Qh   kg 

fraction  house  =  1.9  x  10" 

(14)  Mass  of  houses  in  El  Progreso 

(2,208  houses)(3,840  kg  house"1)  =  8.4  x  10s  kg 

(15)  Seismic  energy  transmitted  to  El  Progreso 

1.17  x  106  Cal  (see  Table  12)x  166  sec-l  =  ]gs   x  1Q4  Cal  sec-! 
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(16)  Energy  dissipated  by  structure  of  El  Progreso  is  given  by  number 
of  structures  x  area  of  structure  x  energy  per  unit  area. 

2,208  (houses  assuming  5  people  per  house) x  6  m2  (per  house 
foundation)  x  3.2  x  10"3Cal  m"2  x  1.66  x  10"2  sec_1=  0.7  Cal  sec-1(10) 

(17)  Depreciation  is  calculated  by  total  storage  of  wave 
energy/time  of  half  life 

5.44  x  1011  Cal  x  1.66  x  10"2  sec"1    (Espinosa  et  a]_.,  1976) 
=  9.04  x  109  Cal  sec"1 

(18)  Construction  and  maintenance  of  buildings  in  El  Progreso 
10%  annual  construction  maintenance  of  stock  assumed 
2.3  x  10  Cal  yr~l   (gravitational )x  3.17  x  10"8yr  sec"1 
=  7.29  x  10"4  Cal  sec"1 

(19)  Damage  95% 

2.3  x  10  Cal  (per  event)  1.66  x  10"2  sec"1  =  3.83  x  103  Cal  sec"1 

(20)  Depreciation  of  El  Progreso 

Based  on  construction  rate  plus  5%  growth 

2.14  x  10  Cal  yr"1  x  3.17  x  10"8  yr   sec'1  =  6.78  x  10_l+  Cal  sec"1 

(21)  Disorder  material  (see  note  13) 
Mass  of  town  times  95%  damage 

5.8  x  10  kg  (per  event)  x  1.66  x  10~2  sec'1  =  9.66  x  10k    kg  sec-1 

(22)  Construction  material 

Mass  of  town  times  construction  rate 

6.13  x  105  kg  yr"1  x  3.17  x  10"8yr  sec"1  =  1.94  x  10'2kg  sec"1 


APPENDIX  II 
DYNAMO  PROGRAMS 
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GAMMA  RADATIQN  MODEL      8/16/78 

*  GAMMA     RADATIQN     MODEL 

L  QI. K=Q1 .J+DT*<R1 . JK-R2. JK-R3 .JK) 

N  01=42450 

A  G2.K=57030-Q1 .K 

R  R1.KL=K  l*J*Ql.K*Q2.K/(  1*-K0*G1  .  K*Q2  .K) 

C  J=2300 

C  K0=3»22E-9 

C  K1=7E-11 

*  RADATION    SECTION 

A  R0.K=X-(R*1E2*6.6E-5)*X 

C  X=744,5 

R  R2.KL=CLIP(RX.K, 0.93. TIME. K) 

A  RX.K=K2*R0.K*Gl.K/(  ( 1  +  K4*Q 1 . K ) * ( 4*3 .  1  4 *R* R* I E4 )  ) 

C  K2=1.3  5E2 

C  K4  =  l .436E-4 

R  R3.KL=K3*Q1.K 

C  K3=  0.0  0  05  183 

C  R=2  0 

PRINT    Q1.Q2.R1.R2.R3.J.R0 

PLOT  01  =  1  ,G2=2(  0.  75  00  0)/Rl=*.R2  =  %,RJ=-«- 

SPEC  DT=. 1/LENGTH=365/PLTPER=7.3/PRTPER=73 

RUN 
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3/01/73 


JOHNSTC 
FAILURE 
Q  l»K=Qh 
Ql=li 36E 
PlaKL=(K 
KO=l,57E 
J=7E-6 
K l=0,  414 
R  2»KL=K2 
K2=7.27E 
°  3oKL  =  K3 
K3=9o87E 

Q3:>K  =  Q3:> 

Q 3=990 

R4jKL=K4 
K4=2, 52  E 
R  6-><L=K6 
K6=0,  l 
Q2-,  K=L-K 
L=i,55E3 
<7=C,  0  52 
MODEL  OF 
Q6J<=063 
•0  6=06  1 
Q6I=2,5E 
R9jKL=K8 
R=l->  *E  12 
K8=2E-2 
P9*KL=KO 
X9=2,  5  5E 
RlOaKL=< 
<1 :=9,6E 
Q7,<=Q  7, 
07=3 El  5 
=>  123  <L=K 
<  12=6,  73 
C4->K  =  QA> 
Q4=l) 8E1 
R  13,  KL=K 
K13=1,35 
P  14,KL=K 
K  14=9,4.5 
P  15s  KL=< 
K  15=4  ,9E 
R16»  KL=K 
K  16=8E-3 
Q5*K=C-< 
K11=3E-7 
C  =  1  >  1 E  6 
SPZC  OT=0.0 
Pf--I  NT  01  ,02 
PLOT  01=1/0 
PLOT  03=3/0 
RUN 


*N    ="LOOD     MAY     31      18e9 

OF     THE     SOUTH     FORK     DAM 
J+DTM  Rl»  JK-P2,  JK-P3  3  JK) 

7 

l*J*Q2aK)/(K,:*ai»K+l  ) 
-  8 


^Ql  -,  K 

—  6 

*Q1?  K*Q3,  K 

J+0T*<R2.JK+F*, JK-P6 ,JK) 

*Q1 ,K*03jK 

-3 

*Q3o  K 

-"«Q1  ,K 


JOHNSTOWN     1389 


=*Q4-.K*Q5»K*Q5oK*a7>K 

-  39 

1  0*0  6oK 

-5 

J+r)-r-*(  _s  12*  JK) 

12*Q43  K*  05  ->  K  -  06  3  K*  Q"?  ,  K 
~ -3Q 

<+D'*(  RU,JK-R13,JK-R15-)J<-3  16.JK  ) 
2 

1 3*0  43<*Q5,<  *0  6a  K*Q7jK 
E  -3  9 

14*0  *-»K»Q5»K.*Q65  K*0~*o  K 
E— 39 
15-04,  K 
-4 
1  6*03:>  K*Qa.,< 

11"0  4,K 


1  /LENGTH=100/PLT°EP=2/P°"r3E9  =  l  0 

.  0"»  .PI   .R2  .  R3  .P4  ,  R&  ,  01  .  Q5  ,  06.  07  ,  PI  6 

2  =  2/02=3/°3=* 
4=4/  Rl  6=  ♦ 
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EARThCLAKE     SIMULATION 


3/1 4/73 


C 

P 

C 
* 
R 
C 
L 
N 
C 
P 
C 

c 

p 
c 

* 
R 

c 

A 

c 

A 
"L 

PP 

PP 
SP 


EARTHQUAKE     SIMULATION 

CI.K=G1.J+DT=(P1.JK-P2.JK-P3.JK) 

C  l=Gl  I 

CI  I— 2.62E12 

R1.KL=K1*(E1-E2)*C2.K 

K  1=3. SCE-2 
E  1  =  2. 12 
E2-C 

RS .KL=K2*Q1  .K 

K2  =  2 . CAE -2 

C2.k=CE-C1 .K 

CE=3E12 

C3.K=G3  ,J  +  CT*(R4,JK-R7,JK-h;>,JK  ) 

C2  =  0 

54.«L=«4*Q1.<*S.«.K 

*4=.2341 

S7.KL=K7*Q2.K*G4.  < 

K7=7.Q35E-S 

Re«KL=K3*C3.K 
RS«KL=CLIP( KS*C3,K.Q3.K/DT .03 .K/DT • K3*Q3« K) 

.<?=1  .043E6 

C*.K=QA .J  +  CT*<  P9. JK-R1  I . JK-P 10. JK J 

C4  =  Q4  I 

C4 1= .21 *E5 

RQ.KL  =  K,;*E2>«QA.K*G5.< 
£2=.22E5 
kc=1 0.94E-13 

R 1 1 .KL=K I 1*Q4.K 

Kl  1=   .1 

RIO.  KL  =  X10*Q3.K*G4.K 
R  1G.KL=CLIP<K  10«C3.K*Q4,K  t34, x/0TfO4, K/0T.K1 0*G3.K«C4  .<  ) 
K  1  0=  •  3  8  3 
5  .K-  =  CT-Q4.K 
T  =  1E6 

».K  =  PULSE!  10  .  . 9  5E-5  .20  ) 
C 1=1/02= 2/ G3=2/: 4= 4/G5=5 

Gi . C2.G3. G4 ,GS 

Rl.R2.R3.R4.R7.Re.R5 
CT=1  ,<3E-7/LENGTH=3.8E-  5/°LTPER- 1 «9E-7/PRTPER=l • SE-7 


C 
C 

c 

DT  ~ 

INT 

INT 

EC 

N 


2  20 


,koLR-0 I sGSCcft  XubLL 


1  l/_.9/76 


• 

CkC 

* 

jyx 

L 

Ul 

N 

Ql 

v. 

Gl 

rl 

ST 

A 

GT 

c 

CT 

K 

til 

A 

El 

c 

£1 

c 

KL 

r» 

h_: 

A 

c2 

c 

E2 

c 

f.2. 

A 

uZ 

A 

Ql 

A 

GRi 

oP_ 

C    0 

PR! 

NT 

PLL 

iT     G 

KwN 

CK-OISORUCK     MC.uc.i_ 
ANG     Mj__/c_.     PuK     FIG*     2».2 
.  K  =  J  1  .Ji-uI*lK  T. JK-nl  .jK.-Kid.JK) 
=  Gl  I 
I=.0l 

•  KL-=Kl*cl  .K*GT.K 

•  <-=G  T  i 
I  =  l  .0 

tKL=<i»-l  .K.*G  l  .  .< 
.*  =  El  i 
I=i  .0 
=  0.2 

•  Kl_-=  *.-.*_.-.  .K*G  l  .K. 

..<=E2i 

i=l  .  0 
=  .C2 

.  K=GT»N—^  l  »  K 

OCT  •  <=K  T  .on  — r<  l  ._K.—i-<_..JK 

o*TH  ,K-=ui  G«_T  •  K/  Gl  .*. 

T=  .  0  l/LcNt_TH=5G/PL.  TPE?t=l  .0/PkTP=.R=2.5 

GI.G2»i-i_'i_iT»G«G/iTi-»>_T.ril  « R ._ 


l=C  ,Ghu*Th  =  i'i 


,T=T  »k  l  =  l  .R2  =  2  (0  .  I  ) 


MJ^OCATALY^IC    ?PCE~-r :?c: 


-•c 


:l        i  l 

MODEL 


/29/7^ 


-5YVA^!~      V22=L     F"Q?     =!G.     3  .- 
"1  .  K  =0  1  .  wl-f-r  *-:  <  r  ^".  JK-=  I  .  w'K-  =  2  .JO 

a  i = g  i ! 

21  i  =  .  :i 

G  T . K L  =  <  1*GT.K*_1  . ^ 
QT.K=_T  I 

gt:=i  . : 

=  1  .  <  L  =  <  1  *  I  1  .  :<  -  C  1  .  < 

< ; =0.2 

3  2  .<!_  =  <  2  -  2  2  .<  -  C  1  .< 
E2  *<-~2  I 
'2  1=1.3 
K2  =  . 02 

02  .<=";-  .x-gi,< 

.:  i  _  2  ~ .  <  =  =  t  .  j  k  -  -  i  .  j  k  -  ~  2  .  j  < 

3=  :«T-H.  <•=?  12  CT  .  K/G  1  .  K 
C     0"*=.0  i  /."^  2--i=-C/^L*^--=l  .  C /=■=>- -2 
N  -     -1  * 02  *  01  GCr  •  G SO*"*-1 »  i     iS  1  * c2 
JT     01=C.  ;  =  Cl*TH=*.  32=2.^-  =  '. r  1  =  1  .=2  =  2 
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